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Preface 


These lecture notes are intented as a straightforward introduction to partial 
differential equations which can serve as a textbook for undergraduate and 
beginning graduate students. 

For additional reading we recommend following books: W. I. Smirnov [21], 
I. G. Petrowski [17], P. R. Garabedian [8], W. A. Strauss [23], F. John [10], 
L. C. Evans [5] and R. Courant and D. Hilbert[4] and D. Gilbarg and N. S. 
Trudinger [9]. Some material of these lecture notes was taken from some of 
these books. 


Chapter 1 


Introduction 


Ordinary and partial differential equations occur in many applications. An 
ordinary differential equation is a special case of a partial differential equa- 
tion but the behaviour of solutions is quite different in general. It is much 
more complicated in the case of partial differential equations caused by the 
fact that the functions for which we are looking at are functions of more 
than one independent variable. 

Equation 

F(z, y(z),y'(2),-.-,y) =0 

is an ordinary differential equation of n-th order for the unknown function 
y(x), where F' is given. 

An important problem for ordinary differential equations is the initial 
value problem 


y(z) = f(x,y(z)) 
y(zo) = Yo, 


where f is a given real function of two variables x, y and xo, yo are given 
real numbers. 


Picard-Lindel6f Theorem. Suppose 
(i) f(x,y) is continuous in a rectangle 


Q = {(z,y) € R®: |x — 2o| <a, |y— yo| < d}. 


(ii) There is a constant K such that |f(x,y)| < K for all (x,y) €Q. 
(ii) Lipschitz condition: There is a constant L such that 


lf (a, yo) — f(x, y1)| < Liye — y1| 


9 


10 CHAPTER 1. INTRODUCTION 


% 


Xo x 


Figure 1.1: Initial value problem 


for all (x, y1), (@, y2)- 


Then there exists a unique solution y € C!(aj—a, 29 +a) of the above initial 
value problem, where a = min(b/K,a). 


The linear ordinary differential equation 
y™ + an—1(a)y—-) +...a1(2)y’ + ap(x)y = 0, 


where a; are continuous functions, has exactly n linearly independent solu- 
tions. In contrast to this property the partial differential uz, +u,, = 0 in R? 
has infinitely many linearly independent solutions in the linear space C?(R?). 


The ordinary differential equation of second order 


y" (x) — f(z, y(x),y'(#)) 


has in general a family of solutions with two free parameters. Thus, it is 
naturally to consider the associated initial value problem 


y" (x) f(x, y(x), y'(@)) 
y(xo) yo, y (to) = 41, 


where yo and yj; are given, or to consider the boundary value problem 


y" (2) f(z, y(x),y'(2)) 
y(zo) = yo. y(z1) = yL- 


Initial and boundary value problems play an important role also in the 
theory of partial differential equations. A partial differential equation for 
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X90 x) x 
Figure 1.2: Boundary value problem 


the unknown function u(a,y) is for example 
va Y, U, Ug, Uy, Ura, Urys aig) = 0, 


where the function F' is given. This equation is of second order. 

An equation is said to be of n-th order if the highest derivative which 
occurs is of order n. 

An equation is said to be linear if the unknown function and its deriva- 
tives are linear in F’. For example, 


a(z,y)Ue + W(x, y)Uy + c(x, y)u = f(x,y), 


where the functions a, b, c and f are given, is a linear equation of first 
order. 

An equation is said to be quasilinear if it is linear in the highest deriva- 
tives. For example, 


A(X, Ys U, Us, Uy Uae + O(L, Y, U, Uz, Uy )Uay + C(L, Y, U, Ux, Uy)Uyy = 0 


is a quasilinear equation of second order. 


1.1 Examples 


1. uy = 0, where u = u(z,y). All functions u = w(a) are solutions. 


2. Ur = Uy, Where u = u(x,y). A change of coordinates transforms this 
equation into an equation of the first example. Set € =ax+y, n=2-y, 
then 
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Assume u € Cl, then 
1 
Uy 3 (ue — Uy). 
If uz = Uy, then v, = 0 and vice versa, thus v = w(€) are solutions for 
arbitrary C!-functions w(). Consequently, we have a large class of solutions 
of the original partial differential equation: u = w(x + y) with an arbitrary 
C1-function w. 


3. A necessary and sufficient condition such that for given C1-functions 


M, N the integral 
Pi 


M(x, y)dx + N(a,y)dy 
Po 
is independent of the curve which connects the points Py with P; in a simply 
connected domain 2 C R? is the partial differential equation (condition of 
integrability) 


M, = Nz 


in 2). 


x 


Figure 1.3: Independence of the path 


This is one equation for two functions. A large class of solutions is given 
by M = ®,, N = %,, where ®(z,y) is an arbitrary C?-function. It follows 
from Gauss theorem that these are all C!-solutions of the above differential 
equation. 


4. Method of an integrating multiplier for an ordinary differential equation. 
Consider the ordinary differential equation 


M(a,y)da + N(x, y)dy = 0 
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for given C!-functions M, N. Then we seek a C!-function p(x, y) such that 
uMdx + wNdy is a total differential, i. e., that (uM), = (uN)z is satisfied. 
This is a linear partial differential equation of first order for pu: 


Mpy — Npz = L(Ne = M,). 


5. Two C!-functions u(z, y) and u(z, y) are said to be functionally dependent 


if 
act ( eee ) = 0, 
Uz Wy 


which is a linear partial differential equation of first order for u if v is a given 
C'-function. A large class of solutions is given by 


u = H(v(2,y)), 
where H is an arbitrary C!-function. 


6. Cauchy-Riemann equations. Set f(z) = u(x, y)+iu(a, y), where z = x+iy 
and u, v are given C1((Q)-functions. Here is Q a domain in R?. If the function 
f(z) is differentiable with respect to the complex variable z then u, v satisfy 
the Cauchy-Riemann equations 


Ug = Vy, Uy = —Uy. 
It is known from the theory of functions of one complex variable that the 
real part wu and the imaginary part v of a differentiable function f(z) are 
solutions of the Laplace equation 

Au=0, Av=0, 


where Au = tg + Uyy- 


7. The Newton potential 
7 1 


ee ye 


is a solution of the Laplace equation in R® \ (0, 0,0), i. e., of 


Uge + Uyy + Uzz = 0. 
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8. Heat equation. Let u(x,t) be the temperature of a point x € 2 at time 
t, where Q C R? is a domain. Then u(z, t) satisfies in Q x [0,00) the heat 
equation 

Ut = kAu, 


where Au = Ug,2, +Uaroxr9 tUrzr3 and k is a positive constant. The condition 
u(x, 0) as uo(2), ren, 


where uo(x) is given, is an initial condition associated to the above heat 
equation. The condition 


u(x,t) =h(a,t), « € 00, t>0, 


where h(x, t) is given is a boundary condition for the heat equation. 

If h(x, t) = g(x), that is, h is independent of t, then one expects that the 
solution u(x,t) tends to a function u(x) if t > oo. Moreover, it turns out 
that v is the solution of the boundary value problem for the Laplace equation 


Av = 0 inQ 
v = g(x) on OO. 


9. Wave equation. The wave equation 


Figure 1.4: Oscillating string 


Ut = c’Au, 


where u = u(x,t), c is a positive constant, describes oscillations of mem- 
branes or of three dimensional domains, for example. In the one-dimensional 
case 
cme: 
Utt = C Ure 


describes oscillations of a string. 
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Associated initial conditions are 
u(x,0) = uo(x), we(x,0) = u(x), 


where ug, uy are given functions. Thus the initial position and the initial 
velocity are prescribed. 
If the string is finite one describes additionally boundary conditions, for 
example 
u(0,t) =0, u(l,t)=0 for allt > 0. 


1.2 Equations from variational problems 


A large class of ordinary and partial differential equations arise from varia- 
tional problems. 
1.2.1 Ordinary differential equations 
Set 
b 
B(v) = [ fle,v(0),0(a) ae 
and for given ug, up € R 
V ={v € C*[a,b] : v(a) = ua, v(b) = up}, 


where —oo < a < b < ow and f is sufficiently regular. One of the basic 
problems in the calculus of variation is 


(P) mingey E(v). 


Euler equation. Let u € V be a solution of (P), then 


+ fue, u(x), ul(2)) = fulesule),u(2)) 


in (a, 6). 


Proof. Exercise. Hints: For fixed ¢ € C?[a,b] with ¢(a) = ¢(b) = 0 and 
real €, |e] < eo, set g(e) = E(u + ed). Since g(0) < g(e) it follows g’(0) = 0. 
Integration by parts in the formula for g'(0) and the following basic lemma 
in the calculus of variations imply Euler’s equation. 
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Figure 1.5: Admissible variations 


Basic lemma in the calculus of variations. Let h € C(a,b) and 


b 
i h(x)o(a) dx =0 


for all 6 € C§(a,b). Then h(x) = 0 on (a,b). 


Proof. Assume h(ao) > 0 for an xo € (a,b), then there is a 6 > 0 such that 
(xo — 6,20 +0) C (a,b) and h(x) > h(xo)/2 on (ao — 6,20 + 0). Set 


a2) = { (52 — |x — aol)? if « € (to ~ 6,20 +8) 
0 if xe (a,b) \ [zo — 6,20 +4] 


Thus ¢ € Cj(a,b) and 


b x xro+6 
i: h(x) (x) da > (zo) ) (x) dx > 0, 


ae ee 


which is a contradiction to the assumption of the lemma. 


1.2.2 Partial differential equations 


The same procedure as above applied to the following multiple integral leads 
to a second-order quasilinear partial differential equation. Set 


E(v) =/ F(a,v,Vv) da, 
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where 2 C R” is a domain, x = (@1,...,%), vu = v(x) : Qe R, and 
Vv = (vx1,---5Uz,). Assume that the function F is sufficiently regular in 
its arguments. For a given function h, defined on OQ, set 


V ={v €C°(Q): v=h on AQ}. 
Euler equation. Let u € V be a solution of (P), then 


“0 
—F, —F,=0 
» Ox; i 


inQ. 


Proof. Exercise. Hint: Extend the above fundamental lemma of the calculus 
of variations to the case of multiple integrals. The interval (xo — 6, 79 +6) in 
the definition of ¢ must be replaced by a ball with center at x and radius 
é. 


Example: Dirichlet integral 


In two dimensions the Dirichlet integral is given by 
D(v) =} (uz +v,) drdy 
Q 


and the associated Euler equation is the Laplace equation Au = 0 in Q. 
Thus, there is natural relationship between the boundary value problem 


Au=0 inQ, w=h on OQ 
and the variational problem 


min D(v). 

But these problems are not equivalent in general. It can happen that the 
boundary value problem has a solution but the variational problem has no 
solution, see for an example Courant and Hilbert [4], Vol. 1, p. 155, where 
h is a continuous function and the associated solution u of the boundary 
value problem has no finite Dirichlet integral. 

The problems are equivalent, provided the given boundary value function 
h is in the class H'/?(0Q), see Lions and Magenes [14]. 
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Example: Minimal surface equation 


The non-parametric minimal surface problem in two dimensions is to 
find a minimizer u = u(x, 22) of the problem 


min f Tp? sng, day 


where for a given function h defined on the boundary of the domain 2Q 


V={veciQ): v=h on ON}. 


Figure 1.6: Comparison surface 


Suppose that the minimizer satisfies the regularity assumption u € C?(Q), 
then wu is a solution of the minimal surface equation (Euler equation) in Q 


0 Un, 6) Uxe 
t =0. 1.1 
Ox1 (an) 0x2 (as uD 


In fact, the additional assumption u € GA) is superfluous since it follows 
from regularity considerations for quasilinear elliptic equations of second 
order, see for example Gilbarg and Trudinger [9]. 

Let Q = R?. Each linear function is a solution of the minimal surface 
equation (1.1). It was shown by Bernstein [2] that there are no other solu- 
tions of the minimal surface quation. This is true also for higher dimensions 
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n < 7, see Simons [19]. If m > 8, then there exists also other solutions which 
define cones, see Bombieri, De Giorgi and Giusti [3]. 

The linearized minimal surface equation over u = 0 is the Laplace equa- 
tion Au = 0. In R? linear functions are solutions but also many other 
functions in contrast to the minimal surface equation. This striking differ- 
ence is caused by the strong nonlinearity of the minimal surface equation. 

More general minimal surfaces are described by using parametric rep- 
resentations. An example is shown in Figure 1.7!. See [18], pp. 62, for 
example, for rotationally symmetric minimal surfaces. 


Figure 1.7: Rotationally symmetric minimal surface 


Neumann type boundary value problems 


Set V = C1(Q) and 


E(v) =[ F(ax,v,Vv) dx — [. g(x,v) ds, 


where F' and g are given sufficiently regular functions and Q C R” is a 
bounded and sufficiently regular domain. Assume wu is a minimizer of E(v) 
in V, that is 

weV: E(u) < E(v) for allv eV, 


1 An experiment from Beutelspacher’s Mathematikum, Wissenschaftsjahr 2008, Leipzig 
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then 


n 


| (>) Fug, (a, u, Vue; + F,(x,u,Vu)d) dx 
) 


i=1 


— | Gu(x,u)d ds =0 
(oh0) 


for all ¢ € C1(Q). Assume additionally u € C?(Q), then u is a solution of 
the Neumann type boundary value problem 


“0 
y Fu, -F. = 0 ind 
; Ox; * 
i=1 

n 

So Fun, Vi — Ju = 0 onA@Q, 
i=1 


where v = (11,...,%) is the exterior unit normal at the boundary 0Q. This 
follows after integration by parts from the basic lemma of the calculus of 
variations. 


Example: Laplace equation 


Set ‘ 
E(v) = 5 | |Vul|? dx -| h(a)u ds, 
2 Jo an 
then the associated boundary value problem is 
Au = 0 inQ 


Ou 
OV 


= h on OQ. 
Example: Capillary equation 
Let 2 C R? and set 


BE) = [ VIFTWeP ae + 5 f 2 de— cosy f vu ds. 
Q Q an 


Here « is a positive constant (capillarity constant) and + is the (constant) 
boundary contact angle, i. e., the angle between the container wall and 
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the capillary surface, defined by v = v(x, 22), at the boundary. Then the 
related boundary value problem is 


div (Tu) = «cu inQ 
y-Tu = cosy on ag, 


where we use the abbreviation 
= Vu 


OT + Fue 


div (Tu) is the left hand side of the minimal surface equation (1.1) and it 
is twice the mean curvature of the surface defined by z = u(21, 22), see an 
exercise. 

The above problem describes the ascent of a liquid, water for example, 
in a vertical cylinder with cross section 2. Assume the gravity is directed 
downwards in the direction of the negative x3-axis. Figure 1.8 shows that 
liquid can rise along a vertical wedge which is a consequence of the strong 
nonlinearity of the underlying equations, see Finn [7]. This photo was taken 


Figure 1.8: Ascent of liquid in a wedge 


from [15]. 
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1.3. Exercises 


1. 


Find nontrivial solutions u of 


Ugy — Uyx =0. 


. Prove: In the linear space C?(R?) there are infinitely many linearly 


independent solutions of Au = 0 in R?. 


Hint: Real and imaginary part of holomorphic functions are solutions 
of the Laplace equation. 


. Find all radially symmetric functions which satisfy the Laplace equa- 


tion in R”\ {0} for n > 2. A function wu is said to be radially symmetric 
if u(x) = f(r), where r = (Sc? 2?)/2. 


Hint: Show that a radially symmetric u satisfies Au = r'~” (ere): 
by using Vu(x) = f’(r)#. 


r 


. Prove the basic lemma in the calculus of variations: Let Q Cc R” bea 


domain and f € C(Q) such that 
if f(a)h(a) dz =0 
Q 


for all h € C2(Q). Then f =0inQ. 


. Write the minimal surface equation (1.1) as a quasilinear equation of 


second order. 


. Prove that a sufficiently regular minimizer in C!(Q) of 


E(v) > F(x,v,Vv) dx -| g(v,v) ds, 
Q 0Q 
is a solution of the boundary value problem 


e -F = 0mQ 


n 
Bit =o = 0 onan, 
i=1 


where v = (11,...,,) is the exterior unit normal at the boundary 0. 


1.3. 


10. 


EXERCISES 23 


. Prove that vy - Tu = cosy on OQ, where 7¥ is the angle between the 


container wall, which is here a cylinder, and the surface S, defined by 
z= u(#1, £2), at the boundary of S, v is the exterior normal at 00. 


Hint: The angle between two surfaces is by definition the angle between 
the two associated normals at the intersection of the surfaces. 


. Let 2 be bounded and assume u € C?(Q) is a solution of 
divTu = CinQ 
V- ons EN = cosy on OQ, 


JV1+ |Vul? 
where C’ is a constant. 
Prove that 


[AQ 
C = — cosy. 
[2| 


Hint: Integrate the differential equation over 2. 


. Assume 2 = Br(0) is a disc with radius R and the center at the origin. 


Show that radially symmetric solutions u(x) = w(r), r = Vx} + 23, 
of the capillary boundary value problem are solutions of 


( rw i 
—— = Krw inOd<r<R 
===] 


w! 


Vl+w? 


cosy ifr=R. 


Remark. It follows from a maximum principle of Concus and Finn [7] 
that a solution of the capillary equation over a disc must be radially 
symmetric. 


Find all radially symmetric solutions of 


( rw : 
—— = Cr in0O<r<R 
as 


w' 


Vvl+w? 


Hint: From an exercise above it follows that 


cosy ifr=R. 


2 
C= Ros? 
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11. Show that div Tu is twice the mean curvature of the surface defined 
by z= u(x, 22). 


Chapter 2 


Equations of first order 


For a given sufficiently regular function F’ the general equation of first order 
for the unknown function u(z) is 


F(x,u, Vu) =0 


in Q € R”. The main tool for studying related problems is the theory of 
ordinary differential equations. This is quite different for systems of partial 
differential of first order. 

The general linear partial differential equation of first order can be writ- 


ten as 
n 


Yo ai(x)ux, + e(x)u = f(a) 

i=1 
for given functions a;, c and f. The general quasilinear partial differential 
equation of first order is 


n 


SS ai(X, u)ue, + c(x, u) = 0. 


i=1 


2.1 Linear equations 


Let us begin with the linear homogeneous equation 
ai (2; y)te + aoe; y)uy = 0. (2.1) 


Assume there is a C!-solution z = u(x,y). This function defines a surface 
S which has at P = (2, y, u(x, y)) the normal 
1 
N= (=tizy ty; 1) 


J1+ |Vul? 
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and the tangential plane defined by 


C— 2 = Ura, y)( — &) + Uy(a,y)(n — y): 


Set p = uy(z,y), g = uy(x,y) and z = u(x,y). The tuple (x,y, z,p,¢q) is 
called surface element and the tuple (x, y, z) support of the surface element. 
The tangential plane is defined by the surface element. On the other hand, 
differential equation (2.1) 


ai(x, y)p + a2(x,y)q =0 


defines at each support (x, y, z) a bundle of planes if we consider all (p, q) sat- 
isfying this equation. For fixed (x, y), this family of planes II(A) = II(,; 2, y) 
is defined by a one parameter family of ascents p(A) = p(\;2,y), q(A) = 
q(A; x,y). The envelope of these planes is a line since 


a(x, y)p(A) + a2(x, y)q(A) = 0, 


which implies that the normal N(A) on II(A) is perpendicular on (a1, ae, 0). 
Consider a curve x(T) = (x(7T), y(7T), 2(7)) on S, let Tx, be the tangential 
plane at xo = ((70), y(70), 2(70)) of S and consider on Ty, the line 


L: l(c) =xo+ox'(m), ER, 


see Figure 2.1. 

We assume LF coincides with the envelope, which is a line here, of the 
family of planes II(A) at (x,y,z). Assume that T;, = II(Ao) and consider 
two planes 


(Ao): 2-2 = (@—20)p(Ao) + (¥ — yo)a(o) 
H(Ao +h): 2-2 = (&-20)p(ro +h) + (y— yo)a(ro + h)- 


At the intersection I(a) we have 
(2 — xo)p(Ao) + (y = yo)a\Ao) = (@ — @o)p(Ao + h) + (y = yo)aXo + Fh). 


Thus, 
a’ (t0)p' (Ao) + y' (70) q' (Ao) = 9. 


From the differential equation 


a1(2(70), ¥(70) PA) + @2(@(70); y(70) a) = 0 
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xX 


Figure 2.1: Curve on a surface 


it follows 
ap’ (Ao) + a2q'(Ao) = 0. 


Consequently 
x'(r) 


(a'(r), y'(7)) = Ga a) 


(a1(x(7), y(7)), a2(2(7), y(7)), 


since 7) was an arbitrary parameter. Here we assume that x’(r) 4 0 and 


ai(2(T), y(T)) #0. 


Then we introduce a new parameter t by the inverse of 7 = T(t), where 


t(T) = i: a) ds. 


» M(x(s), y(s)) 


It follows x(t) = ai(x,y), y’(t) = ao(az, y). We denote x(r(t)) by x(t) again. 
Now we consider the initial value problem 


a’ (t) = ai(x,y), y'(t) = a(x, y), x(0) = 2X0, y(0) = Vo- (2.2) 


From the theory of ordinary differential equations it follows (Theorem of 
Picard-Lindelof) that there is a unique solution in a neighbouhood of t = 0 
provided the functions a1, a2 are in C!. From this definition of the curves 
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(x(t), y(t)) is follows that the field of directions (a1 (20, yo), @2(Xo, yo)) defines 
the slope of these curves at (x(0), y(0)). 


Definition. The differential equations in (2.2) are called characteristic 
equations or characteristic system and solutions of the associated initial value 
problem are called characteristic curves. 


Definition. A function ¢(, y) is said to be an integral of the characteristic 


system if o(x(t), y(t)) = const. for each characteristic curve. The constant 
depends on the characteristic curve considered. 


Proposition 2.1. Assume ¢ € C! is an integral, then u = $(a,y) is a 
solution of (2.1). 


Proof. Consider for given (29, yo) the above initial value problem (2.2). 
Since $(x(t), y(t)) = const. it follows 


yx! st dyy’ =0 
for |t| < to, to > 0 and sufficiently small. Thus 


bx (Xo, Yo)a1(Z0, Yo) + Sy(Xo, Yo)a2(ro, Yo) = 0. 


Remark. If ¢(z,y) is a solution of equation (2.1) then also H(¢(z, y)), 
where H(s) is a given C!-function. 


Examples 


1. Consider 
aUz + aQUy = 0, 


where a1, Gg are constants. The system of characteristic equations is 
a! = 41, y = a2. 
Thus the characteristic curves are parallel straight lines defined by 


x=ayt+A, y=agt+B, 
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where A, B are arbitrary constants. From these equations it follows that 


P(@, y) = agr — ary 


is constant along each characteristic curve. Consequently, see Proposi- 
tion 2.1, wu = agx — ayy is a solution of the differential equation. From 
an exercise it follows that 


u = H(agx — ary), (2.3) 


where H(s) is an arbitrary C1-function, is also a solution. Since u is constant 
when aga — a,y is constant, equation (2.3) defines cylinder surfaces which 
are generated by parallel straight lines which are parallel to the (x, y)-plane, 
see Figure 2.2. 


Z 


eer 


Xx 


Figure 2.2: Cylinder surfaces 


2. Consider the differential equation 
LUg + yuy = 0. 
The characteristic equations are 


We a 
T=2X,Y =Y; 
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and the characteristic curves are given by 
x= Ae’, y= Be’, 


where A, B are arbitrary constants. Thus, an integral is y/x, x 4 0, and for 
a given C-function the function u = H(2/y) is a solution of the differential 
equation. If y/x = const., then u is constant. Suppose that H’(s) > 0, 
for example, then u defines right helicoids (in German: Wendelflachen), see 
Figure 2.3 


Figure 2.3: Right helicoid, a2 < 2? + y? < R? (Museo Ideale Leonardo da 
Vinci, Italy) 


3. Consider the differential equation 
YUg — Tuy = 0. 


The associated characteristic system is 


If follows 
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or, equivalently, 

d 

qe +y") =0, 

which implies that «2? + y? = const. along each characteristic. Thus, rota- 
tionally symmetric surfaces defined by u = H(a? + y), where H' 0, are 
solutions of the differential equation. 


4. The associated characteristic equations to 
ayUg + bruy = 0, 

where a, 0 are positive constants, are given by 
x’ =ay, y = be. 


It follows baa’ — ayy’ = 0, or equivalently, 


d 
qn — ay’) =0. 


Solutions of the differential equation are u = H(bax? — ay”), which define 
surfaces which have a hyperbola as the intersection with planes parallel to 
the (x, y)-plane. Here H(s) is an arbitrary C'-function, H’(s) 4 0. 


2.2 Quasilinear equations 


Here we consider the equation 
ai (X,Y, UUs + a2(X,Yy, U)Uy = a3(x, y, U). (2.4) 
The inhomogeneous linear equation 
ai (@, YU + a2(a, y)Uy = a3(x,y) 


is a special case of (2.4). 


One arrives at characteristic equations 7’ = a1, y’ = ao, 2’ = a3 
from (2.4) by the same arguments as in the case of homogeneous linear 
equations in two variables. The additional equation z’ = a3 follows from 


Z(t) = piaja(r) + aQ)y'(7) 
= pa, + qa2 

= 43, 

see also Section 2.3, where the general case of nonlinear equations in two 

variables is considered. 
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2.2.1 A linearization method 


We can transform the inhomogeneous equation (2.4) into a homogeneous 
linear equation for an unknown function of three variables by the following 
trick. 

We are looking for a function w(z, y, u) such that the solution u = u(a, y) 
of (2.4) is defined implicitly by ~(az, y,u) = const. Assume there is such a 
function 7 and let u be a solution of (2.4), then 


Vx + Wut, = 0, Wy + WyUy = 0. 
Assume W, 4 0, then 


ig Pe t= by 
tu? 9 the 
From (2.4) we obtain 
ay (2, Y; Z) Wy a ao(x, Yy; Z) py aS a3(x, y; Z)y = 0, (2.5) 


where z := u. 
We consider the associated system of characteristic equations 


tly == ae ase) 
y(t) = a2(x, y, Z) 
oi) = aa, weh 


One arrives at this system by the same arguments as in the two-dimensional 
case above. 


Proposition 2.2. (i) Assume w € C!, w = w(2,y,z), is an integral, i. 
e., it is constant along each fired solution of (2.5), then pW = w(a,y, z) is a 
solution of (2.5). 

(ii) The function z = u(x, y), implicitly defined through w(x, u, z) = const., 
is a solution of (2.4), provided that w, 4 0. 

(iii) Let z = u(x, y) be a solution of (2.4) and let (x(t), y(t)) be a solution of 


a’ (t) = ay(z,y, u(x, y)), y'(t) = a2(x,y, u(x, y)), 


then z(t) := u(x(t), y(t)) satisfies the third of the above characteristic equa- 
tions. 


Proof. Exercise. 
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2.2.2 Initial value problem of Cauchy 


Consider again the quasilinear equation 
(x) ay (2, Y, U)Ux + a(x, y, U)Uy = a3(x, y, U). 
Let 
T: x=20(s), y=yols), z= 20(s), $1 < 8 < 82, —CO < 81 < 82 < +00 


be a regular curve in R? and denote by C the orthogonal projection of I 
onto the (x, y)-plane, i. e., 


Initial value problem of Cauchy: Find a C!-solution u = u(x,y) of 
(x) such that u(xo(s), yo(s)) = zo(s), % €., we seek a surface S defined by 
z=u(x,y) which contains the curve T. 


Z 


Figure 2.4: Cauchy initial value problem 


Definition. The curve - is said to be noncharacteristic if 


Xo(8)a2(xo(s), ¥o(s)) — yo(s)a1(w0(s), yols)) F 0. 


Theorem 2.1. Assume aj, a2, a2 € C! in their arguments, the initial data 
x0, Yo, 20 € C![s1, 82] andT is noncharacteristic. 
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Then there is a neighbourhood of C such that there exists exactly one 
solution wu of the Cauchy initial value problem. 


Proof. (i) Existence. Consider the following initial value problem for the 
system of characteristic equations to (x): 


w(t) = "a (29,2) 
y(t) = a2(x, Y; z) 
z'(t) = a3 (a, Y; z) 


with the initial conditions 


y(s,0) = yo(s) 
2(s,0) = z0(s). 


Let x = x(s,t), y = y(s,t), z = z(s,t) be the solution, s; < s < sg, |t] <7 
for an 7 > 0. We will show that this set of strings sticked onto the curve 
I, see Figure 2.4, defines a surface. To show this, we consider the inverse 
functions s = s(x,y), t = t(z,y) of x = x(s,t), y = y(s,t) and show that 
z(s(x,y), t(, y)) is a solution of the initial problem of Cauchy. The inverse 
functions s and ¢ exist in a neighbourhood of t = 0 since 


Ca 
* 8G,1) La 


Ls Lt 
Ys Yt 


— r(s)az — Yo(s)a1 #0, 
t=0 


and the initial curve [ is noncharacteristic by assumption. 
Set 


u(x, y) = 2(s(z,y), H(z, y)), 
then wu satisfies the initial condition since 
u(x, y)|t=0 = 2(s,0) = 20(s). 


The following calculation shows that wu is also a solution of the differential 
equation (x). 


G1, + aotty = 1(2sSe + Zitz) + d2(ZsSy + Zty) 
= 23(015q + G2Sy) + 2(aite + aaty) 
= 25(Sr%t + Syyt) + 24(tete + tyYe) 


= a3 
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since 0 = s; = 8,24 + Syy and 1 = t; = tra, + tyys. 


(ii) Uniqueness. Suppose that v(x, y) is a second solution. Consider a point 
(z’,y’) in a neighbourhood of the curve (xo(s), y(s)), $1 —€ < 8 < sg +6, 
€ > 0 small. The inverse parameters are s’ = s(2’,y’), t! = t(a’,y’), see 
Figure 2.5. 


(x’,y’) 


(X9(s’),Y (8’) 


Figure 2.5: Uniqueness proof 


Let 
Ar 2) Hale tye t) 20) = 205, B) 
be the solution of the above initial value problem for the characteristic dif- 
ferential equations with the initial data 


x(s’,0) = x0(s’), y(s’,0) = yo(s’), z(s’,0) = 20(s’). 


According to its construction this curve is on the surface S defined by u = 
u(a,y) and u(a’,y’) = z(s’,t’/). Set 


V(t) = v(x(t), y(t)) — 2), 
then 
W(t) = vga! + vyy! — 2! 
= £1 + Vya2 — a3 = 0 


and 

w(0) = v(a(s', 0), y(s’,0)) — z(s’,0) = 0 
since v is a solution of the differential equation and satisfies the initial con- 
dition by assumption. Thus, 2(t) = 0, i. e., 


v(ax(s’,t), y(s’,t)) — z(s’,t) =0. 
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Set t = 7’, then 
u(a’,y’) — z(s’, t') = 0, 


which shows that v(2’,y’) = u(2’, y’) because of z(s’, t’) = u(a’, y’). 


Remark. In general, there is no uniqueness if the initial curve T is a 
characteristic curve, see an exercise and Figure 2.6 which illustrates this 
case. 


Xx 


Figure 2.6: Multiple solutions 


Examples 
1. Consider the Cauchy initial value problem 
Ug + Uy = 0 
with the initial data 
ro(s) = 8, yo(s) =1, z0(s) is a given C'-function. 


These initial data are noncharacteristic since yja1— «paz = —1. The solution 
of the associated system of characteristic equations 
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with the initial conditions 


x(s,0) = xo(s), y(s,0) = yo(s), z(s,0) = zo(s) 


is given by 
x=t4+20(s), y=t+yo(s), z = 2(s), 


i. €., 


e=t+s, y=tt+1, z= 2%(s). 


It follows s = x—y+1, t=y-—1 and that u = 20(4 — y +1) is the solution 
of the Cauchy initial value problem. 


2. A problem from kinetics in chemistry. Consider for x > 0, y > 0 the 
problem 


Ug + Uy = (koe + kp) (1— wu) 
with initial data 
u(x,0) =0, x > 0, and u(0,y) = uo(y), y > 0. 


Here the constants k; are positive, these constants define the velocity of the 
reactions in consideration, and the function uo(y) is given. The variable x 
is the time and y is the hight of a tube, for example, in which the chemical 
reaction takes place, and u is the concentration of the chemical substance. 
In contrast to our previous assumptions, the initial data are not in Ct. 


The projection C; UC, of the initial curve onto the (x, y)-plane has a corner 
at the origin, see Figure 2.7. 
y 
X=y 
2, 
CG 
OQ) 
C x 


Figure 2.7: Domains to the chemical kinetics example 
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The associated system of characteristic equations is 
o'(t) =1, ¥(®) =1, 2) = (hoe + ke) (12). 


It follows x =t+c,, y = t+ cy, with constants c;. Thus the projection 
of the characteristic curves on the (x, y)-plane are straight lines parallel to 
y = x. We will solve the initial value problems in the domains ; and Qo, 
see Figure 2.7, separately. 


(i) The initial value problem in Qy. The initial data are 
xo(s) = 8, yo(s) =0, zo(0) =0, s > 0. 


It follows 
r=x2(s,t)=t+s, y=y(s,t) =t. 


Thus 
z'(t) = (koe ™ +9) + ko)(1 — z), 2(0) = 0. 


The solution of this initial value problem is given by 


ki ki 
z(s,t) =1— exp (Ferber Cf Beh) 


Consequently 


k 
Ul (x, y) = 1- exp (2 e hie kay koko» J 


is the solution of the Cauchy initial value problem in 9). If time x tends to 


co, we get the limit 


lim u1(2,y) = 1—e7*2. 
«w—-0oO 


(ii) The initial value problem in Qe. The initial data are here 


xo(s) = 0, yo(s) = 8, z0(0) = uo(s), s > 0. 


It follows 


c= a2(s,t)=t, y=y(s,t)=tts. 


Thus 
z'(t) = (koe ™* + kg)(1 — z), z(0) =O. 
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The solution of this initial value problem is given by 


k k, 
z(s,t) =1-— (1— uo(s)) exp (Fem kot 2) : 


Consequently 


k; 
ua(es9) = 1 (= aly ~2)) exp (Fe-® — hye ~ 2 ) 


is the solution in Qo. 
If ¢ = y, then 


Ki k; 
wiley) = 1 exp (Fert koa: 2) 


RO jus k 
ue(x,y) = 1 = (1=ua(0))exp (Fe RU baa 2). 


If uo(0) > 0, then ws < wz if x = y, i. e., there is a jump of the concentration 
of the substrate along its burning front defined by x = y. 


Remark. Such a problem with discontinuous initial data is called Riemann 
problem. See an exercise for another Riemann problem. 
The case that a solution of the equation is known 


Here we will see that we get immediately a solution of the Cauchy initial 
value problem if a solution of the homogeneous linear equation 


ay (@, y)Uz + ao(z, y)Uy = 0 


is known. 
Let 
r0(s), yo(s), z0(s), 81 <8 < 82 
be the initial data and let u = ¢(2, y) be a solution of the differential equa- 
tion. We assume that 
x(0(8), Yo(s))2o(8) + dy(xo0(s), Yo(s))¥o(s) FO 
is satisfied. Set g(s) = $(xo0(s), yo(s)) and let s = h(g) be the inverse 


function. 


The solution of the Cauchy initial problem is given by uo (h( (2, y))). 
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This follows since in the problem considered a composition of a solution is 
a solution again, see an exercise, and since 


uo (h(P(#0(s); Yo(s))) = uo(h(g)) = uo(s). 


Example: Consider equation 
Ug + Uy = 0 
with initial data 
xo(s) = 8, yo(s) =1, uo(s) is a given function. 
A solution of the differential equation is $(x, y) = x — y. Thus 
$((%0(s), yo(s)) =s—1 


and 
ug(d +1) = u(x —y + 1) 


is the solution of the problem. 


2.3. Nonlinear equations in two variables 
Here we consider equation 


Pe, 2,2, q) = 0, (2.6) 


where z = u(x,y), p = Ux(2,y), Y = Uy(z, y) and F € C? is given such that 
FR +B? #0. 

In contrast to the quasilinear case, this general nonlinear equation is 
more complicated. Together with (2.6) we will consider the following system 
of ordinary equations which follow from considerations below as necessary 
conditions, in particular from the assumption that there is a solution of 


(2.6). 


a(t) = F, (2.7) 
y(t) = Fy (2.8) 
Z(t) = phyt+qFy (2.9) 
p(t) = —Fr- Fup (2.10) 
q(t) = —F,—Fuq (2.11) 
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Definition. Equations (2.7)—(2.11) are said to be characteristic equations 
of equation (2.6) and a solution 


(x(t), yd), 2(¢), p(t), a(t)) 


of the characteristic equations is called characteristic strip or Monge curve. 


Figure 2.8: Gaspard Monge (Panthéon, Paris) 


We will see, as in the quasilinear case, that the strips defined by the char- 
acteristic equations build the solution surface of the Cauchy initial value 
problem. 

Let z = u(x, y) be a solution of the general nonlinear differential equa- 
tion (2.6). 

Let (0, yo, Zo) be fixed, then equation (2.6) defines a set of planes given 
by (0, Yo, 20,2, q); i. €., planes given by z = v(x, y) which contain the point 
(x0, Yo, 29) and for which v, = p, vy = q at (xo, yo). In the case of quasilinear 
equations these set of planes is a bundle of planes which all contain a fixed 
straight line, see Section 2.1. In the general case of this section the situation 
is more complicated. 


Consider the example 


py + ¢ = FCs), (2.12) 
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where f is a given positive function. Let E be a plane defined by z = v(x, y) 
and which contains (9, yo, Zo). Then the normal on the plane F directed 
downward is 


1 
N= ? 9 iL > 
Tae” al) 


where p = Uz(X0, Yo); ¢ = Vy(Xo, yo). It follows from (2.12) that the normal 
N makes a constant angle with the z-axis, and the z-coordinate of N is 
constant, see Figure 2.9. 


(a) 


xX 


Figure 2.9: Monge cone in an example 


Thus the endpoints of the normals fixed at (xo, yo, 20) define a circle 
parallel to the (x, y)-plane, i. e., there is a cone which is the envelope of all 
these planes. 


We assume that the general equation (2.6) defines such a Monge cone at 
each point in R°. Then we seek a surface S which touches at each point its 
Monge cone, see Figure 2.10. 

More precisely, we assume there exists, as in the above example, a one 
parameter C!-family 


P(A) = PAS #952), WA) = Aix, y, 2) 
of solutions of (2.6). These (p(A), q(A)) define a family II(A) of planes. 
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BD 
DoS 


ra 


» 


Figure 2.10: Monge cones 


Let 
x(7) = (x(r), y(7), 2(7)) 
be a curve on the surface S which touches at each point its Monge cone, 
see Figure 2.11. Thus we assume that at each point of the surface S the 
associated tangent plane coincides with a plane from the family II(A) at 


this point. Consider the tangential plane T;, of the surface S at x9 = 
(x(70), y(70), 2(70)). The straight line 


l(c) = xo + 0x’ (70), —0<a<o, 


is an apothem (in German: Mantellinie) of the cone by assumption and is 
contained in the tangential plane T;,, as the tangent of a curve on the surface 
S. It is defined through 

x'(7|) =V(o). (2.13) 


The straight line l(c) satisfies 
I3(a) — 20 = (l(a) — 0)p(Ao) + (l2(7) — Yo)aQo); 


since it is contained in the tangential plane T;, defined by the slope (p, q). 
It follows 


13(0) = p(Ao)i (a) + q(Ao)lo(a). 


44 CHAPTER 2. EQUATIONS OF FIRST ORDER 


Xx 


Figure 2.11: Monge cones along a curve on the surface 


Together with (2.13) we obtain 
z'(7) = p(Ao)2" (rT) + q(Ao)y'(7). (2.14) 


The above straight line 1 is the limit of the intersection line of two neigh- 
bouring planes which envelopes the Monge cone: 


z— 2% = (x—20)p(Ao) + (y — yo)a(ro) 
z— 2 = (x-—29)p(Ao +h) + (y— yo)a(ro +h). 


On the intersection one has 
(x — x0)p(A) + (y — yo)a(Ao) = (a — to)p(Ao + h) + (y — yo)a(Ao + fh). 
Let h — 0, it follows 
(x — xo)p"(Ao) + (y — yo)q'(Ao) = 9. 
Since « =1,(c), y = lo() in this limit position, we have 
P' (Ao) (a) + g'(Ao)lo(a) = 0, 
and it follows from (2.13) that 


p'(Ao)a"(r) + d'(Ao)y'(7) = 0. (2.15) 
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From the differential equation F(x, yo, 20, P(A), g(A)) = 0 we see that 
Fyp' (A) + Fqq'(A) = 0. (2.16) 


Assume x'(7) # 0 and F, £4 0, then we obtain from (2.15), (2.16) 


ull) _ Fy 
a'(7>) Fp’ 
and from (2.14) (2.16) that 
Zl (m F 
s a) =prg - 
x! (79) Fy 


It follows, since 7 was an arbitrary fixed parameter, 


x(7) = (e'(r).v().20) 
o'(r).a'(r) Fa) (+ 054) ) 
a(n) 


= (Fp, Fq, pF p + qF a), 
Pp 


i. e., the tangential vector x’(r) is proportional to (Fp, Fy, pF» + qFy). Set 


where F = F(x(r), y(t), 2(7), p(A(7)), ¢(A(7))).. Introducing the new pa- 
rameter ¢ by the inverse of r = r(t), where 


we obtain the characteristic equations (2.7)—(2.9). Here we denote x(r7(t)) 
by x(t) again. From the differential equation (2.6) and from (2.7)-(2.9) 
we get equations (2.10) and (2.11). Assume the surface z = u(x, y) under 
consideration is in C?, then 


Fy + Fp + pps + Fopy 
Ff, + Fept+ x(t) py a y'(t)py 
FF, + Fp +p (t) 


) (dx = Py) 


oOo O&O 
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since p = p(x, y) = p(x(t), y(t)) on the curve x(t). Thus equation (2.10) of 
the characteristic system is shown. Differentiating the differential equation 
(2.6) with respect to y, we get finally equation (2.11). 


Remark. In the previous quasilinear case 


Fe; Y,2,P, q) = ai(2, Y, z)p + a2(a, Yy, z)q =~ a3(2, Y, z) 


the first three characteristic equations are the same: 


a’ (t) = ai(z,y, 2), y'(t) — a(x, y, z), z(t) = a3(x, y, z). 


The point is that the right hand sides are independent on p or q. It follows 
from Theorem 2.1 that there exists a solution of the Cauchy initial value 
problem provided the initial data are noncharacteristic. That is, we do not 
need the other remaining two characteristic equations. 

The other two equations (2.10) and (2.11) are satisfied in this quasilin- 
ear case automatically if there is a solution of the equation, see the above 
derivation of these equations. 


The geometric meaning of the first three characteristic differential equa- 
tions (2.7)-(2.11) is the following one. Each point of the curve 
A: (a(t), y(t), z(t)) corresponds a tangential plane with the normal direc- 
tion (—p, —q, 1) such that 


2'(t) = p(t)2'(t) + a(d)y'(t). 


This equation is called strip condition. On the other hand, let z = u(z,y) 
defines a surface, then z(t) := u(«(t), y(t)) satisfies the strip condition, where 
p= Uy and q = uy, that is, the ’scales” defined by the normals fit together. 


Proposition 2.3. F(z,y,z,p,q) is an integral, i. e., it is constant along 
each characteristic curve. 


Proof. 


d 

ar (at) wt), 2(t), pt), a(t) = Fee Fyy Bie Pep Fyd 
= Fy, Fy t+ Fyfg + pF Fp + oF Fy 

—Fy fe — FFp — Fyky — Fy F4 


0. 
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Corollary. Assume F(x, yo, 20, Po, go) = 0, then F’ = 0 along characteristic 
curves with the initial data (x0, yo, 20, Po; qo): 


Proposition 2.4. Let z = u(x,y), u € C?, be a solution of the nonlinear 
equation (2.6). Set 


20 = u(xo, Yo; ) Po= Ux (Xo, Yo), do = Uy (Xo, Yo): 


Then the associated characteristic strip is in the surface S, defined by z = 
u(a,y). Thus 


2(t) = u(a(t),y(4)) 
p(t) = ur(a(t), y(t)) 
q(t) = uy(2(t), y@), 


where (a(t), y(t), z(t), p(t), q(t)) is the solution of the characteristic system 
(2.7)-(2.11) with initial data (xo, yo, 20; Po, qo) 


Proof. Consider the initial value problem 
v(t) = F,(2,y,u(x,y), Ur(z,y), Uy(x,y)) 
y(t) Fa(x, y, u(2, y), Ue, y), Uy(x, y)) 
with the initial data (0) = xo, y(0) = yo. We will show that 
(x(t), y(t), w(a(t), y(t)), Ua (w(t), y(t), Uy (a(t), y(t))) 


is a solution of the characteristic system. We recall that the solution exists 
and is uniquely determined. 


Set z(t) = u(x(t), y(t)), then (x(t), y(t), z(t)) C S, and 
z'(t) = ua" (t) + uyy'(t) = UrFp + Uy Fy. 
Set p(t) = ux(x(t), y(t)), a(t) = uy (x(t), y(t), then 
P(t) = UseFy + UryFy 
C(t) = Uye Fy + Uyy Fy. 


Finally, from the differential equation F(x, y, u(x, y), Ue(2, y), Uy(x, y)) = 0 
it follows 

p(t) = -F,-Fup 

q(t) = —F,- Fug. 
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2.3.1 Initial value problem of Cauchy 
Let 


z= x9(s), y= yo(s), z= 20(s), p= pols), g=490(s), $1 < 8 < $2, (2.17) 


be a given initial strip such that the strip condition 


2(8) = po(s)x9(s) + qo(s)¥o(s) (2.18) 


is satisfied. Moreover, we assume that the initial strip satisfies the nonlinear 
equation, that is, 


F'(xo(s), yo(s), 20(8), Po(s), go(s)) = 0. (2.19) 


Initial value problem of Cauchy: Find a C?-solution z = u(a,y) of 
F(x,y,2,p,q) = 0 such that the surface S defined by z = u(x,y) contains 
the above initial strip. 


Similar to the quasilinear case we will show that the set of strips de- 
fined by the characteristic system which are sticked at the initial strip, see 
Figure 2.12, fit together and define the surface for which we are looking at. 


Definition. A strip (2(7), y(7), z(7), p(T), g(7)), T1 < T < 72, is said to be 
noncharacteristic if 


w'(r)Fa(a(r), (7), 2(7),P(7), (7) -¥! (7) Fp(x(7), y(); 2(7),P(7), a7) F 0. 


Theorem 2.2. For a given noncharacteristic initial strip (2.17), x0, Yo, Zo € 

C? and po, qo € C! which satisfies the strip condition (2.18) and the dif- 
ferential equation (2.19) there exists exactly one solution z = u(x,y) of 

the Cauchy initial value problem in a neighbourhood of the initial curve 
(xo(s), yo(s), 20(s)), «. @, z= u(x, y) is the solution of the differential equa- 

tion (2.6) and u(xo(s), yo(s)) = 20(8), Ue(ao(s), yo(s)) = pols), Uy(wo(s), yo(s)) = 
qo(s). 


Proof. Consider the system (2.7)—(2.11) with initial data 
z(s,0) = x9(s), y(s, 0) = yo(s), z(s,0) = z0(s), p(s, 0) = po(s), q(s, 0) = qo(s). 


We will show that the surface defined by x = x(s,t), y(s,t) is the surface 
defined by z = u(x,y), where wu is the solution of the Cauchy initial value 
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x 
Figure 2.12: Construction of the solution 
problem. It turns out that u(x, y) = z(s(z,y),t(x,y)), where s = s(x, y), 


t = t(x, y) is the inverse of x = x(s,t), y = y(s,t) in a neigbourhood of t = 0. 
This inverse exists since the initial strip is noncharacteristic by assumption: 


(x,y) 
ost O(s,t) lt=0 


= Looky — yorg x 0. 


Set 
P(x,y) = p(s(a,y), t(z,y)), Q(x, y) = a(s(a, y), H(z, y)). 
From Proposition 2.3 and Proposition 2.4 it follows F(x, y,u, P,Q) =0. We 


will show that P(x,y) = uz(x,y) and Q(z,y) = uy(x,y). To see this, we 
consider the function 


h(s,t) = Zs — P@s — Ws. 
One has 
h(s,0) = 29(s) — po(s)2°o(s) — qo(s)yo(s) = 0 
since the initial strip satisfies the strip condition by assumption. In the 


following we will find that for fixed s the function h satisfies a linear ho- 
mogeneous ordininary differential equation of first order. Consequently, 
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h(s,t) = 0 in a neighbourhood of t = 0. Thus the strip condition is also sat- 
isfied along strips transversally to the characteristic strips, see Figure 2.18. 
Thaen the set of ’scales” fit together and define a surface like the scales of 
a fish. 


From the definition of h(s,t) and the characteristic equations we get 


hi(s,t) 


Zst — Pts — GtYs — PXst — CYst 


0 
= 75 (% PLE — Gye) + Pate + Gai — UYs — Pits 


= (pr, + qys) Fz + Futs + Fyzs + Fops + Fos. 


Since F'(x(s,t), y(s, t), 2(s,t), p(s, t), ¢(s,t)) = 0, it follows after differentia- 
tion of this equation with respect to s the differential equation 


hy = —F-Yh. 
Hence h(s,t) = 0, since h(s,0) = 0. 


Thus we have 


Zs = PXs 7 Ws 
St = Prev Wt 
Zs = Ugls + UyYs 
Ze = Unt + UyVe- 


The first equation was shown above, the second is a characteristic equation 
and the last two follow from z(s,t) = u(x(s,t), y(s,t)). This system implies 
(P = Ug) Zs + (Q = Uy) Ys = 0 
(P —uz)te + (Q—uy)ye = 0. 


It follows P = ug and Q = uy. 
The initial conditions 


u(a(s,t), y(s, ¢)) = z(s(a,y ; ( ,Y)) = z(s,t) 
Ux(x(s, t), y(s,t = p(s(x,y ) ; )) = p(s,t 
Uy(a(s,t),y(s,t)) = q(s(x,y),t(x,y)) = a(s,t) 
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The uniqueness follows as in the proof of Theorem 2.1. 


Example. A differential equation which occurs in the geometrical optic is 


uz tu, = f(x,y), 


where the positive function f(z, y) is the index of refraction. The level sets 
defined by u(x, y) = const. are called wave fronts. The characteristic curves 
(a(t), y(t)) are the rays of light. If n is a constant, then the rays of light are 
straight lines. In R® the equation is 


u +u, +u2 = f(x,y, z). 


Thus we have to extend the previous theory from R?2 to R", n> 3. 


2.4 Nonlinear equations in R” 
Here we consider the nonlinear differential equation 


where 
L=(%1,...,%n), z=ulz): OCR" HR, p=Vu. 


The following system of 2n +1 ordinary differential equations is called char- 
acteristic system. 


at) = Ver 
2) Spey ok 
p(t) = —-V,F — F.p. 
Let 
xo(s) = (o1(S),.--,2on(s)), 8 = (S1,---,Sn—1), 


be a given regular (n-1)-dimensional C?-hypersurface in R”, i. e., we assume 


Here s € D is a parameter from an (n — 1)-dimensional parameter domain 
D. 

For example, x = xo(s) defines in the three dimensional case a regular 
surface in R?. 
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Assume 
zo(s): Dt R, po(s) = (po1(s),.--,Pon(s)) 


are given sufficiently regular functions. 
The (2n + 1)-vector 


(xo(s), 20(s), po(s)) 


is called initial strip manifold and the condition 


1=1,...,n—1, strip condition. 
The initial strip manifold is said to be noncharacteristic if 


Fy, Fy, Fp, 
Ox01 Ox02 2. 9X0n 
det Osi Osi Osi Zz 0 
9 
0x01 0x02 Oxon 
OSn—1 OSn—1 O8n—1 


where the argument of F,, is the initial strip manifold. 


Initial value problem of Cauchy. Seek a solution z = u(x) of the 
differential equation (2.20) such that the initial manifold is a subset of 
{(x, u(x), Vu(x)): «ce Qh. 


As in the two dimensional case we have under additional regularity as- 
sumptions 


Theorem 2.3. Suppose the initial strip manifold is not characteristic and 
satisfies differential equation (2.20), that is, F(xo(s), zo(s), po(s)) = 0. Then 
there is a neighbourhood of the initial manifold (xo(s), zo(s)) such that there 
exists a unique solution of the Cauchy initial value problem. 


Sketch of proof. Let 
x=x(s,t), z= 2(s,t), p=p(s,t) 


be the solution of the characteristic system and let 
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be the inverse of x = x(s,t) which exists in a neighbourhood of t = 0. Then, 
it turns out that 


z= u(x) := 2($1(@1,...,Un),---; $n—1(41,---,n), t(@1,.-., Ln)) 


is the solution of the problem. 


2.5 Hamilton-Jacobi theory 
The nonlinear equation (2.20) of previous section in one more dimension is 


F(x, se+5Un,Un4+1,%,Pl1,--- Pn; Pn+1) = 0. 


The content of the Hamilton'-Jacobi? theory is the theory of the special 
case 


F = pati t A(a,...,Un,Un41,P1,---;Pn) = 9, (2.21) 


i. e., the equation is linear in py, and does not depend on z explicitly. 
Remark. Formally, one can write equation (2.20) 
F(a1,...,%n,U,Uz,,--+,Ux,) = 0 
as an equation of type (2.21). Set 241 =u and seek u implicitely from 
O(@1,---,;Ln,2n41) = const., 


where ¢ is a function which is defined by a differential equation. 
Assume ¢z,,,, #0, then 


0 = F(a1,...,%n,U, Uay,-+-5 Urn) 
Per Pin 
ome Prns1 


= Gigi igi ieee ar): 


= PUL. see Un, Unt; 


) 


Suppose that G¢,,, # 0, then 


Pants a A(x, ggg Lest 9p Or inves Pinsi)- 


1Hamilton, William Rowan, 1805-1865 
? Jacobi, Carl Gustav, 1805-1851 
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The associated characteristic equations to (2.21) are 


Da) = ap, ey kee ake 
n+l n 

ett) — So Fn = S > piHy, + Pri 
l=1 l=1 


n 


= S> piHy, — 
t=] 


Pn+1(T) = —Fy., — Pepsi 
= — Feng, 
p(T) = —Fy, F.pr 


Set t := @n+41, then we can write partial differential equation (2.21) as 


u + A(2,t, Vou) =0 (2:22) 
and 2n of the characteristic equations are 
a(t) = V>pH (cz, t, p) (2.23) 
p(t) = -VeH(2,t,p). (2.24) 
Here is 
ed (x4,. oe Saas p= (P1,--+5Pn)- 
Let x(t), p(t) be a solution of (2.23) and (2.24), then it follows p/,,,(#) and 


z'(t) from the characteristic equations 


Patilt) = —Hi 
2'(t) p:-VpH — H. 


Definition. The function H(z,t,p) is called Hamilton function, equa- 
tion (2.21) Hamilton-Jacobi equation and the system (2.23), (2.24) canonical 
system to H. 


There is an interesting interplay between the Hamilton-Jacobi equation 
and the canonical system. According to the previous theory we can con- 
struct a solution of the Hamilton-Jacobi equation by using solutions of the 
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canonical system. On the other hand, one obtains from solutions of the 
Hamilton-Jacobi equation also solutions of the canonical system of ordinary 
differential equations. 


Definition. A solution ¢(a;x,t) of the Hamilton-Jacobi equation, where 
a = (a1,...,@,) is an n-tuple of real parameters, is called a complete integral 
of the Hamilton-Jacobi equation if 


det (¢x,a, rai a 0. 


Remark. If u is a solution of the Hamilton-Jacobi equation, then also 
u + const. 


Theorem 2.4 (Jacobi). Assume 
u = o(a;x,t) +c, c=const., PE C? in its arguments, 
is a complete integral. Then one obtains by solving of 
bi = $a; (a; 2, t) 


with respect to x; = x1(a,b,t), where bj i=1,...,n are given real constants, 
and then by setting 
i= be (a; Tay bit); 2) 


a 2n-parameter family of solutions of the canonical system. 


Proof. Let 
xi(a,b;t), L=1,...,n, 


be the solution of the above system. The solution exists since ¢ is a complete 
integral by assumption. Set 


Di Get) = Oe (GAO, Et) em, Py 


We will show that x and p solves the canonical system. Differentiating ¢,, = 
b; with respect to ¢ and the Hamilton-Jacobi equation ¢; + H(x,t, Vid) = 0 
with respect to a;, we obtain fori =1,...,n 


” Ox 
Pta; el Ss Papa oa = 0 
k=1 


I 
2° 


n 
Pta; oF S- bxzpa;H pr, 


k=1 
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Since ¢ is a complete integral it follows for k = 1,...,n 
Ox, 
Ee — Hon 


Along a trajectory, i. e., where a, b are fixed, it is 2%" = x/,(t). Thus 
y Ot k 


rh (t)= Hy, . 


Now we differentiate p;(a,b;t) with respect to t and ¢ + H(2,t, Vio) = 0 
with respect to x;, and obtain 


Dit) = Paste + > Prey 2h (t) 
k= 


0 = Dajt sir S- a,x, Hy, =F lz, 
k= 
n 


0 = dajt alr S- Cie tO) +E Hy, 


k= 


It follows finally that p/(t) = —Hz,. 
Example: Kepler problem 


The motion of a mass point in a central field takes place in a plane, 
say the («,y)-plane, see Figure 2.13, and satisfies the system of ordinary 
differential equations of second order 


aE) = et = Vig 

where 

k2 

[72 4 yr 

Here we assume that k? is a positive constant and that the mass point is 
attracted of the origin. In the case that it is pushed one has to replace U 
by —U. See Landau and Lifschitz [12], Vol 1, for example, for the related 
physics. 

Set 


U(x, y) = 


and 
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7 (x().y(®) 


(U,.U,) 


0) 


Figure 2.13: Motion in a central field 


then 
a(t) = Hp, y(t) = Hq 
p(t) = —Hz, q(t) =—Hy. 
The associated Hamilton-Jacobi equation is 
Ke 


Ve 


gi + 5(br + 346) = —- (2.25) 


Now we will seek a complete integral of (2.25) by making the ansatz 
oy = —a=const. dg = —2 = const. (2.26) 


and obtain from (2.25) that 


r 2 2 
o=+/ 2a 4 ak e: dp + c(t, 0). 
TO p p 


From ansatz (2.26) it follows 


c(t, 0) = —at — 36. 
Therefore we have a two parameter family of solutions 


b = (a, 8; 4,7, t) 
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of the Hamilton-Jacobi equation. This solution is a complete integral, see 
an exercise. According to the theorem of Jacobi set 


ga = —to, $8 = —%. 
Then 


eae [ dp 
o= : 


The inverse function r = r(t), r(0) = ro, is the r-coordinate depending on 


time t, and 
Tr dp 
9 =8 | ; as 
ie eae 


Substitution 7 = p~! yields 


1/r 
9-% = -6 [ al 
1/ro / 2a + 2k27 — B27 


2 
tl oid n( 2a! 
= — arcsin SF peTeT Ts + arcsin =e ae . 
1+ 2 14 er 


Set 
bil 
6; = 0 +arcsin | “22 
1+ 20? 
and 
B? 3 2a 32 
p= }2” CAS 1 + kA ’ 
then 
Po 
6— 4, =~ aresin (# ; ) , 
€ 
It follows 
r=r(60) P 


1 — e? sin(@ — 61)’ 
which is the polar equation of conic sections. It defines an ellipse if0 < « < 1, 
a parabola if ¢ = 1 and a hyperbola if ¢ > 1, see Figure 2.14 for the case 
of an ellipse, where the origin of the coordinate system is one of the focal 
points of the ellipse. 


For another application of the Jacobi theorem see Courant and Hilbert [4], 
Vol. 2, pp. 94, where geodedics on an ellipsoid are studied. 
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1+ 2 


Figure 2.14: The case of an ellipse 


2.6 Exercises 


1. 


Suppose u: R? + R is a solution of 
a(x, y)Ug + O(a, y)Uy = 0. 


Show that for arbitrary H € C1 also H(u) is a solution. 


. Find a solution u # const. of 


Ug + Uy = 0 
such that 
graph(u) := {(x,y,z) € R®: z=u(x,y), (x,y) € R*} 


contains the straight line (0,0,1) + s(1,1,0), sER. 


. Let (x,y) be a solution of 


a1 (@, YU + a2(a, y)Uy = 0. 


Prove that level curves Sc := {(@,y): (x,y) = C = const.} are 
characteristic curves, provided that V¢ 4 0 and (a1, a2) 4 (0,0). 


60 


10. 


11. 
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. Prove Proposition 2.2. 


. Find two different solutions of the initial value problem 


Ug + Uy = 1, 


where the initial data are xo(s) = s, yo(s) = 8, zo(s) =. 


Hint: (xo, yo) is a characteristic curve. 


. Solve the initial value problem 


LU + YUy = U 


with initial data xo(s) = s, yo(s) = 1, zo(s), where zo is given. 


. Solve the initial value problem 


—EUy + YUy = ru’, 


x9(s) = 8, yo(s) =1, 20(s) =e. 


. Solve the initial value problem 


UU, + Uy = 1, 


xro(s) = 8, yo(s) = s, z2o(s) = 8/2 if0<s<1. 


. Solve the initial value problem 


UUg + UUy = 2, 
xo(s) = 8, yo(s) = 1, zo(s) =14+sif0<s<1. 


Solve the initial value problem u2 + u, = 1+ 2 with given initial data 


xo(s) =0, yo(s) =s, uo(s) =1, po(s) =1, qo(s) =0, -co<s< om. 
Find the solution (2, y) of 

(2 —Y)tiz + 2yuy = 32 
such that the surface defined by z = ®(x, y) contains the curve 


C: xo(s)= 8, yo(s)=1, zo(s) =0, s ER. 


2.6. 


12. 


13. 


14. 


15. 


16. 


17. 


18. 
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Solve the following initial problem of chemical kinetics. 
Ug + Uy = (koe + ky) (1- u)’, z>0, y>0 


with the initial data u(x, 0) = 0, u(0, y) = uo(y), where up, 0 < up < 1, 
is given. 
Solve the Riemann problem 
Ur, +Uro = O 
u(t1,0) = (21) 


in Q1 = {(21, 22) € R?: 21 > x2} and in Q2 = {(a1, 272) € R?: 21 < 
x2}, where 


with constants uw ~ Ur. 


Determine the opening angle of the Monge cone, i. e., the angle be- 
tween the axis and the apothem (in German: Mantellinie) of the cone, 
for equation 

iG + ic = f(x,y, u), 


where f > 0. 


Solve the initial value problem 


where xo(0) = acos6, yo(@) = asin@, z9(0) = 1, po(O) = cos@, 
go(@) = sin@ if 0 < 6 < 2m, a= const. > 0. 


Show that the integral ¢(a,(;6,1r,t), see the Kepler problem, is a 
complete integral. 


a) Show that S= /az+Vl—ay+8,a, BER, 0<a<l,isa 
complete integral of S, — ,/1—.$? =0. 
b) Find the envelope of this family of solutions. 


Determine the length of the half axis of the ellipse 


P 


= e* sin(@ — 6p) 


»,O<e<l. 
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19. Find the Hamilton function H(«,p) of the Hamilton-Jacobi-Bellman 
differential equation if h = 0 and f = Ax + Ba, where A, B are 
constant and real matrices, A: R™ +> R"”, B is an orthogonal real 
n x n-Matrix and p € R” is given. The set of admissible controls is 
given by 


U={aeER”": Soa} <1}. 
i=1 


Remark. The Hamilton-Jacobi-Bellman equation is formally the Hamilton- 
Jacobi equation uz; + H(z, Vu) = 0, where the Hamilton function is 
defined by 


H(x,p) := min (f(x, a) -p+h(a,a)), 


f(x,a@) and h(x,q) are given. See for example, Evans [5], Chapter 10. 


Chapter 3 


Classification 


Different types of problems in physics, for example, correspond different 
types of partial differential equations. The methods how to solve these 
equations differ from type to type. 

The classification of differential equations follows from one single ques- 
tion: Can we calculate formally the solution if sufficiently many initial data 
are given? Consider the initial problem for an ordinary differential equa- 
tion y/(x) = f(x, y(x)), y(xo) = yo. Then one can determine formally the 
solution, provided the function f(x,y) is sufficiently regular. The solution 
of the initial value problem is formally given by a power series. This formal 
solution is a solution of the problem if f(z, y) is real analytic according to 
a theorem of Cauchy. In the case of partial differential equations the re- 
lated theorem is the Theorem of Cauchy-Kowalevskaya. Even in the case 
of ordinary differential equations the situation is more complicated if y’ is 
implicitly defined, i. e., the differential equation is F(x, y(x), y’(x)) = 0 for 
a given function F’. 


3.1 Linear equations of second order 

The general nonlinear partial differential equation of second order is 
F(a,u, Du, D?u) = 0, 

where z € R",u: QCR" HR, Du= Vu and D?u stands for all second 


derivatives. The function F is given and sufficiently regular with respect to 
its 2n +1-+n? arguments. 
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In this section we consider the case 


S- a’*(x)ue,e, + f(a,u, Vu) = 0. (3.1) 
ik=1 


’ 


The equation is linear if 
f= 0'(w)us, + e(x)u+ d(z). 
i=l 


Concerning the classification the main part 


n 


ye a’ (x) Uric, 


i,k=1 
plays the essential role. Suppose u € C?, then we can assume, without 
restriction of generality, that a’* = a’, since 


n 


n 
S- a a, = oS (a) tages 


i,k=1 i,k=1 


where i 
(a’*)* _ 5 (a"" +4 ae", 


Consider a hypersurface S in R” defined implicitly by x(x) = 0, Vy 4 0, 
see Figure 3.1 
Assume u and Vu are given on S. 


Problem: Can we calculate all other derivatives of u on S by using differ- 
ential equation (3.1) and the given data? 


We will find an answer if we map S onto a hyperplane Sp by a mapping 


Nie So XB Ba) 
MN = A(T1,.-.,en), t=1,...,n-1, 


for functions A; such that 


O(A1, ice! An) 


act O(a1,..-,2n) 


#0 


in Q Cc R”. It is assumed that y and A, are sufficiently regular. Such a 
mapping \ = A(z) exists, see an exercise. 
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X3 


Va RY 


ay 


Figure 3.1: Initial manifold S 
The above transform maps S onto a subset of the hyperplane defined by 
An = 0, see Figure 3.2. 


We will write the differential equation in these new coordinates. Here we 
use Einstein’s convention, i. e., we add terms with repeating indices. Since 


where x = (%,...,%n) and A = (Aq,..., An), we get 


Or 
Ux; = WV); Ox,’ (3.2) 
- =. > Ori OAL . OY; 

Oh EN Ox; OLE re Ox; OLK 

Thus, differential equation (3.1) in the new coordinates is given by 
OX; OA 
aa) Bag + terms known on So = 0. 

Since vy, (A1,---,An—1,0), k =1,...,n, are known, see (3.2), it follows that 
Ud, dp 6 = 1,...,n2—-1, are known on Sp. Thus we know all second derivatives 


¥,,; On So with the only exception of v),,),,- 
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Ay 


Figure 3.2: Transformed flat manifold So 


We recall that, provided v is sufficiently regular, 
Wag npAlip io 4 Appi) 


is the limit of 


Uy, (Ar, 2 ADEA, Al41) oan An—1; 0) = Un, (AL, say ATs AL+1) aesecs An—1, 9) 
h 


ash — 0. 
Thus the differential equation can be written as 


3 a)* (ar) Orn OAn Ud» An = terms known on So. 
Tk OL; Orr 


It follows that we can calculate v,,),, if 


n 


Sa (2)Xa:Xa; #0 (3.3) 


t,j=1 


on S. This is a condition for the given equation and for the given surface S. 
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Definition. The differential equation 


n 


S- BP (2) NaN, =0 


t,j=1 


is called characteristic differential equation associated to the given differen- 
tial equation (3.1). 

If vy, Vx # 0, is a solution of the characteristic differential equation, then 
the surface defined by x = 0 is called characteristic surface. 


Remark. The condition (3.3) is satisfied for each y with Vx 4 0 if the 
quadratic matrix (a‘/(x)) is positive or negative definite for each x € Q, 
which is equivalent to the property that all eigenvalues are different from 
zero and have the same sign. This follows since there is a A(x) > 0 such 
that, in the case that the matrix (a’’) is poitive definite, 


n 


de a4 (#) GG; = AC@YICP 


4,j=1 


for all ¢ € R”. Here and in the following we assume that the matrix (a’’) is 
real and symmetric. 


The characterization of differential equation (3.1) follows from the signs of 
the eigenvalues of (a’)(a)). 


Definition. Differential equation (3.1) is said to be of type (a,G,y) at 
x € © if a eigenvalues of (a’)(x) are positive, 3 eigenvalues are negative 
and ¥y eigenvalues are zero (a+ 8+y =n). 

In particular, equation is called 

elliptic if it is of type (n,0,0) or of type (0,7, 0), i. e., all eigenvalues are 
different from zero and have the same sign, 

parabolic if it is of type (n—1, 0, 1) or of type (0,n—1, 1), i. e., one eigenvalue 
is zero and all the others are different from zero and have the same sign, 
hyperbolic if it is of type (n — 1,1,0) or of type (1,n — 1,0), i. e., all 
eigenvalues are different from zero and one eigenvalue has another sign than 
all the others. 


68 CHAPTER 3. CLASSIFICATION 


Remarks: 


1. According to this definition there are other types aside from elliptic, 
parabolic or hyperbolic equations. 


2. The classification depends in general on x € Q. An example is the 
Tricomi equation, which appears in the theory of transsonic flows, 


YUge + Uyy = 0. 


This equation is elliptic if y > 0, parabolic if y = 0 and hyperbolic for y < 0. 


Examples: 
1. The Laplace equation in R? is Au = 0, where 
Aw := Ugg + Uyy + Uzz- 


This equation is elliptic. Thus for each manifold S given by {(z,y,z) : 
x(x, y,z) = 0}, where x is an arbitrary sufficiently regular function such 
that Vy £4 0, all derivatives of wu are known on S, provided u and Vu are 
known on S. 


2. The wave equation ut = Urx + Uyy + Uzz, Where u = u(z,y,z,t), is 
hyperbolic. Such a type describes oscillations of mechanical structures, for 
example. 


3. The heat equation uz = Ugg tUyy+uzz, where u = u(x, y, Z,t), is parabolic. 
It describes, for example, the propagation of heat in a domain. 


4. Consider the case that the (real) coefficients a’? in equation (3.1) are 
constant. We recall that the matrix A = (a’) is symmetric, i. e., A? = A. 
In this case, the transform to principle axis leads to a normal form from 
which the classification of the equation is obviously. Let U be the associated 
orthogonal matrix, then 
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Here is U = (21,...,2n), where z,1 =1,...,n, is an orthonormal system of 
eigenvectors to the eigenvalues ;. 
Set y =U? x and v(y) = u(Uy), then 


n n 
s a tee = Se AVyiy;- (3.4) 
i=l 


t,j=1 


3.1.1 Normal form in two variables 
Consider the differential equation 
A(X, Y)Uge + 20(@, y)Ury + C(x, y)Uyy + terms of lower order =0 (3.5) 
in Q C R?. The associated characteristic differential equation is 
ax? + 2bX2Xy + x, =0. (3.6) 


We show that an appropriate coordinate transform will simplify equation (3.5) 
sometimes in such a way that we can solve the transformed equation explic- 
itly. 

Let z = (x,y) be a solution of (3.6). Consider the level sets {(a,y) : 
(x,y) = const.} and assume dy # 0 at a point (xo, yo) of the level set. 
Then there is a function y(a) defined in a neighbourhood of xo such that 
o(«, y(x)) = const. It follows 


which implies, see the characteristic equation (3.6), 
ay’? — 2by’ +ce=0. (3.7) 


Then, provided a 4 0, we can calculate jz := y’ from the (known) coefficients 


a, b and c: ; 
M12 = = (0 + Vb? — ac) ‘ (3.8) 


These solutions are real if and only of ac — b2 <0. 
Equation (3.5) is hyperbolic if ac — b? < 0, parabolic if ac — b? = 0 and 
elliptic if ac— b? > 0. This follows from an easy discussion of the eigenvalues 


of the matrix 
a b 
bc}? 


see an exercise. 
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Normal form of a hyperbolic equation 


Let ¢ and w are solutions of the characteristic equation (3.6) such that 


¥=M = a“ 
dy 
Y= p2 = = 
Py’ 


where ji; and j1g are given by (3.8). Thus ¢ and w are solutions of the linear 
homogeneous equations of first order 


Pe + br (@, y) by = 0 (3.9) 
letounaes a (3.10) 


Assume ¢$(2, y), w(x, y) are solutions such that V¢ 4 0 and Vw 4 0, see an 
exercise for the existence of such solutions. 

Consider two families of level sets defined by ¢(z, y) = a and w(a, y) = , 
see Figure 3.3. 


O(% y)= a 1 


P(x, y)= 10 2 


YWXY= B 1 


WXy)= Bo 


Figure 3.3: Level sets 


These level sets are characteristic curves of the partial differential equa- 
tions (3.9) and (3.10), respectively, see an exercise of the previous chapter. 


Lemma. (i) Curves from different families can not touch each other. 


(ii) bay _ by Pex # 0. 
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Proof. (i): 


/ 


2 
¥o— Y= be-M = awe ac # 0. 


(ii): 


Proposition 3.1. The mapping € = ¢(2,y), 7 = V(a,y) transforms equa- 
tion (3.5) into 
Ven = lower order terms, (3.11) 


where v(€,n) = u(x(E, 7), y(E,7))- 


Proof. The proof follows from a straightforward calculation. 


Up = Vehy + Unde 
Uy = Veby + UnYy 
Une = vech> + vendre Une lower order terms 
Ury = veehrby + Ven baby + bye) + Unn ey + lower order terms 
Uyy = vege + 2venbyPy Um lower order terms. 
Thus 
AUge + 2uey + CUyy = AVEe + 2ZVE_ + -YUnn + L.0.t., 
where 
a: = adgz + wdbrdy + c¢3 
B: = abzipy + (bay + bye) + ehytby 
ys = ap? + bray + cy}. 


The coefficients @ and y are zero since ¢ and w are solutions of the charac- 
teristic equation. Since 


ay — ? = (ac _ b°) (bay = dbyWn)”, 


it follows from the above lemma that the coefficient G is different from zero. 


Example: Consider the differential equation 


Ugg — Uyy = 0. 
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The associated characteristic differential equation is 


Xz — Xy =0. 


Since 44 = —1 and p2 = 1, the functions ¢ and w satisfy differential equa- 
tions 


da + by = 0 
V2 — Py = 0. 


Solutions with Vé 4 0 and Vw #0 are 


Thus the mapping 
f=@£-y, n=rt+y 


leads to the simple equation 


ven (€, n) = 0. 


Assume v € C? is a solution, then ve = fi(€) for an arbitrary C’ function 
fi(&). It follows 


€ 
v(E,n) = | fila) da + 9(n), 


where g is an arbitrary C? function. Thus each C?-solution of the differen- 
tial equation can be written as 


(x) v(m) = f(g) + 9(n), 
where f, g € C?. On the other hand, for arbitrary C?-functions f, g the 


function (x) is a solution of the differential equation vg, = 0. Consequently 
each C?-solution of the original equation ur, — Uyy = 0 is given by 


u(xz,y) = f(w@—y) + g9(@+y), 


where f, g € C?. 
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3.2  Quasilinear equations of second order 


Here we consider the equation 
n 
SS a (a, u, Vu)ue;«; + (x, u, Vu) = 0 (3.12) 
i,j=l 
in a domain 2 C R”, where u: Q+> R. We assume that a = a", 
As in the previous section we can derive the characteristic equation 


n 


S> al (x, u, Vu) Xai X2; =0. 
ij=l 


In contrast to linear equations, solutions of the characteristic equation de- 
pend on the solution considered. 


3.2.1 Quasilinear elliptic equations 


There is a large class of quasilinear equations such that the associated char- 
acteristic equation has no solution y, Vy # 0. 
Set 
U = {(2,z,p): c€0, zER, pe R"}. 


Definition. The quasilinear equation (3.12) is called elliptic if the matrix 
(a (x, z,p)) is positive definite for each (a, z,p) € U. 


Assume equation (3.12) is elliptic and let A(x, z,p) be the minimum and 
A(a,z,p) the maximum of the eigenvalues of (a’/), then 


0< Nz, z,p)|¢|? S ye a) (a, 2, D)GCj < A(x, z,p)|¢|? 


igj=l 
for all ¢ € R”. 


Definition. Equation (3.12) is called uniformly elliptic if A/X is uniformly 
bounded in U. 


An important class of elliptic equations which are not uniformly elliptic 
(nonuniformly elliptic) is 


“0 Ue, 
“2 + lower order terms = 0. (3.13) 
» Ox; (= + =| 
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The main part is the minimal surface operator (left hand side of the minimal 
surface equation). The coefficients a are 


ae _1/2 Dip; 
a") (x, z,p) = ( + Ip!) / (5; — tt) 5 


6;; denotes the Kronecker delta symbol. It follows that 


1 1 


(1+ |p|2)°/?’ (1 + |p|2)/? 


Thus equation (3.13) is not uniformly elliptic. 

The behaviour of solutions of uniformly elliptic equations is similar to 
linear elliptic equations in contrast to the behaviour of solutions of nonuni- 
formly elliptic equations. Typical examples for nonuniformly elliptic equa- 
tions are the minimal surface equation and the capillary equation. 


3.3. Systems of first order 


Consider the quasilinear system 


S 7 A*(a, u)ttu, + d(x, u) = 0, (3.14) 
k=1 


where A’ are m x m-matrices, sufficiently regular with respect to their ar- 
guments, and 


UL U1,ap by 
Uu= ; » Ux, = ; , b= 


Um Um,rp bm 


We ask the same question as above: can we calculate all derivatives of u 
in a neighbourhood of a given hypersurface S in R” defined by x(x) = 0, 
Vx #0, provided u(x) is given on S? 

For an answer we map S onto a flat surface Sp by using the mapping 
A = A(x) of Section 3.1 and write equation (3.14) in new coordinates. Set 
v(A) = u(a(A)), then 


n 
S- A¥(x,u)Xx,Ud, = terms known on So. 
k=1 
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We can solve this system with respect to v),,, provided that 
n 
det (>: A¥ (a, orn) #0 
k=1 
on S. 
Definition. Equation 
n 
det (>: A* (x, ox) = 0 
k=1 


is called characteristic equation associated to equation (3.14) and a surface 
S: y(x) = 0, defined by a solution x, Vx 4 0, of this characteristic equation 
is said to be characteristic surface. 


Set 
n 
C(ax,u,¢) = det (>: A* (x, v6 

for ¢ € R”. 

Definition. (i) The system (3.14) is hyperbolic at (x,u(x)) if there is a 
regular linear mapping ¢ = Qn, where 7 = (m,---,%m-—1,4), Such that there 
exists m real roots Kp = Kp (@,Uu(@),M,---,Mm—1), k= 1,...,m, of 

D(z, u(x), My+++9%Mn-1s k) =0 


for all (m1,.-.,%,-1), where 


D(a, u(x), Mss +5Mn-1; kK) = C(x, u(x), 2, Qn). 


(ii) System (3.14) is parabolic if there exists a regular linear mapping ¢ = Qn 
such that D is independent of «, i. e., D depends on less than n parameters. 


(iii) System (3.14) is elliptic if C(x, u,¢) = 0 only if ¢ = 0. 


Remark. In the elliptic case all derivatives of the solution can be calculated 
from the given data and the given equation. 
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3.3.1 Examples 


1. Beltrami equations 


Wuz — bug — cy = 0 (3.15) 
Wuy +avz+bvy = 0, (3.16) 


where W, a, b, c are given functions depending of (x,y), W 4 0 and the 


matrix 
a b 
bc 
is positive definite. 


The Beltrami system is a generalization of Cauchy-Riemann equations. 
The function f(z) = u(x, y) + iv(x,y), where z = x + iy, is called a qua- 
siconform mapping, see for example [9], Chapter 12, for an application to 
partial differential equations. 


Set 
1 WwW —b 2 0 = C 
A= (4 a tee W b/}* 


Then the system (3.15), (3.16) can be written as 
Al Ur 4 Ae Uy = 0 
Ux Vy 0 : 


55) Wika: BG ekp 
COM O=| We ali +06 


Thus, 


= W (aC? + 2b61¢2 + ¢¢3), 


which is different from zero if ¢ 4 0 according to the above assumptions. 
Thus the Beltrami system is elliptic. 


2. Maxwell equations 


The Maxwell equations in the isotropic case are 


croty H = AE +c, (3.17) 
c rot, E — pH, (3.18) 
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where 

E = (e1, €2,€3)" electric field strength, e; = e;(a,t), « = (x1, 22, 73), 
H = (hj,ho,h3)’ magnetic field strength, h; = h;(z,t), 

c speed of light, 

A specific conductivity, 

€ dielectricity constant, 

js magnetic permeability. 

Here c, A, € and py are positive constants. 


Set po = Xt, Pi = Xx,;, 1 = 1,...3, then the characteristic differential equa- 
tion is 


epo/c 0 0 0 P3 —P2 
0 epo/e 0 —P3 0 PL 
0 0  epo/c pz Pi De lacy 
0 —p3 pe ppo/c 0 0 
D3 0 —pi 0 ppo/e 0 
—p2 pi 0 0 0  — upo/e 


The following manipulations simplifies this equation: 

(i) multiply the first three columns with ppo/c, 

(ii) multiply the 5th column with —p3 and the the 6th column with po and 
add the sum to the 1st column, 

(iii) multiply the 4th column with p3 and the 6th column with —p, and add 
the sum to the 2th column, 

(iv) multiply the 4th column with —p2 and the 5th column with p; and add 
the sum to the 3th column, 

(v) expand the resulting determinant with respect to the elements of the 
6th, 5th and 4th row. 


We obtain 
q+Pi Pip, — P1p3 
Pip2 q+p3; pops | =0, 
Pips P2P3.,- G+: D3 
where 
au 
Ce —aPO — g 
with g? := Dy + p 2+ p3. The evaluation of the above equation leads to 
(q+ 9°) =0, i. 


ELL 
x? (xd — [Vexl?) = 
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It follows immediately that Maxwell equations are a hyperbolic system, see an 
exercise. There are two solutions of this characteristic equation. The first 
one are characteristic surfaces S(t), defined by x(x,t) = 0, which satisfy 
xz = 0. These surfaces are called stationary waves. The second type of 
characteristic surfaces are defined by solutions of 


EU 

“9 Xt = |Vexl?. 

c 
Functions defined by y = f(n-x—Vt) are solutions of this equation. Here is 
f(s) an arbitrary function with f’(s) 4 0, n is a unit vector and V = c/,/ep. 
The associated characteristic surfaces S(t) are defined by 

x(a,t) = f(n-x—Vt)=0, 

here we assume that 0 is in he range of f: Rt R. Thus, S(t) is defined 
by n-2—Vt=c, where c is a fixed constant. It follows that the planes S(t) 


with normal n move with speed V in direction of n, see Figure 3.4 


X9 


S(t) 
d(t) = 


SO) 
Figure 3.4: d'(t) is the speed of plane waves 
V is called speed of the plane wave S(t). 


Remark. According to the previous discussions, singularities of a solution 
of Maxwell equations are located at most on characteristic surfaces. 


A special case of Maxwell equations are the telegraph equations, which 
follow from Maxwell equations if div FE = 0 and div H = 0,i. e«., FE and 
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HT are fields free of sources. In fact, it is sufficient to assume that this 
assumption is satisfied at a fixed time tg only, see an exercise. 
Since 

rot, rot, A = grad, div; A—A,A 


for each C?-vector field A, it follows from Maxwell equations the uncoupled 
system 


€ » 
A,E = En + SF 


€ Xr 
ee: ee Hy + SH 


3. Equations of gas dynamics 
Consider the following quasilinear equations of first order. 
1 
ut (v- Vz) v+-Vep =f (Euler equations). 
p 


Here is 
v = (v1, v2, v3) the vector of speed, vj = v;(x, t), v = (1, 22, 23), 
p pressure, p = (a, ¢), 


p density, o = p(e,2), 
f = (fi, fo, fg) density of the external force, f; = f;(x, t), 
(vu: Vz)v = (uv: Vev1, U+ Vev2,U + Vezv3))?. 


The second equation is 
prtvu:-Vept+p dive v=0 (conservation of mass). 
Assume the gas is compressible and that there is a function (state equation) 
P=P(p), 
where p'(p) > 0 if p > 0. Then the above system of four equations is 
n+ (v- V+ o(e)Ve = f (3.19) 
pte divu+u-Vp = 0, (3.20) 


where V = Vz and div = divyg, i. e., these operators apply on the spatial 
variables only. 
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The characteristic differential equation is here 


Z 0 0 FPXa 
d 

0 x 0 SP Xe 

O. 0) A ops 


’ 
dt 


PXe, PXep PXog 


where 


Evaluating the determinant, we get the characteristic differential equation 


(2) (2) -Herrar)=0 ay 


This equation implies consequences for the speed of the characteristic sur- 
faces as the following consideration shows. 

Consider a family S(t) of surfaces in R? defined by y(x,t) = c, where 
x € R? and cis a fixed constant. As usually, we assume that Vzy 4 0. One 
of the two normals on S(t) at a point of the surface S(t) is given by, see an 
exercise, 


n= Vex 
[Wane 
Let Qo € S(to) and let Q1 € S(t1) be a point on the line defined by Qo + sn, 


where n is the normal (3.22 on S(to) at Qo and to < t1, t1 — tp small, see 
Figure 3.5. 


(3.22) 


S(t) 


S(ty) 


Figure 3.5: Definition of the speed of a surface 
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Definition. The limit 
ait 1Q1 — Qol 
im 


tito §=t, — to 


P 


is called speed of the surface S(t). 
Proposition 3.2. The speed of the surface S(t) is 


Pa, (3.23) 


Proof. The proof follows from y(Qo, to) = 0 and x(Qo + dn, tp + At) = 0, 
where d = |Q1 — Qo| and At = ty — to. 


Set vp, := u-n which is the component of the velocity vector in direction 
n. From (3.22) we get 


Definition. V := P — vpn, the difference of the speed of the surface and the 
speed of liquid particles, is called relative speed. 


Figure 3.6: Definition of relative speed 


Using the above formulas for P and vy, it follows 


Xt CFV EX = 1 dx 
|Vxx| |Vrx| |Vxx| dt” 


V=P-v,= 


Then, we obtain from the characteristic equation (3.21) that 


V7\Vexl? (V7|Vexl? — p’(e)|Vexl?) = 0. 
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An interesting conclusion is that there are two relative speeds: V = 0 or 
V? = p'(p). 


Definition. \/p’/(p) is called speed of sound . 


3.4 Systems of second order 


Here we consider the system 


n 


>a A™ (x, u, Vu)Uz,2, + lower order terms = 0, (3.24) 
kjl=1 
where A*™ are (m x m) matrices and u = (u1,...,Um)’. We assume A* = 


A'®, which is no restriction of generality provided u € C? is satisfied. As in 
the previous sections, the classification follows from the question whether or 
not we can calculate formally the solution from the differential equations, 
if sufficiently many data are given on an initial manifold. Let the initial 
manifold S be given by x(x) = 0 and assume that Vy 4 0. The mapping 
x = (A), see previous sections, leads to 


n 
y AP sone: U\»An = terms known on S, 
kl=1 


where v(A) = u(x())). 


The characteristic equation is here 
n 
det S- AP yes =0. 
kl=1 


If there is a solution x with Vy 4 0, then it is possible that second derivatives 
are not continuous in a neighbourhood of S. 


Definition. The system is called elliptic if 
det | $7 AM G.¢ | #0 
kl=1 


for all ¢ € R", ¢ £0. 
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3.4.1 Examples 


1. Navier-Stokes equations 


The Navier-Stokes system for a viscous incompressible liquid is 


1 
ut(v-Va)u = -—-Vep+yAcv 
p 
div, v = O, 


where p is the (constant and positive) density of liquid, 
+ is the (constant and positive) viscosity of liquid, 
v = v(a,t) velocity vector of liquid particles, 2 € R° or in R?, 
p = p(x, t) pressure. 
The problem is to find solutions v, p of the above system. 


2. Linear elasticity 


Consider the system 
Oru 
POP 
Here is, in the case of an elastic body in R®, 
u(a,t) = (u(x,t), u2(a, t), u3(x, t)) displacement vector, 
f(«,t) density of external force, 
p (constant) density, 
A, - (positive) Lamé constants. 


= pA,ut (A+ p)V2(dive u) + f. (3.25) 


The characteristic equation is det C = 0, where the entries of the matrix 
C’ are given by 


Cig = (A+ W)XaiXaj + ig (u|Vexl? — px?) - 
The characteristic equation is 


((A + 21) Vexl? — px?) (ulVexl? — px?)” = 0. 


It follows that two different speeds P of characteristic surfaces S(t), defined 
by x(a, t) = const., are possible, namely 


A+ 2 
Py = rela and Py = i 
p p 


We recall that P = —x/|Vezx\. 
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3.5 Theorem of Cauchy-Kovalevskaya 


Consider the quasilinear system of first order (3.14) of Section 3.3. Assume 
an initial manifolds S is given by v(x) = 0, Vx 4 0, and suppose that x is 
not characteristic. Then, see Section 3.3, the system (3.14) can be written 
as 


n-1 
a a’(x,U)Uy, + b(a, u) (3.26) 
i=l 
U(a1,.+-)2a-1,0) = f(e1,.-+)2n-1) (3.27) 
Here is u = (ui,-.-,Um)2, b = (b1,---, bn)? and a’ are (mxm)-matrices. We 


assume a’, b and f are in C® with respect to their arguments. From (3.26) 
and (3.27) it follows that we can calculate formally all derivatives D°u in a 
neigbourhood of the plane {x : x, = 0}, in particular in a neighbourhood 
of 0 € R”. Thus we have a formal power series of u(x) at x = 0: 


1 
u(a) ~ ai uO)a". 


For notations and definitions used here and in the following see the appendix 
to this section. 


Then, as usually, two questions arise: 
(i) Does the power series converge in a neighbourhood of 0 € R"? 


(ii) Is a convergent power series a solution of the initial value problem (3.26), 
(3.27)? 


Remark. Quite different to this power series method is the method of 
asymptotic expansions. Here one is interested in a good approximation of 
an unknown solution of an equation by a finite sum See 9 i(x) of functions 
¢;. In general, the infinite sum )77°) ¢;(x) does not converge, in contrast 
to the power series method of this section. See [15] for some asymptotic 
formulas in capillarity. 


Theorem 3.1 (Cauchy-Kovalevskaya). There is a neighbourhood of 0 € R” 
such there is a real analytic solution of the initial value problem (3.26), 
(8.27). This solution is unique in the class of real analytic functions. 
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Proof. The proof is taken from F. John [10]. We introduce u— f as the new 
solution for which we are looking at and we add a new coordinate u* to the 
solution vector by setting u*(x) = x». Then 


Uz, = 1, uy, =0, K=1,...,.n—1, w*(a1,...,Fn-1,0) =0 


and the extended system (3.26), (3.27) is 


U1,an Ula; by 
. n=l ma 0 
u - 2, ( 0. 0 ) Uw x b 
M,Ln i=l MX; m 
em Un, 1 
where the associated initial condition is u(x1,...,2%p—-1,0) = 0. The new u 
is u = (u1,...,Um)*, the new a! are a*(11,...,%n—1,U1,..., Um, u*) and the 
new 6 is b= (21,...,2n—1,U1,---,Um,U*)?. 


Thus we are led to an initial value problem of the type 


n—-1 N 
tiie, =) YY cig ue, + O,(z), P= 1. N (3.28) 
i=1 k=1 
Ue) = Vat t= 0, (3.29) 
where j =1,...,N and z = (a,...,%n-1,t1,---, UN). 


The point here is that a, and b; are independent of z,. This fact 
simplifies the proof of the theorem. 

From (3.28) and (3.29) we can calculate formally all D’u;. Then we 
have formal power series for w;: 


uj(a) ~ So Dae, 


Qa 


where i 
cl) = iP us (0). 


We will show that these power series are (absolutely) convergent in a neigh- 
bourhood of 0 € R”, i. e., they are real analytic functions, see the appendix 
for the definition of real analytic functions. Inserting these functions into 
the left and into the right hand side of (3.28) we obtain on the right and on 
the left hand side real analytic functions. This follows since compositions 
of real analytic functions are real analytic again, see Proposition A7 of the 
appendix to this section. The resulting power series on the left and on the 
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right have the same coefficients caused by the calculation of the derivatives 
D°u;(0) from (3.28). It follows that u;(x), 7 = 1,...,n, defined by its 
formal power series are solutions of the initial value problem (3.28), (3.29). 


Pa ee a 
NOR Otnn aA 


Lemma A. Assume u€ C™ in a neighbourhood of 0 € R”. Then 


Set 


Druj0) = Pa (ai Ob(0), 


where |3|, \y| < |a| and Py are polynomials in the indicated arguments with 
nonnegative integers as coefficients which are independent of a’ and of 


b. 


Proof. It follows from equation (3.28) that 
D,Duj(0) = Py(d7a',,(0), db; (0), D°ux(0)). (3.30) 
Here is D, = 0/Ox,, and a, (, 7, 6 satisfy the inequalities 
IS], ly Sled, [6] < lal +1, 


and, which is essential in the proof, the last coordinates in the multi-indices 
a = (a1,...,Qn), 6 = (01,...,6n) satisfy bn <Q» since the right hand side 
of (3.28) is independent of x. Moreover, it follows from (3.28) that the 
polynomials P, have integers as coefficients. The initial condition (3.29) 
implies 

Du; (0). = 0, (3.31) 
where a = (a1,...,Qn—1,0), that is, a, = 0. Then, the proof is by induction 
with respect to a,. The induction starts with a, = 0, then we replace 
D°ux(0) in the right hand side of (3.30) by (3.31), that is by zero. Then it 
follows from (3.30) that 


D®uj(0) = Pa(d?ai, (0), d%;(0), D?up(0)), 


where @ = (Q1,...,Q@n—1, 1). 


Definition. Let f = (fi,..., fm), F = (fi,..., Fm), fi = fix), Fi = Fi(z), 
and f, Fe C%. We say f is majorized by F if 


|\D°#,(0)| <D°R.(0), FH 14m 
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for all a. We write f << F, if f is majorized by F’. 


Definition. The initial value problem 


n-1 N 

Us, by S Ai,(2)U ke, + B;(2) (3.32) 
i=1 k=1 

U;(c) = 0 ifa,=0, (3.33) 


j=l,...,N, A‘, B; real analytic, ist called majorizing problem to (3.28), 
(3.29) if | | 
Qin << Aj, and bj << Bj. 


Lemma B. The formal power series 
1 D@ a 
Ss al u;(O)x ) 
a 


where Du; (0) are defined in Lemma A, is convergent in a neighbourhood of 
0 € R” if there exists a majorizing problem which has a real analytic solution 
U in x =0, and 

|D®uj(0)| < D°U;(0). 


Proof. It follows from Lemma A and from the assumption of Lemma B that 


Pr, (|a?a4(0)],1a%5(0)1) 


Pa (|d°Aig(0)], |” B;(0)|) = D°U;(0). 


|D%u;(0)| 


IA 


IA 


The formal power series 
1 
2s aie ug (O)a®, 
a 
is convergent since 
1 1 
Y /P uj (0)2"] < ) {D2U;(0)|2"| 
a a 


The right hand side is convergent in a neighbourhood of x € R” by assump- 
tion. 


Lemma C. There is a majorising problem which has a real analytic solution. 
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Proof. Since ai,(z), b;(z) are real analytic in a neighbourhood of z = 0 it 
follows from Proposition A5 of the appendix to this section that there are 


positive constants M and r such that all these functions are majorized by 


Mr 


T— 21 — «1.7 ZN4n-1 


Thus a majorizing problem is 


Mr n-1 N 
Ujuen = r—21—... U, —...—Un (1+ 53%. 


Ln-— 
nd fat REI 


Usa) SS 0S ay = 0) 


GS Vgc voy IN a 
The solution of this problem is 


Ue Gie hee a) OV (Be ay ts eg PIN, 


where V(s,t), s = 2% +...+%n-1, t = Zp, is the solution of the Cauchy 
initial value problem 


Mr 


V(s,0) = 0. 


which has the solution, see an exercise, 


V(9,8) = 5 (r s—V(r—s)? InMNrt) . 


This function is real analytic in (s,t) at (0,0). It follows that U;(a) are also 
real analytic functions. Thus the Cauchy-Kovalevskaya theorem is shown. 
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Examples: 


1. Ordinary differential equations 


Consider the initial value problem 


y(t) = f(x,y(z)) 
y(zo) = Yo, 


where zo € R and yo € R” are given. Assume f(x,y) is real analytic in a 
neighbourhood of (x0, yo) € Rx R”. Then it follows from the above theorem 
that there exists an analytic solution y(x) of the initial value problem in a 
neighbourhood of x9. This solution is unique in the class of analytic func- 
tions according to the theorem of Cauchy-Kovalevskaya. From the Picard- 
Lindelof theorem it follows that this analytic solution is unique even in the 
class of C'-functions. 


2. Partial differential equations of second order 
Consider the boundary value problem for two variables 


Uyy = f (X,Y, U, Ur, Uy, Ure, Ury) 
u(v,0) = (2) 
th, (2,0) = ap(a). 


We assume that ¢, w are analytic in a neighbourhood of « = 0 and that f 
is real analytic in a neighbourhood of 


(0, 0, (0), $'(0), (0), H'(0)). 


There exists a real analytic solution in a neigbourhood of 0 € R? of the above 
initial value problem. 


In particular, there is a real analytic solution in a neigbourhood of 0 € R? 
of the initial value problem 
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The proof follows by writing the above problem as a system. Set p = Uz, 
q = Uy, T = Urge, & = Ugy, t = Uyy, then 


t= f(z, y,U,D, 4,7, 8). 
Set U = (u, p, q,7, 8, t)*, b — (q, 0, t, 0,0, fy a0 fuqt fie and 


00 0 0 0 0O 
00 1 0 0 0 
00 0 0 0 0O 
aad Ie rear ae 
00 0 0 0 1 
0: Oe ge 


Then the rewritten differential equation is the system Uy = AU; + 6 with 
the initial condition 


U(2, 0) = (d(2), ¢'(x), v(z), ¢' (x), y'(x), fo(x)) ) 
where fo(x) = f(a, 0, (2), d(x), v(x), o" (x), w'(z)). 


3.5.1 Appendix: Real analytic functions 
Multi-index notation 


The following multi-index notation simplifies many presentations of formu- 
las. Let x = (%71,...,%n) and 


wu: QCR"R (or R”™ for systems). 


The n-tuple of nonnegative integers (including zero) 


a = (Q1,...,Qn) 
is called multi-index. Set 
lal = ay +...+an 
a! = aylag!+...+ an! 
eo = ep as? s..0 a0,” (or amonom) 
) 
De es = 
A Ox; 
BD. 2 fOeiue De 
Du = (Diu,...,Dyu) = Vu = grad u 
a (ean a2 an — ale 
pD® = D™ps..... D* 


ans a a An’ 
Of) Olson. OLR" 
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Define a partial order by 
a> if and only if a; > 6; for all z. 


Sometimes we use the notations 


where 0, 1 € R”. 


Using this multi-index notion, we have 


1. 
a! 
(@+y%= DP aaeey’, 


where x, y € R” anda, (, ¥ are multi-indices. 


2. Taylor expansion for a polynomial f(x) of degree m: 


here is D® f(0) := (D° f(x)) |x=0. 


3. Let « = (21,...,2%,) and m > 0 an integer, then 
m! 
m _ BeOS 
(ay +... 4+ ay)” = S- ote 
jaj=m 
A. 


5. Leibniz’s rule: 


Df= FiD*NDa): 
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a! 
Dox — grt 8 ifa > B, 
(a — )! 


D°x* = 0 otherwise. 


7. Directional derivative: 


m 
a fle+ty) = Eh wpe tip’, 


lal|=m 
where xz, y € R” andteR. 


8. Taylor’s theorem: Let u € C™+! in a neighbourhood N(y) of y, then, if 
x € N(y), 


where 


1 1 
SS — 4)Maq(m+1) 
Rin al (1—s)™o"™*D (4) at, 


where ®(¢) = u(y +t(x — y)). It follows from 7. that 


1 
Rn =(m+1) >> * (| (1—2)D*u(y + t(@ — y)) ar) Gay 


0 


9. Using multi-index notation, the general linear partial differential equation 
of order m can be written as 


S- Ao(2)D°u = f(#) in Qc R”. 


|a|<m 
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Power series 


Here we collect some definitions and results for power series in R”. 


Definition. Let cg € R (or € R™). The series 


foe) 


Tash | oa 


m=0 jalj=m 
is said to be convergent if 


foe) 


SS cal = s S- Ica 


m=0 \l|a|=m 


is convergent. 


Remark. According to the above definition, a convergent series is abso- 
lutely convergent. It follows that we can rearrange the order of summation. 


Using the above multi-index notation and keeping in mind that we can 
rearrange convergent series, we have 


10. Let x € R”, then 


LS 
Qa 


I| 
1s 
i 
Me 
8 
> 
Nae 


i=1 \a;=0 
~ 1 
(1 —21)(1—22)-...-(1— 2p) 
1 
(1-2) 


provided |a;| < 1 is satisfied for each 7. 
11. Assume x € R” and |x1| + |va|+...+|@n| <1, then 


Jol! a 7 lah wy 
ar? = Ie aie 


Qa j=0 joal=j 


I 
——~ 
8 

a 
+ 
+ 
8 
3 
YS 
w 
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12. Let x € R”, |x| < 1 for all i, and @ is a given multi-index. Then 


a! yo-B — pe 1 
dX a3) a= a 


B! 


13. Let x € R® and |21|+...+|2,| <1. Then 


LS lol! e-8 _ pb 1 


dn (2 — 8)! b= 2p. ay, 


[3|! 


(1-2, —...—2,) +14 * 


Consider the power series 
S- Cot (3.34) 
a 


and assume this series is convergent for a z € R”. Then, by definition, 


p= > lCa||z*| < 00 


a 


and the series (3.34) is uniformly convergent for all « € Q(z), where 
Q(z): |x| <|z| for all 7. 


Thus the power series (3.34) defines a continuous function defined on Q(z), 
according to a theorem of Weierstrass. 

The interior of Q(z) is not empty if and only if z; 4 0 for all i, see 
Figure 3.7. For given x in a fixed compact subset D of Q(z) there is a 
q, 0<q<1, such that 


|x;| < q|z;| for all 2. 


Set 


fa@= S- Cae. 
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— QZ) 


| Z 


Figure 3.7: Definition of D € Q(z) 


Proposition Al. (i) In every compact subset D of Q(z) one has f € 
C®™(D) and the formal differentiate series, that is )~, D'cqx%, is uniformly 
convergent on the closure of D and is equal to D® f. 


(ii) 
|D° f(x)| < M||!r—4l in D, 


where 


r = (1—q) min |2(|. 


Proof. See F. John [10], p. 64. Or an exercise. Hint: Use formula 12. where 
x is replaced by (q,...,q). 


Remark. From the proposition above it follows 


2,2 = D*s (0). 


Definition. Assume f is defined on a domain 2 C R”, then f is said to be 
real analytic in y € Q if there are ca € R and if there is a neighbourhood 
N(y) of y such that 


f(a) = J cal — y) 
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for all x € N(y), and the series converges (absolutely) for each « € N(y). 
A function f is called real analytic in Q if it is real analytic for each y € 2. 
We will write f € C’(Q) in the case that f is real analytic in the domain 2. 
A vector valued function f(x) = (fi(x),..., fm) is called real analytic if 
each coordinate is real analytic. 


Proposition A2. (i) Let fe C*(Q). Then f ¢ C°(Q). 


(ii) Assume f € C’(Q). Then for each y € Q there exists a neighbourhood 
N(y) and positive constants M, r such that 


Fle) =I (D*F()\(e - w)" 


Qa 


for all x € N(y), and the series converges (absolutely) for each x € N(y), 
and 
|D® F(a)| < M[Bltr"41. 


The proof follows from Proposition Al. 


An open set 2 € R” is called connected if Q is not a union of two nonempty 
open sets with empty intersection. An open set 2 € R” is connected if and 
only if its path connected, see [11], pp. 38, for example. We say that Q 
is path connected if for any x, y € © there is a continuous curve y(t) € Q, 
0<t< 1, with y(0) = # and 7(1) = y. From the theory of one complex 
variable we know that a continuation of an analytic function is uniquely 
determined. The same is true for real analytic functions. 


Proposition A3. Assume f € C%(Q) and (Q is connected. Then f is 
uniquely determined if for one z € Q all D° f(z) are known. 


Proof. See F. John [10], p. 65. Suppose g, h € C¥(Q) and D%g(z) = D*h(z) 
for every a. Set f =g—h and 

Q, = {xE€0Q: D* f(z) =0 for all a}, 

Q. = {xE: D° f(x) £0 for at least one a}. 


The set Q2 is open since D° f are continuous in 2. The set 9, is also open 
since f(a) = 0 in a neighbourhood of y € 0). This follows from 


Fe) = (DFW\(e- wy). 


Qa 


3.5. THEOREM OF CAUCHY-KOVALEVSKAYA 97 


Since z € 1, i. e., OQ, 4D, it follows N2 = 0. 


It was shown in Proposition A2 that derivatives of a real analytic function 
satisfy estimates. On the other hand it follows, see the next proposition, 
that a function f € C™ is real analytic if these estimates are satisfied. 


Definition. Let y € 2 and M, r positive constants. Then f is said to be 
in the class Cyz,-(y) if f € C® in a neighbourhood of y and if 


|D° f(y)| < M|B|br“? 
for all £. 


Proposition A4. f € C%(Q) if and only if f € C™®(Q) and for every 
compact subset S C Q there are positive constants M, r such that 


f€Currly) for ally € S. 


Proof. See F. John [10], pp. 65-66. We will prove the local version of the 
proposition, that is, we show it for each fixed y € 2. The general version 
follows from Heine-Borel theorem. Because of Proposition A3 it remains to 
show that the Taylor series 


Y pOWe-w 


Qa 


converges (absolutely) in a neighbourhood of y and that this series is equal 


to. f(z). 
Define a neighbourhood of y by 


Nay) = {x EQ: |av—ylt...+ en —- Yyn| < dh, 


where d is a sufficiently small positive constant. Set ®(t) = f(y+t(x— y)). 
The one-dimensional Taylor theorem says 


j-l 
fie) =e) => 20) + r;, 
k=0 


where 
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From formula 7. for directional derivatives it follows for x € Na(y) that 


i aot) = YS D*Fly t+ te —w)(e—w)* 


lol=j 


From the assumption and the multinomial formula 3. we get for0 <t< 1 


1 d |a|! —|a| Ay 
a2] < M gr |(z- y)*| 
lol=J 
= Mr¥? (lt. — yi] +... + [en — onl)? 


IA 
= 
oN 
31a 
NS 


Choose d > 0 such that d <r, then the Taylor series converges (absolutely) 
in Na(y) and it is equal to f(x) since the remainder satisfies, see the above 


estimate, 
1 1 sates j 
= a | (1 —t)?-'@I(t) dt| <M le 
Ge Dl ga r 


We recall that the notation f << F (f is majorized by F’) was defined in 
the previous section. 


Ir5| 


Proposition A5. (i) f = (fi,.-.,fm) € Cur(0) if and only if f << 
(®,...,®), where 


(ii) f € Cu,r(0) and f(0) = 0 if and only if 
f << (®-M,...,®-—M), 
where 


M(a1+...+2n) 
r—-21,—...— ky 


O(x) = 


Proof. 
D° (0) = Mal!r!l. 
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Remark. The definition of f << F implies, trivially, that D°f << D°F. 


The next proposition shows that compositions majorize if the involved func- 
tions majorize. More precisely, we have 


Proposition A6. Let f, F: R” +> R™ and g, G maps a neighbourhood of 
0 € R™ into R?. Assume all functions f(x), F(x), g(u), G(u) are in C™, 
f(0) = F(O) =0, f << F andg << G. Then g(f(x)) << G(F(2)). 


Proof. See F. John [10], p. 68. Set 

h(x) = g(f(x)), H(a) = G(F(@)). 
For each coordinate hy of h we have, according to the chain rule, 

D*hy(0) = Po(6°g(0), D7 F;(0)), 
where Py are polynomials with nonnegative integers as coefficients, P, are 
independent on g or f and 6 := (0/0uj,...,0/Oum). Thus, 
|D“hx(0)| P.(|5"91(0)|, |D F;(0)|) 

P.,(5°G(0), D’F;(0)) 
= D°H,(0). 


< 
< 


Using this result and Proposition A4, which characterizes real analytic func- 
tions, it follows that compositions of real analytic functions are real analytic 
functions again. 


Proposition A7. Assume f(x) and g(u) are real analytic, then g(f(a)) is 
real analytic if f(x) is in the domain of definition of g. 


Proof. See F. John [10], p. 68. Assume that f maps a neighbourhood of 
y € R” in R™ and g maps a neighbourhood of v = f(y) in R™. Then 
f € Curry) and g € C,,,p(v) implies 

h(x) = g(f (x)) € Cu,pr /(mM +p) (y). 


Once one has shown this inclusion, the proposition follows from Proposi- 
tion A4. To show the inclusion, we set 


h(y +z) :=g9(f(yt+2)) = o(vt+ flyt 2) — f(z)) = 9° (f*(2)), 
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where v = f(y) and 


g*(u): = gv+u)€C,,(0) 
Ff(y+2)— f(y) € Cu,r(0). 


In the above formulas v, y are considered as fixed parameters. From Propo- 
sition A5 it follows 


where 
M 
B(x) = - 
r Ly v2 aia Ln 
U(u) = Lp 
Pp— 2X1 — x2 In 


From Proposition A6 we get 


h(y +t) << (x(@),---+x(#)) = G(F), 


where 
2 Lp 
@) = 5 a@@)—M) 
_ pp(r — 21 —...— 2p) 
pr—(p+mM)(21+...+2n) 
ypr 
pr—(p+mM)(a1+...+ 2p) 
upr/(p+mM) 


pr/(p+mM) — (a1 +..-2n) 


See an exercise for the ”<<”-inequality. 
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Exercises 


. Let x: R” — R in C!, Vx 40. Show that for given xq € R” there 


is in a neighbourhood of xo a local diffeomorphism \ = ©®(x), ® : 
(@1,..-,2n) #2 (A1,..-,An), such that Ap = y(z). 


. Show that the differential equation 


a(x, Y)Ure + 2b(x, y)Ury + c(@, y)Uyy + lower order terms = 0 


is elliptic if ac — b? > 0, parabolic if ac — b? = 0 and hyperbolic if 
ac—b? <0. 


. Show that in the hyperbolic case there exists a solution of ¢; + pidy = 


0, see equation (3.9), such that Vo 4 0. 


Hint: Consider an appropriate Cauchy initial value problem. 


. Show equation (3.4). 


. Find the type of 


Lu := 2tgy + 2Ugy + 2Uyy = 0 


and transform this equation into an equation with vanishing mixed 
derivatives by using the orthogonal mapping (transform to principal 
axis) z = Uy, U orthogonal. 


. Determine the type of the following equation at (x,y) = (1,1/2). 


Lu := @Ugg + 2YUgy + 2ZyUyy = O. 


. Find all C?-solutions of 


Ure — 4Ugy + Uyy = 0. 


Hint: Transform to principal axis and stretching of axis lead to the 
wave equation. 


. Oscillations of a beam are described by 
1 
Wy — Rt =~ 0 


I 
2S 


On — pPUy 
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where o stresses, w deflection of the beam and E, p are positive con- 
stants. 

a) Determine the type of the system. 

b) Transform the system into two uncoupled equations, that is, w, o 
occur only in one equation, respectively. 

c) Find non-zero solutions. 


9. Find nontrivial solutions (Vy 4 0) of the characteristic equation to 
tUee — Uyy = f(x,y, u, Vu), 
where f is given. 
10. Determine the type of 
Ung — LUyg + Uyy + 3Ugz = 22, 
where u = u(x, y). 
11. Transform equation 
Ure + (1 y?)tay = 0, 
u = u(2,y), into its normal form. 
12. Transform the Tricomi-equation 
Yee + Uyy = 9, 
u=u(2,y), where y < 0, into its normal form. 


13. Transform equation 
tUre — YyUyy =0, 


u = u(x, y), into its normal form. 


14. Show that 
1 1 


(A+ [py (+ Ip?) 
are the minimum and maximum of eigenvalues of the matrix (a’), 


where 
( [pI ) ay 1+ [pl 


15. Show that Maxwell equations are a hyperbolic system. 
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16. 


17. 


18. 


19. 


20. 


21. 


22: 


23. 


24. 
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Consider Maxwell equations and prove that div E = 0 and div H = 0 
for all t if these equations are satisfied for a fixed time to. 


Hint. div rot A = 0 for each C?-vector field A = (Ai, Ag, A3). 


Assume a characteristic surface S(t) in R® is defined by x(a, y, z,t) = 
const. such that x; = 0 and x, 4 0. Show that S(t) has a nonparamet- 
ric representation z = u(x, y,t) with uz = 0, that is S(t) is independent 
of t. 

Prove formula (3.22) for the normal on a surface. 

Prove formula (3.23) for the speed of the surface S(t). 

Write the Navier-Stokes system as a system of type (3.24). 


Show that the following system (linear elasticity, stationary case of (3.25) 
in the two dimensional case) is elliptic 


pAut (A+) grad(div u) + f =0, 


where u = (w1,u2). The vector f = (fi, f2) is given and \, mw are 
positive constants. 


Discuss the type of the following system in stationary gas dynamics 
(isentrop flow) in R?. 


pPuug + pty + apr = 0 
PUVg + PUy + a’ py = 9 
plus + vy) + Upset upy = 0. 


Here are (u,v) velocity vector, p density and a = \/p’(p) the sound 
velocity. 


Show formula 7. (directional derivative). 


Hint: Induction with respect to m. 


Let y = y(x) be the solution of: 


y'(z) = f(x,y(x)) 
y(o) Yos 
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25. 


26. 


27. 


28. 


29. 
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where f is real analytic in a neighbourhood of (xo, yo) € R*. Find the 
polynomial P of degree 2 such that 

y(o) = P(w — a9) + O(a — 20)) 
as L—> Xo. 
Let u be the solution of 

Au = 1 
UELO)) Sake Oy 0: 
Find the polynomial P of degree 2 such that 
u(a,y) = P(w,y) + O((#? + y?)*””) 

as (x,y) — (0,0). 
Solve the Cauchy initial value problem 

y= ST (lt N-1%,) 

V(s,0) = 0. 

Hint: Multiply the differential equation with (r — s — NV). 
Write A?u = —u as a system of first order. 
Hint: A?u = A(Au). 
Write the minimal surface equation 


6) Us 6) Uy 


T 
Ox 1+ u2 a2 Oy fy bee ue 


as a system of first order. 


Hint: vy := Ue] 4/1 + u2 + ud, ve := uy/,/1+ u2 + ub. 


Let f: Rx R™” — R” be real analytic in (29, yo). Show that a real 
analytic solution in a neighbourhood of xo of the problem 


y' (a) f(x,y) 
y(Xo) Yo 


exists and is equal to the unique C![xo — €, xo + €|-solution from the 
Picard-Lindelof theorem, € > 0 sufficiently small. 
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105 
30. Show (see the proof of Proposition A7) 
pp(r— 21 —...— 2p) we pr 
pr—(p+mM)(a1+...+2n) pr —(p+mM)(a1+...+2n) 


Hint: Leibniz’s rule. 


Chapter 4 


Hyperbolic equations 


Here we consider hyperbolic equations of second order, mainly wave equa- 
tions. 


4.1 One-dimensional wave equation 


The one-dimensional wave equation is given by 


1 
ptt — Uae = 0, (4.1) 


where u = u(x,t) is a scalar function of two variables and c is a positive 
constant. According to previous considerations, all C?-solutions of the wave 
equation are 


u(x,t) = f(a + ct) + g(x — ct), (4.2) 
with arbitrary C?-functions f and g 


The Cauchy initial value problem for the wave equation is to find a C?- 
solution of 


1 

are —Ure = 0 
u(z,0) = a(x) 
u(z,0) = G(a), 


where a, 3 € C?(—00, 00) are given. 
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Theorem 4.1. There exists a unique C?(R x R)-solution of the Cauchy 
initial value problem, and this solution is given by d’Alembert’s' formula 


u(a,t) = 


a(e+ct)t+a(x—ct) 1 sete 
5 5c | B(s) ds. (4.3) 


—ct 


Proof. Assume there is a solution u(2,t) of the Cauchy initial value problem, 
then it follows from (4.2) that 


u(z,0) = f(x) +9(x) = a(z) 
cf’ (x) — cg'(x) = B(z). (4.5) 


a 
ey 
wo 


ut(x, 0) 


From (4.4) we obtain 
f(a) + g(x) = (2), 
which implies, together with (4.5), that 


a) = Led+ ere 
fia) — Le)= Blade 
Then 
fa) = 4S foe) astro, 
ge) = SO _ > [ as) as + Ce 


The constants C1, C2 satisfy 
Ci + C2 = f(x) + g(x) — a(x) = 0, 


see (4.4). Thus each C?-solution of the Cauchy initial value problem is given 
by d’Alembert’s formula. On the other hand, the function u(a,t) defined by 
the right hand side of (4.3) is a solution of the initial value problem. 


Corollaries. 1. The solution u(x,t) of the initial value problem depends 
on the values of a at the endpoints of the interval [x — ct, x + ct] and on the 
values of 3 on this interval only, see Figure 4.1. The interval [x — ct, x + ct] 
is called domain of dependence. 


'd’Alembert, Jean Babtiste le Rond, 1717-1783 
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(Xg,to) 


x+ct=const. 


J 


Xo—Ch Xo Xo +cty x 


Figure 4.1: Interval of dependence 


2. Let P be a point on the x-axis. Then we ask which points (x,t) need 
values of a or @ at P in order to calculate u(x,t)? From the d’Alembert 
formula it follows that this domain is a cone, see Figure 4.2. This set is 
called domain of influence. 


t 


x—ct=const. 


Figure 4.2: Domain of influence 


4.2 Higher dimensions 


Set 


u=un,—-CAu, A=A,=07/dr}+...4+ 07/022, 
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and consider the initial value problem 


u = 0 inR”xR 4.6) 
ute OU) = fH) 4.7) 
ur(@,0) = g(x), (4.8) 


where f and g are given C?(R7)-functions. 
By using spherical means and the above d’Alembert formula we will 
derive a formula for the solution of this initial value problem. 


Method of spherical means 


Define the spherical mean for a C?-solution u(x,t) of the initial value prob- 


lem by 
1 
M(r,t) = iy Bes u(y,t) dSy, (4.9) 


where 
wy = (2n)"/? /T(n/2) 


is the area of the n-dimensional sphere, wpr”~! is the area of a sphere with 
radius r. 

From the mean value theorem of the integral calculus we obtain the 
function u(x,t) for which we are looking at by 


u(a,t) = lim M(r,t). (4.10) 


Using the initial data, we have 


M(r,0) = Atl fly) dS, =: F(r) (4.11) 
M0) = Sot fy, WO y=), (412) 


which are the spherical means of f and g. 


The next step is to derive a partial differential equation for the spherical 
mean. From definition (4.9) of the spherical mean we obtain, after the 
mapping € = (y — x)/r, where x and r are fixed, 


1 
M(r,t) = — | u(x+r€,t) dS. 
By (0) 


Wn, 
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It follows 


M,(r,t) (0) 
i=1 


7 =e, Say, (y, t)&i Ay. 


B,.(«) j=1 


| . 

— yy, (a + r€, t)&: dS 

wn Jap, yi ( €,t)g § 
1 


Integration by parts yields 
1 n 
ed Uyiy:(y, t) dy 
Wr” : es X 


since € = (y—2)/r is the exterior normal at 0B,(a). Assume uw is a solution 
of the wave equation, then 


1 


ro M, = i utt(y,t) dy 
CoWn B,(«) 


1 r 
—. ——— t d . 
2a, | ee ae 


The previous equation follows by using spherical coordinates. Consequently 


1 
(r°-1 My)» = , ur(y,t) dSy 
dB, (a) 


Cw, 


rr! 6 1 
Cc ot? (si ss u(y, t) is, 


n-1 
r 
= Mit. 
C2 


Thus we arrive at the differential equation 
(r”-1M,.)> = cor Mi, 


which can be written as 


n—1l 


Mpr + M, = c 7 Mu. (4.13) 


This equation (4.13) is called Euler-Poisson-Darboux equation. 
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4.2.1 Case n=3 


The Euler-Poisson-Darboux equation in this case is 
(rM) pp = c 7 (rM)it. 


Thus rM is the solution of the one-dimensional wave equation with initial 
data 


(rM)(r,0) =rF(r) (rM)(r,0) = rG(r). (4.14) 
From the d’Alembert formula we get formally 


(r+ ct)F(r + ct) + (r — ct) F(r — ct) 
2r 


r+ct 
+ : / EG(E) dé. (4.15) 


2cr r—ct 


M(r,t) = 


The right hand side of the previous formula is well defined if the domain 
of dependence [x — ct, xz + ct] is a subset of (0,00). We can extend F' and 
G to Fo and Go which are defined on (—co, 00) such that rFo and rGp are 
C?(R)-functions as follows. Set 


Fir) : r>0 
Fo(r) = f(z) : r=0 
F(-r) : r<0o 


The function Go(r) is given by the same definition where F' and f are re- 
placed by G and g, respectively. 


Lemma. rFo(r), rGo(r) € C?(R?). 


Proof. From definition of F(r) and G(r), r > 0, it follows from the mean 
value theorem 


Jim, FP) = f(x), lim, G(r) = g(a). 


Thus rFfo(r) and rGo(r) are C(R)-functions. These functions are also in 
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C1(R). This follows since Fy and Gp are in C'(R). We have, for example, 


q 1 . 
F'(r) = a - dS¢ 


i 
S 
t 
mi 


Peo fy (Ey be 


= — oe, ny aS, 
ie yj (2) 8B, (0) 7g GMs 


= 0. 


Then, rFo(r) and rGo(r) are in C?(R), provided F” and G” are bounded as 
r — +0. This property follows from 


ve --| fury; (@ + EEE; aASe. 
8B, (0) ~> YY; fp WOE 
Thus 
F" (+0) fu; x) [ nn; dS¢. 
a Begg 


We recall that f,g € C?(R?) by assumption. 


The solution of the above initial value problem, where F and G are replaced 
by Fo and Go, respectively, is 
(r + ct) Fo(r + ct) + (r — ct)Fo(r — ct) 

2r 


EGo(E) dé. 


Mo (r, t) 
1 r+ct 


26T Jr neg 


Since Fp and Go are even functions, we have 


ct—r 
J, Sole) ag =0, 


Thus 


(r + ct) Fo(r + ct) — (ct — r)Fo(ct — r) 


Mo(r, t) = Di 


ct+r 


cae &Go(E) dé, (4.16) 


2cr J ot—p 
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t a | 
r—ct ct-r ct4+r 


Figure 4.3: Changed domain of integration 
see Figure 4.3. For fixed t > 0 and 0 <r < ct it follows that Mo(r,t) is the 


solution of the initial value problem with given initially data (4.14) since 
Fo(s) = F(s), Go(s) = G(s) if s > 0. Since for fixed t > 0 


u(a,t) = lim Mo(r, t), 
r—0 
it follows from d’Hospital’s rule that 


u(x,t) = ctF’(ct) + F(ct) +tG(ct) 
= © (iF (ct)) + #G(el). 


Theorem 4.2. Assume f € C%(IR?) and g € C?(R*) are given. Then 
there exists a unique solution u € C?(R® x [0,co)) of the initial value prob- 
lem (4.6)-(4.7), where n = 3, and the solution is given by the Poisson’s 
formula 


1 Of1 
es Arc? Ot Gan fe) is 
1 


Proof. Above we have shown that a C?-solution is given by Poisson’s for- 
mula. Under the additional assumption f € C®? it follows from Poisson’s 
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formula that this formula defines a solution which is in C?, see F. John [10}, 
p. 129. 


Corollary. From Poisson’s formula we see that the domain of dependence 
for u(x, to) is the intersection of the cone defined by |y — x| = c|t — to| with 
the hyperplane defined by t = 0, see Figure 4.4 


t 


(x,t 0) 


7 ly—xl=cl tty | 


Figure 4.4: Domain of dependence, case n = 3 


4.2.2 Casen=2 


Consider the initial value problem 


Ure + Vyy = Cun (4.18) 
v(z,y,0) = f(x,y) (4.19) 
ve(x,y,0) = g(x,y), (4.20) 


where f € C?, g € C?. 

Using the formula for the solution of the three-dimensional initial value 
problem we will derive a formula for the two-dimensional case. The following 
consideration is called Hadamard’s method of decent. 

Let v(x, y,t) be a solution of (4.18)-(4.20), then 


u(x, yy a, t) = (a, Y, t) 
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is a solution of the three-dimensional initial value problem with initial data 
f(x,y), g(x,y), independent of z, since u satisfies (4.18)-(4.20). Hence, since 
u(x, y, z,t) = ula, y,0,t) + uz(2, y, 6z,t)z,0<6< 1, and u, = 0, we have 


v(x, y,t) = u(a, y, 0, t). 


Poisson’s formula in the three-dimensional case implies 


1 afl 
t) = elie d 
v(x, y, ) Arc? Ot ( les f(E,n) 5) 
1 i 
in) dS. 4.21 
ord ee g(§n) (4.21) 


E 


Figure 4.5: Domains of integration 


The integrands are independent on ¢. The surface S' is defined by 
x(€,7, 0 := (€-2)? + (n—-y)? + C — 7t? = 0. Then the exterior normal n 
at S isn = Vx/|Vx| and the surface element is given by dS = (1/|n3|)d&dn, 
where the third coordinate of n is 


fet — (€= 2)? —(m—y)? 
‘3 ct , 


ng = 


4.3. INHOMOGENEOUS EQUATION 117 


The positive sign applies on St, where ¢ > 0 and the sign is negative on S~ 
where ¢ < 0, see Figure 4.5. We have $= StUS-. 
Set p= /(€—2)? + (n — y)?. Then it follows from (4.21) 


Theorem 4.3. The solution of the Cauchy initial value problem (4.18)- 
(4.20) is given by 


ee ee f(é,n) 
v(z,y,t) = aces Jato d&dn 
1 g(&,) dédn. 


4+. | 
ANC JBaley) Vet — p? 


Corollary. In contrast to the three dimensional case, the domain of depen- 
dence is here the disk Be,(xo, yo) and not the boundary only. Therefore, see 
formula of Theorem 4.3, if f, g have supports in a compact domain D C R?, 
then these functions have influence on the value v(z, y, t) for all time t > T, 
T sufficiently large. 


4.3. Inhomogeneous equation 


Here we consider the initial value problem 


u = w(z,t) onzeR",teR (4.22) 
u(z,0) = f(x) (4.23) 
u(x,0) = g(x), (4.24) 


where Ou := uz — c2Au. We assume fe CG. gE C? and w € C!, which 
are given. 

Set wu = uz + ue, where uj is a solution of problem (4.22)-(4.24) with 
w := 0 and uz is the solution where f = 0 and g = 0 in (4.22)-(4.24). Since 
we have explicit solutions u; in the cases n = 1, n = 2 and n = 3, it remains 
to solve 


u = w(z,t) onze R”,teR (4.25) 
u(z,0) = 0 (4.26) 
ur(z,0) = 0. (4.27) 


The following method is called Duhamel’s principle which can be considered 
as a generalization of the method of variations of constants in the theory of 
ordinary differential equations. 
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To solve this problem, we make the ansatz 


t 
u(a,t) = | u(a,t,s) ds, (4.28) 
0 
where v is a function satisfying 
v=0 for all s (4.29) 
and 
v(x, s,s) =0. (4.30) 


From ansatz (4.28) and assumption (4.30) we get 
t 
ut = v(a,t,t) +f u(x, t, s) ds, 
0 


= iE uz(2, t, 8). (4.31) 
0 


It follows w:(x,0) = 0. The initial condition u(x,t) = 0 is satisfied because 
of the ansatz (4.28). From (4.31) and ansatz (4.28) we see that 


t 
Ut = w(e.tt) + f v(x, t, 8) ds, 
0 
t 
A;u = i A,v(a,t, s) ds. 
0 


Therefore, since u is an ansatz for (4.25)-(4.27), 


t 
un —CA,uU = wu(estst) + f( v) (x,t, s) ds 
0 


w(a,t). 


Thus necessarily vu;(x,t,t) = w(ax,t), see (4.29). We have seen that the 
ansatz provides a solution of (4.25)-(4.27) if for all s 


v=0, v(2,8,5)=0, u(x, 5,5) = w(z, 8). (4.32) 


Let v*(x,t,s) be a solution of 


v=0, v(x,0,s)=0, u(x,0,s) = w(z, 8), (4.33) 


then 
v(a, t, s) = v* (x,t — 8,8) 
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is a solution of (4.32). In the case n = 3, where v* is given by, see Theorem 
4.2, 
1 
t dSe. 
ae Ancet i a 


Then 


v(a,t,s) = v*(a,t— s,s) 
1 


——— ie ey aS 
4rc?(t — s) Vas (6,8) de 


from ansatz (4.28) it follows 


t 

h v(az, t, s) ds 

= w(§, 8) 

= dSeds. 
ine Bees) (a) . —s 


Changing variables by 7 = c(t — s) yields 


7 w(t 1/6) 
u(z,t) = in| Vea ie = dS¢dt 


2 Baa 
Arc? Bet (x) 
where r = |x — &|. 


Formulas for the cases n = 1 and n = 2 follow from formulas for the as- 
sociated homogeneous equation with inhomogeneous initial values for these 
cases. 


u(x,t) 


Theorem 4.4. The solution of 


u=w(2,t), u(#,0)=0, u(r,0) =0 


where w € C1, is given by: 


u(e,t) = 5 f 5 wet ae 


Arc? 


Case n = 3: 


where r = |x — &|, = (41, 22,73), € = (£1, €9, &3). 
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es er i w(E,T) 
ne = Ze | ( how i i) i 


1 t xz+c(t—T) 
u(x,t) = ma (es w(€,T) i) dt. 


Remark. The integrand on the right in formula for n = 3 is called retarded 
potential. The integrand is taken not at t, it is taken at an earlier time 
t—r/c. 


4.4 A method of Riemann 


Riemann’s method provides a formula for the solution of the following 
Cauchy initial value problem for a hyperbolic equation of second order in 
two variables. Let 


S71 x=x(t),y = y(t), ty St< to, 


be a regular curve in R?, that is, we assume 2, y € C1[ty, t2] and 2?+y”? # 0. 
Set 

Lu := Uny + a(@, y)ug + b(x, yy + c(z, y)u, 
where a, b€ C! andc, f € C ina neighbourhood of S. Consider the initial 
value problem 


Iu = f(x,y) (4.34) 
uo(t) = u(x(t), y(t)) (4.35) 
po(t) = ue(x(t), y(t)) (4.36) 
qo(t) = Uy(x(t), y(t), (4.37) 


where f € C in a neighbourhood of S and uo, po, go € C! are given. 


We assume: 
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(i) up(t) = po(t)a’ (t) + go(t)y'(t) (strip condition), 


(ii) S is not a characteristic curve. Moreover assume that the characteristic 
curves, which are lines here and are defined by x = const. and y = const., 
have at most one point of intersection with S, and such a point is not a 
touching point, i. e., tangents of the characteristic and S are different at 
this point. 


We recall that the characteristic equation to (4.34) is y2xy = 0 which 
is satisfied if yz(a,y) = 0 or xy(z,y) = 0. One family of characteris- 
tics associated to these first partial differential of first order is defined by 
z(t) =1, y'(t) = 0, see Chapter 2. 


Assume u, v € C! and that uzy, Vey exist and are continuous. Define the 
adjoint differential expression by 


Mv = Ury — (av) 2 — (bv)y + cv. 


We have 


2(vLu —uMv) = (ugv — veut 2buv)y + (uyv — vyu+ 2auv)e. (4.38) 
Set 


P = —(ugu— tzu + 2buv) 


Q = Uy — vyut 2auv. 
From (4.38) it follows for a domain Q € R? 
2 | (vLu—uMv) dxdy = | (—Py+ Qy) dady 
Q 2 
= ¢ Pdz + Qdy, (4.39) 


where integration in the line integral is anticlockwise. The previous equation 
follows from Gauss theorem or after integration by parts: 


| (—P) + Q,) dedy = | (—Pn2 + Qn1) ds, 
Q dQ 


where n = (dy/ds, —dx/ds), s arc length, (x(s), y(s)) represents 00. 
Assume u is a solution of the initial value problem (4.34)-(4.37) and 
suppose that v satisfies 
Mv=0 in. 


122 CHAPTER 4. HYPERBOLIC EQUATIONS 


P=(Xy.Vq) 


Xx 


Figure 4.6: Riemann’s method, domain of integration 


Then, if we integrate over a domain 2 as shown in Figure 4.6, it follows 
from (4.39) that 


2 | uf dedy = | Pax+Qdy+ [ Pax-+Qay+ [ Pdx+Qdy. (4.40) 
Q BA AP PB 


The line integral from B to A is known from initial data, see the definition 


of P and Q. 
Since 
UgV — VzU + 2buv = (uv), + 2u(bu — vz), 
it follows 
| Pdx+Qdy = -| ((uv)2 + 2u(bu — vz)) dx 
AP AP 


—(uv)(P) + (uv) (A) — i: 2u(bu — vz) da. 


By the same reasoning we obtain for the third line integral 


i. Pdz+Qdy = i, ((uv)y + 2u(av — vy)) dy 
PB PB 


= (uv)(B) — (uv)(P) +f 2u(av — vy) dy. 


PB 
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Combining these equations with (4.39), we get 
2u(P)u(P) = i; (uzv — Uy + 2buv) dx — (uyv — vyu+ 2auv) dy 
BA 


+u(A)v(A) + u(B)v(B) + 2 fu(be —v,) dx 


+2 i: u(av — vy) dy — 2 f fe dady. (4.41) 


Let v be a solution of the initial value problem, see Figure 4.7 for the defi- 
nition of domain D(P), 


y 


P=(X,.Yq) c 


D(P) 


Figure 4.7: Definition of Riemann’s function 


Mv = 0 in D(P) (4.42) 
byv—v, = 0 onCy (4.43) 
avu—vy = 0 onC, (4.44) 

wP). =: 41; (4.45) 


Assume v satisfies (4.42)-(4.45), then 
2u(P) = u(A)v(A) + u(B)v(B) - 2 | fu dxdy 
Q 


= | (UgY — Vz + 2buv) dx — (uyv — vyu + 2auv) dy, 
BA 
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where the right hand side is known from given data. 


A function v = v(2,y;20, yo) satisfying (4.42)-(4.45) is called Riemann’s 
function. 


Remark. Set w(z,y) = v(x, y; Xo, yo) for fixed xo, yo. Then (4.42)-(4.45) 
imply 


w(x, yo) exp ( b(7, yo) ar) on Cj, 


xO 


y 
w(to,y) = exp € a(Xo, T) ar) on C. 
yl 


0 


Examples 
1. Ury = f(x,y), then a Riemann function is v(x, y) = 1. 
2. Consider the telegraph equation of Chapter 3 
Epluy = Agu — Amur, 
where u stands for one coordinate of electric or magnetic field. Introducing 
u=w(z,te™, 
where & = —A/(2e), we arrive at 


2 2 
c 
Sp 
Wtt an aW 4e2 


Stretching the axis and transform the equation to the normal form we get 
finally the following equation, the new function is denoted by u and the new 
variables are denoted by x, y again, 


Ugy + cu =0, 
with a positive constant c. We make the ansatz for a Riemann function 
V(x, Y; Zo, Yo) = w(s), 8 = (x&—2x0)(y — Yo) 


and obtain 
sw’ +w'+cw =0. 
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Substitution o = 4cs leads to Bessel’s differential equation 
a7 2"(a) +. a2'(c) + 072(c) = 0, 


where z(a) = w(a?/(4c)). A solution is 


Jo(7) = Jo ( y/4e(a = x0)(y — yo)) 
which defines a Riemann function since Jo(0) = 1. 


Remark. Bessel’s differential equation is 


xy" (x) + ay’ (a) + («” — n*)y(a) = 0, 
where n € R. If n € NU {0}, then solutions are given by Bessel functions. 
One of the two linearly independent solutions is bounded at 0. This bounded 
solution is the Bessel function J,,(x) of first kind and of order n, see [1], for 
example. 


4.5 Initial-boundary value problems 


In previous sections we looked at solutions defined for all « € R" andt € R. 
In this and in the following section we seek solutions u(x,t) defined in a 
bounded domain 2 C R” and for all ¢ € R and which satisfy additional 
boundary conditions on 0. 


4.5.1 Oscillation of a string 


Let u(x, t), x € [a,b], t € R, be the deflection of a string, see Figure 1.4 from 
Chapter 1. Assume the deflection occurs in the (x, w)-plane. This problem 
is governed by the initial-boundary value problem 


tet, by = tee a yt) con (0,2) (4.46) 
GeO) SF (ey) (4.47) 
ur(z,0) = g(x) (4.48) 
a0,t) = at) =O: (4.49) 


Assume the initial data f, g are sufficiently regular. This implies compati- 
bility conditions f(0) = f(l) = 0 and g(0) = g(l). 
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Fourier’s method 


To find solutions of differential equation (4.46) we make the separation of 
variables ansatz 
u(x, t) = v(x) w(t). 


Inserting the ansatz into (4.46) we obtain 
v(x)w"'(t) = o"(x)w(t), 


or, if v(x)w(t) 4 0, 


It follows, provided v(a)w(t) is a solution of differential equation (4.46) and 


v(x)w(t) # 9, 
w" (t) 


w(t) 


= const. =: —X 


and 


since x, t are independent variables. 
Assume v(0) = v(l) = 0, then v(x) w(t) satisfies the boundary condi- 
tion (4.49). Thus we look for solutions of the eigenvalue problem 


=o" (2): =. Ava) att (0;1) (4.50) 
v(0) = v(l) =0, (4.51 


and associated eigenfunctions are 


_ (Tt 
Un = sin (Gna) : 
Solutions of 
—w" (t) = Anw(t) 


sin(/Ant), cos(+/ Ant). 


W(t) = ap cos(1/Ant) + Bp sin(/Ant), 


are 


Set 
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where @, By, € R. It is easily seen that w,,(t)v,(«) is a solution of differential 
equation (4.46), and, since (4.46) is linear and homogeneous, also (principle 
of superposition) 


N 
un = >) wn(t)on(2) 
n=1 


which satisfies the differential equation (4.46) and the boundary condi- 
tions (4.49). Consider the formal solution of (4.46), (4.49) 


u(a,t) = + (an cos(1/Ant) + Bn sin(/\nt)) sin (nv) : (4.52) 
n=1 


”Formal” means that we know here neither that the right hand side con- 
verges nor that it is a solution of the initial-boundary value problem. For- 
mally, the unknown coefficients can be calculated from initial conditions (4.47), 
(4.48) as follows. We have 


u(z,0) = y Op Sin(\/Anx) = f(a). 
n=1 


Multiplying this equation by sin(/A;,2) and integrate over (0,1), we get 


l l 
an | sin? (Anz) dx =| f(x) sin(/A,x) dx. 
0 0 


We recall that 


l 
i 
if sin(/Anx) sin(Ap2) dx = 3 onk- 
0 


ae = ° i, f(a) sin (%e) dic. (4.53) 


By the same argument it follows from 


Then 


ur(x, 0) = S- Bn An sin(/Anx) = g(x) 
n=1 


that 


D7 k 
Be= al g(x) sin (7c) dx. (4.54) 


Under additional assumptions f € C§(0,1), g € C§(0,1) it follows that 
the right hand side of (4.52), where an, By, are given by (4.53) and (4.54), 
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respectively, defines a classical solution of (4.46)-(4.49) since under these 
assumptions the series for u and the formal differentiate series for uz, wie, 
Ug, Ure Converges uniformly on0 <2 <1,0<t<T,0<T < o fixed, see 
an exercise. 


4.5.2 Oscillation of a membrane 


Let Q C R? be a bounded domain. We consider the initial-boundary value 
problem 


ux(z,t) = A,uwu nQxR, (4.55) 
u(aj0) = Fe), 2 EO, (4.56) 
u(z,0) = g(x), rEQ, (4.57) 
u(z,t) = 0 ondQxR. (4.58) 


As in the previous subsection for the string, we make the ansatz (separation 
of variables) 

u(x,t) = w(t)v(2) 
which leads to the eigenvalue problem 


—Av = dM inQ, (4.59) 

v = 0 ondQ. (4.60) 

Let Ap are the eigenvalues of (4.59), (4.60) and v, a complete associated 
orthonormal system of eigenfunctions. We assume 2 is sufficiently regular 
such that the eigenvalues are countable, which is satisfied in the following 


examples. Then the formal solution of the above initial-boundary value 
problem is 


Ce S- (an cos(/Ant) + Bn sin(/\nt)) arvenge 


n=1 


where 
Qn = | f(x)vn(a) da 
Be = S| Been) de. 
Remark. In general, eigenvalues of (4.59), (4.59) are not known explicitly. 


There are numerical methods to calculate these values. In some special 
cases, see next examples, these values are known. 
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Examples 
1. Rectangle membrane. Let 
Q = (0,a) x (0,0). 


Using the method of separation of variables, we find all eigenvalues of (4.59), (4.60) 
which are given by 


2 
Akl = k,l=1,2,... 
a 


and associated eigenfunctions, not normalized, are 


nie wk : al 
ug (2) = sin 7 sin | 5-22 } . 


2. Disk membrane. Set 
Q={reR’?: si+2? < R71. 


In polar coordinates, the eigenvalue problem (4.59), (4.60) is given by 


= ((rurde + uo = du (4.61) 
u(R,6) = 0, (4.62) 


here is u = u(r,6) := v(rcos6,rsin@). We will find eigenvalues and eigen- 
functions by separation of variables 


u(r, 0) = v(r)q(8), 


where v(R) = 0 and q(@) is periodic with period 2m since u(r, @) is single 
valued. This leads to 


Dividing by vq, provided vq 4 0, we obtain 


1 ((rv'(r)) 1g") 
a ots we) ) = cee 


which implies 
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Thus, we arrive at the eigenvalue problem 


Lq(@) 
q(@ + 2m). 


| 
iQ 
Q oS 
, erhii S ee aa 
ss 
Sa 
| | 


It follows that eigenvalues 4 are real and nonnegative. All solutions of the 
differential equation are given by 


q(@) = Asin(,/“@) + B cos(/18), 
where A, B are arbitrary real constants. From the periodicity requirement 
Asin(./10) + Bcos(,/u6) = Asin(./u(6 + 27)) + Bcos(/pu(6 + 27)) 
it follows? 
sin(./p7) (A cos(./u8 + fun) — Bsin(,/u8 + /um)) = 0, 


which implies, since A, B are not zero simultaneously, because we are look- 
ing for q not identically zero, 


sin(,/pi77) sin(.,/w8 + 6) = 0 
for all @ and a 6 = 6(A, B, 4). Consequently the eigenvalues are 
Ln =n’, T= Og lyese% 


Inserting q"(8)/q(0) = —n? into (4.63), we obtain the boundary value prob- 
lem 


ru" (r) + rv (r) + (Ar? —n?)v = 0 on (0,R) (4.64) 
UCR) =. 30 (4.65) 
sup |v(r)| < o. (4.66) 

ré(0,R) 


sinzgz—siny = 2 cos 2 sin 


- & : 
cosx — COSY = 2sin 2 sin 
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where z > 0. Solutions of this differential equations which are bounded at 
zero are Bessel functions of first kind and n-th order J,,(z). The eigenvalues 
follows from boundary condition (4.65), i. e., from Jn(WAR) = 0. Denote 
by Tz the zeros of J,(z), then the eigenvalues of (4.61)-(4.61) are 


sa = (8) 


and the associated eigenfunctions are 
In(V Ankr) sin(nd), a be ere 
In(/ Ankr) cos(n8), n=0,1,2)..2 
Thus the eigenvalues Apx are simple and Ajz,, n > 1, are double eigenvalues. 


Remark. For tables with zeros of J,(a) and for much more properties of 
Bessel functions see [25]. One has, in particular, the asymptotic formula 


Rey= eu (costa nelD nis) LO (+)) 


as x — oo. It follows from this formula that there are infinitely many zeros 


Obs gla): 


4.5.3. Inhomogeneous wave equations 


Let 2 c R” be a bounded and sufficiently regular domain. In this section 
we consider the initial-boundary value problem 


ut: = Lut f(a,t) nQOxR (4.67) 

u(v,0) = (x) rEQ (4.68) 

u(t,0) = (2) cen (4.69) 

u(z,t) = 0 fore € OQ andte€ R”, (4.70) 

where u = u(x,t), © = (%1,...,%n), f, ¢, W are given and L is an elliptic 


differential operator. Examples for LD are: 
1. L = 6? /0x?, oscillating string. 


2. L = Ay,, oscillating membrane. 
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where aI = a!’ are given sufficiently regular functions defined on 2. We 
assume L is uniformly elliptic, that is, there is a constant v > 0 such that 


So a Ge; > IC)? 


t,j=1 


for all ¢ € QO and ¢ € R”. 


4, Let wu = (t1,...,Um) and 


n fa) P 
Lu = 2 Ape (A? (x) tz, ) 5 
ij=l 9 


where Av = AJ’ are given sufficiently regular (m x m)-matrices on Q. We 
assume that L defines an elliptic system. An example for this case is the 
linear elasticity. 


Consider the eigenvalue problem 


—-Iv = dv inQ (4.71) 
v = 0 on AQ. (4.72) 


Assume there are infinitely many eigenvalues 
0< AY <AV<... SOO 


and a system of associated eigenfunctions v1, v2,... which is complete and 
orthonormal in L?(Q). This assumption is satisfied if Q is bounded and if 
OQ is sufficiently regular. 

For the solution of (4.67)-(4.70) we make the ansatz 


foe) 


u(x,t) = S> vp(x)we(t), (4.73) 


k=1 
with functions w,;(t) which will be determined later. It is assumed that all 
series are convergent and that following calculations make sense. Let 


co 


f(z,t) = Ce(t) vp (2) (4.74) 
k=1 
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be Fourier’s decomposition of f with respect to the eigenfunctions vg. We 
have 


cp (t) =[ f(a, t)ug(x) da, (4.75) 


which follows from (4.74) after multiplying with v;(2) and integrating over 
Q. 


Set 
(6,04) = fA baoglae) de, 
then 
dc) = S°(d,ve)ve(x) 
k=1 


= 
& 

I 
Rage: 


(), Ue)R(2) 


are Fourier’s decomposition of ¢ and w, respectively. 
In the following we will determine w;(t), which occurs in ansatz (4.73), 
from the requirement that u = v;,(a7)w,(t) is a solution of 


ute = Lu + cp (t)vg(x) 
and that the initial conditions 
we(0) = (9, ve), wh,(O) = (eb, ve) 
are satisfied. From the above differential equation it follows 


wz (t) = —rpwp(t) + cr(t). 


Thus 
wre(t) = ap cos(/Agt) + by sin(V/Agt) (4.76) 
1 u . 
+ | Ce(T) sin(V/Ag(t —7)) dr, 
where i 
ap = (, UK), be = IW, abana 


Summarizing, we have 
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Proposition 4.2. The (formal) solution of the initial-boundary value prob- 
lem (4.67)-(4.70) is given by 


foe) 


u(x,t) = S> vp(x)we(t), 


k=1 
where vz is a complete orthonormal system of eigenfunctions of (4.71), (4.72) 
and the functions wz are defined by (4.76). 
The resonance phenomenon 
Set in (4.67)-(4.70) ¢ = 0, ~ = 0 and assume that the external force f is 
periodic and is given by 


f(a, t) = Asin(wt)un (a), 


where A, w are real constants and v,, is one of the eigenfunctions of (4.71), (4.72). 
It follows 


Celt) | f(a, t)ug(x) dx = Adjx sin(wt). 
2 
Then the solution of the initial value problem (4.67)-(4.70) is 


= Avn(z) an wT) sin —T)) dr 
u(x,t) = a (wr) sin(\/Aa(t—7)) a 


Aon Oy (+ sin(/ut) — sin(t) | 


provided w 4 VAn. It follows 


A sin(V/Ant) 
u(x,t) > DJign(*) (Me = rw) 


if ww — VA. The right hand side is also the solution of the initial-boundary 
value problem if w = V/Xn.- 
Consequently |u| can be arbitrarily large at some points x and at some 


times t if w = WAn. The frequencies VA, are called critical frequencies at 
which resonance occurs. 


A uniqueness result 


The solution of of the initial-boundary value problem (4.67)-(4.70) is unique 
in the class C?(Q x R). 
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Proof. Let uz, ug are two solutions, then u = ug — uj, satisfies 


ut = Lu nQxR 

0) = 0 rE 

0) = Owe 0 

t) = 0 forx€ OQ andt eR”. 


As an example we consider Example 3 from above and set 


n 


E(t) -[ (s” a’) (x)Ue,Ue,; + ugut) dx. 


tj=l 


Then 


n 


Bey! = 2f (Se a’) (x)Uy;Uejt + Ueuet) da 
ij=l 
n 


= 2 | (S” a) (x)ue,uin;) dS 
0Q ij=l 
+2 f uz(—Lu + uzt) dx 
Q 
= 0. 
It follows E(t) = const. From u(x,0) = 0 and u(x,0) = 0 we get E(0) = 0. 
Consequently E(t) = 0 for all t, which implies, since L is elliptic, that 


u(x,t) = const. on Q x R. Finally, the homogeneous initial and boundary 
value conditions lead to u(z,t) =0 on 2 x R. 
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4.6 Exercises 


1. Show that u(a,t) € C?(R?) is a solution of the one-dimensional wave 
equation 
Ute = CUee 


if and only if 
u(A) + u(C) = u(B) + u(D) 


holds for all parallelograms ABC D in the (x, t)-plane, which are bounded 
by characteristic lines, see Figure 4.8. 


t 


Figure 4.8: Figure to the exercise 


2. Method of separation of variables: Let vz(x) be an eigenfunction to 
the eigenvalue of the eigenvalue problem —v"(x) = Av(x) in (0,1), 
v(0) = v(l) = 0 and let w,(t) be a solution of differential equation 
—w"(t) = A,w(t). Prove that vz(x)wz(t) is a solution of the partial 
differential equation (wave equation) uz = Ure. 


3. Solve for given f(a) and yz € R the initial value problem 


Upt+Ug + UUereere = O in Rx Ry 
u(z,0) = f(z). 


4. Let S := {(a,t); t = yx} be spacelike, i. e., |y| < 1/c?) in (z, t)-space, 
x = (%1,%2,73). Show that the Cauchy initial value problem Ou = 0 
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with data for u on S can be transformed using the Lorentz-transform 


x1 — crt 
TP 


into the initial value problem, in new coordinates, 


t— yx 


ia 


/ / / 
» ty = %Q, C3 = 13, 6 = 


t= 


u = 0 
u(a’,0) = f(a’) 
uy(2’,0) = g(a’). 


Here we denote the transformed function by u again. 


5. (i) Show that 
as THY gin (7 
u(x,t) : 2 on c08( i t) sin ( i r) 


is a C?-solution of the wave equation ut = Ure if |r| < c/n, where 
the constant c is independent of n. 


(ii) Set 
; ™m 
Qn i= | f(x)sin( —ax) dz. 
f Seen (Fe) 
Prove |an| < c/n*, provided f € C§(0, 1). 


6. Let 2 be the rectangle (0,a) x (0,0). Find all eigenvalues and associ- 
ated eigenfunctions of —Au = Au in Q, u=0 on AN. 


Hint: Separation of variables. 


7. Find a solution of Schrédinger’s equation 
he 
thu, = ———Azw +V(x2)y in R” xR, 
2m 
which satisfies the side condition 
n 
jf weePae=1, 
R 
provided FE € R is an (eigenvalue) of the elliptic equation 


2 
Aut so (E —~V(z))u=0 in R” 
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10. 


11. 


12. 


CHAPTER 4. HYPERBOLIC EQUATIONS 


under the side condition fy |u|?dz =1,u: R°% C. 
Here is w : R” x Rt C, h Planck’s constant (a small positive con- 
stant), V(x) a given potential. 


Remark. In the case of a hydrogen atom the potential is V(x) = 
—e/|x|, e is here a positive constant. Then eigenvalues are given by 
En = —me*/(2h?n?), n € N, see [22], pp. 202. 


. Find nonzero solutions by using separation of variables of uz = A,u 


in Q x (0,00), u(x,t) = 0 on OO, where C is the circular cylinder 
0 = {(zi, 29,23) € R” > at +02 < R?, 0 < ag < A}. 


. Solve the initial value problem 


3Utt = Athn:55 = 0 
u(z,0) = sing 
u(xv,0) = 1. 


Solve the initial value problem 


Utt — Cure = x’, t>0, cER 
u(z,0) = 2 
u(z,0) = 0. 


Hint: Find a solution of the differential equation independent on t, 
and transform the above problem into an initial value problem with 
homogeneous differential equation by using this solution. 


Find with the method of separation of variables nonzero solutions 
u(x,t), O<a<1,0<t<ow, of 


Ut — Ure FU=0, 
such that u(0,t) = 0, and u(1,t) = 0 for all t € [0, 00). 
Find solutions of the equation 
un — CUre = Mu, A = const. 
which can be written as 


u(z,t) = f(x? —c*#”) = f(s), s:= a? — ct? 


4.6. 


13. 


14. 
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with f(0) = K, K a constant. 


Hint: Transform equation for f(s) by using the substitution s := z?/A 
with an appropriate constant A into Bessel’s differential equation 


2 f"(@2)+2f' 2) +2? =r )\f ao, 250 


with n = 0. 


Remark. The above differential equation for u is the transformed tele- 
graph equation (see Section 4.4). 


Find the formula for the solution of the following Cauchy initial value 
problem uzy = f(x,y), where S: y = ax +b, a > 0, and the initial 
conditions on S are given by 


u = axt By+7, 
Uz = a, 
Uy = 8B, 

a, b, a, @, y constants. 

Find all eigenvalues ps of 


—q"(0) = pq() 
q(@) = q(0+2z). 


Chapter 5 


Fourier transform 


Fourier’s transform is an integral transform which can simplify investigations 
for linear differential or integral equations since it transforms a differential 
operator into an algebraic equation. 


5.1 Definition, properties 


Definition. Let f € Cj(R"), s=0,1,.... The function f defined by 


Fl) = (em)? | efx) der, (5.1) 


n 


where € € R”, is called Fourier transform of f, and the function g given by 


Gar) = (2n)-"? | elf 9(€) a€ (5.2) 


n 


is called inverse Fourier transform, provided the integrals on the right hand 
side exist. 


From (5.1) it follows by integration by parts that differentiation of a func- 
tion is changed to multiplication of its Fourier transforms, or an analytical 
operation is converted into an algebraic operation. More precisely, we have 


Proposition 5.1. 
Do f(g) = dle F(), 


where |a| < s. 
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The following proposition shows that the Fourier transform of f decreases 
rapidly for || — oo, provided f € C§(IR”). In particular, the right hand 
side of (5.2) exists for g:= f if f € OpT(R®), 


Proposition 5.2. Assume g € C)(IR”), then there is a constant M = 


M(n, s,g) such that 
M 


Proof. Let € = (&1,...,&n) be fixed and let 7 be an index such that |&;| = 
max, ||. Then 


4 1/2 
eS (>: a) < Vnilé;| 
k=1 


which implies 


(+e) = 3 Ql 


8 
Ss ey intel 
k=0 
Z 2878/2 Sy |E*. 
la|<s 


This inequality and Proposition 5.1 imply 


(L+IE))*1HE)] < 2%n*? SY |(HE)°G(O)| 


la|<s 


yr S- | |D%g(ax)| dx =: M. 
R” 


lja|<s 


IA 


The notation inverse Fourier transform for (5.2) is justified by 
Theorem 5.1. f= f and f= f. 
Proof. See [27], for example. We will prove the first assertion 


(ony? fb Fle) dé = f) (5.3 
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here. The proof of the other relation is left as an exercise. All integrals 
appearing in the following exist, see Proposition 5.2 and the special choice 
of g. 


(i) Formula 


[, sfc’ ac =f auyre+y) ay (65.4) 


Rn 


follows by direct calculation: 


ff wo (emr? fern dy) o€ ag 


= (Qn)? ‘a ; ( i age? a) f(y) dy 


= [Gu 2)Fy) ay 
= [dwt +y) ay 


(ii) Formula 


(anynl? fo WEg(e6) dé = e-"Glu/e) (55) 
for each ¢ > 0 follows after substitution z = e€ in the left hand side of (5.1). 


iii) Equation 
(iii) Eq 


[, sepFee ae= fi auto + ev) ay (5.6) 


follows from (5.4) and (5.5). Set G(€) := g(e€), then (5.4) implies 


n~ 


G()fee§ d= | Gly) f(x+y) dy. 
R” R” 


Since, see (5.5), 
(onyn? fee g(ee) dé 


e "G(y/e), 


2 
= 
| 
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we arrive at 


[, sOF® = 


Letting « — 0, we get 


Fle) > dé = flax g ; : 
a(0) {Flee ag = fle) fat) ay (5.7) 
Set 
g(2) = ele? /2 
then 
[5 ay = em? (5.8) 


Since g(0) = 1, the first assertion of Theorem 5.1 follows from (5.7) and (5.8). 
It remains to show (5.8). 


(iv) Proof of (5.8). We will show 


Gly): = (2n)-/? i elt /2.-iee ge 
ewe. 
The proof of 
[ el"? dy = (2m)"/? 


is left as an exercise. Since 
et elle 
2 2 


(a | 5) ) (Fa | 5) - 


i PA ee ped Se 2 
i. el? /29-iY de = i an? sa WP/2 ay 
n nm 


—lay|2 “2 
ce chin | oda 


= anag lohe f eT dy 


it follows 
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where 
x .Y 
) = + .—_. 


V2 V2 


Consider first the one-dimensional case. According to Cauchy’s theorem we 


have 
¢ er dn = 0, 
C 


where the integration is along the curve C' which is the union of four curves 
as indicated in Figure 5.1. 


Im n 
ix 
Bo 
Cc, C 4 
-R e R Ren 
Figure 5.1: Proof of (5.8) 
Consequently 


2 1 R 2 2. 2 
e” dj= = / e* /? dx -{ e” dy -{ e” dn. 
[. V2 —R C2 C4 


It follows 
lim er dn = Jn 
R-0oo C3 

since 


lim e? dn=0, k=2, 4. 
R—-0oo Cr 


The case n > 1 can be reduced to the one-dimensional case as follows. Set 


— +i4 =( ) 
SS FS = pees ‘ 
7 V2 V2 mM Tn 


where 
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From dn = dn, ...dm and 


n 
2 2. 2 
et —e hia — | ent 
i=1 


n 
| eo” dn = IL / et dm, 
- =r 7 


it follows 


where for fixed y 


zx 
Tj={zeC: gb jE -wo < 2) < +00}. 


v2 V2 


There is a useful class of functions for which the integrals in the definition 
of f and f exist. 
For u € C™(IR”) we set 


aja) = mare, (sup ((1 + [eP)!/2|D%u(e)))). 


a: jal<k \ Rn 


Definition. The Schwartz class of rapidly degreasing functions is 


S(R”") = {ue C™(R") : a3.4(u) < co for any j,k Ee NU {O}}. 


This space is a Frechét space. 
Proposition 5.3. Assume u € S(R"), then U and u € S(R"). 


Proof. See [24], Chapter 1.2, for example, or an exercise. 


5.1.1 Pseudodifferential operators 


The properties of Fourier transform lead to a general theory for linear partial 
differential or integral equations. In this subsection we define 


a 


Dy, = = til PRON 
i Oxy’ ’ gies 


and for each multi-index a@ as in Subsection 3.5.1 


D® = D®,,, D&, 
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Thus 
1 ala 


D* = ; 
ilol Oat! ...Oxn” 


Let 


p(t,D):= S> aa(x)D*, 


ja|<m 


be a linear partial differential of order m, where dq are given sufficiently 
regular functions. 

According to Theorem 5.1 and Proposition 5.3, we have, at least for 
wu € S(R”), 


u(t) = (2m)? i el” (e) dl, 


n 


which implies 


n 


D° u(x) = Cas ael el SE°R(€) dé. 
Consequently 
ple, D)u(x) = (2n)-"?? f olp(a, ya(e) a, (5.9) 


where 


p(a,€) = So aala)é*. 


ja|<m 


The right hand side of (5.9) makes sense also for more general functions 
p(x,€), not only for polynomials. 


Definition. The function p(x, &) is called symbol and 


(Pu)(w) = (2n)-"/? | el Spl, €)A(€) dé 


n 


is said to be pseudodifferential operator. 


An important class of symbols for which the right hand side in this defini- 
tion of a pseudodifferential operator is defined is S™ which is the subset of 
p(x, €) € C*(Q x R”) such that 


|D2DEep(a,£)| < Cka,a(p) (1 + [él 


for each compact K CQ. 
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Above we have seen that linear differential operators define a class of pseu- 
dodifferential operators. Even integral operators can be written (formally) 
as pseudodifferential operators. Let 


(Pu)(x) = be K(x, y)u(y) dy 


be an integral operator. Then 
(Pyle) = Qn? f Kaw) fe terate) ag 
= enn f oet(f olesx(e,y) av) a) 
Then the symbol associated to the above integral operator is 


p(a,é) =f 8K (4) dy 
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5.2 Exercises 
1. Show 
| en WP? dy = (an)"/2, 
2. Show that u € S(R”) implies @, wu € S(R”). 
3. Give examples for functions p(,€) which satisfy p(x,€) € S™. 


4. Find a formal solution of Cauchy’s initial value problem for the wave 
equation by using Fourier’s transform. 


Chapter 6 


Parabolic equations 


Here we consider linear parabolic equations of second order. An example is 
the heat equation 


up = a’Au, 


where u = u(a,t), « € R°, t > 0, and a? is a positive constant called 
conductivity coefficient. The heat equation has its origin in physics where 
u(x,t) is the temperature at x at time t, see [20], p. 394, for example. 


Remark 1. After scaling of axis we can assume a = 1. 


Remark 2. By setting t := —t, the heat equation changes to an equation 
which is called backward equation. This is the reason for the fact that 
the heat equation describes irreversible processes in contrast to the wave 
equation Ou = O which is invariant with respect the mapping t > —t. 
Mathematically, it means that it is not possible, in general, to find the 
distribution of temperature at an earlier time ¢t < to if the distribution is 
given at to. 


Consider the initial value problem for u = u(x,t), u € C™(R" x Ri), 


u = Au inzxeR", t>0, (6.1) 


where ¢ € C(R”) is given and A= A,. 
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6.1 Poisson’s formula 


Assume wu is a solution of (6.1), then, since Fourier transform is a linear 
mapping, 
Ut — Au = 0. 


From properties of the Fourier transform, see Proposition 5.1, we have 


provided the transforms exist. Thus we arrive at the ordinary differential 
equation for the Fourier transform of u 


du 


Q2~ 
pid a) 
qt ere =o, 


where € is considered as a parameter. The solution is 
f(g, t) = (je 


since u(€,0) = (€). From Theorem 5.1 it follows 


Wat) SQ)? bee KteE* de 
R” 


en" f o)( fete ae) ay 


K(a,y,t) = (20) i eif(a—u)1E07* ag. 


n 


Set 


By the same calculations as in the proof of Theorem 5.1, step (vi), we find 


EO CG ut) oe (Ant) —7/2@7 lel? /4t, (6.3) 
Thus we have 
1 2 
t)= = See ar oF 6.4 
ee) = Grape f, z (6.1) 


Definition. Formula (6.4) is called Poisson’s formula and the function K 
defined by (6.3) heat kernel or fundamental solution of the heat equation. 
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Figure 6.1: Kernel K(2,y,t), p = |x — yl, ti < ta 


Proposition 6.1 The kernel K has following properties: 
(i) K(a,y,t) € C@(R” x R" x R4), 

(ii) (0/0t — A)K(a,y,t) =0, t>0, 

(itt) K(x,y,t) > 0, t > 0, 

(0) fon K@y3t) dy =1, ve Rt > 0, 


(v) For each fited 6 > 0: 


lim K(a2,y,t) dy =0 
t— 0 YR"\Bs(2) 
t>0 
uniformly for x € R”. 


Proof. (i) and (iii) are obviously, and (ii) follows from the definition of K. 
Equations (iv) and (v) hold since 


| Kay a).dy. = i: (4rt)—7/2—l-¥P/4 dy 
R”\ Bs (x) R"\ Bs (x) 


gon? ‘| eM an 
B"\ Bs, /ap(0) 
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by using the substitution y = a + (4t)!/2n. For fixed 6 > 0 it follows (v) and 
for 6 := 0 we obtain (iv). 


Theorem 6.1. Assume ¢ € C(R") and suppn |¢(x)| < oo. Then u(z,t) 
given by Poisson’s formula (6.4) is in C°(R” x Rz), continuous on R” x 
[0,co) and a solution of the initial value problem (6.1), (6.2). 


Proof. It remains to show 


lim u(z,t) = (2). 
r—€ 
t— 0 


Since ¢ is continuous there exists for given ¢ > 0 a6 = 6(€) such that |¢(y)— 


g E+5 +28 


Figure 6.2: Figure to the proof of Theorem 6.1 


o(€)| < € if |y —€| < 26. Set M := supgn |6(y)|. Then, see Proposition 6.1, 
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It follows, if |x — €| < 6 and t > 0, that 
lu(a,t) - 4(6)| < | oy Keewr8) 16) ~ 96) dy 
5 (x 


5 i K(x,y,t) |o(y) — 6(€)| dy 
R"\B;(«) 


< | K(x,y,t) |o(y) — o(6)| dy 
Bo5(x) 
42M K(a,y,t) dy 
R”\ Bs (a) 
< ef K(z,y,t)dy+2M K(x,y,t) dy 
Rn R"\ Bs (x) 
< Qe 


if 0 <t < to, to sufficiently small. 
Remarks. 1. Uniqueness follows under the additional growth assumption 
|u(a,t)| < Mel” in Dr, 


where M and a are positive constants, see Proposition 6.2 below. 
In the one-dimensional case, one has uniqueness in the class u(x,t) > 0 in 
Dr, see [10], pp. 222. 


2. u(x,t) defined by Poisson’s formula depends on all values ¢(y), y € R”. 
That means, a perturbation of ¢, even far from a fixed x, has influence to the 
value u(x,t). In physical terms, this means that heat travels with infinite 
speed, in contrast to the experience. 


6.2 Inhomogeneous heat equation 
Here we consider the initial value problem for u = u(x,t), u€ C°(R" x R4), 


u—-Au = f(x,t) nz eR", t>0, 
u(a, 0) (2), 


where ¢ and f are given. From 


nee 


u, — Au = f(z,t) 
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we obtain an initial value problem for an ordinary differential equation: 
du = 

dt 
u(g,0) = 


Nn 


(¢,t) 
(¢). 


Pa = f 
e 
The solution is given by 


iE, t) = e PG (eg) + | oe PC=1) FE, 7) dr. 


0 


Applying the inverse Fourier transform and a calculation as in the proof of 
Theorem 5.1, step (vi), we get 


(at) = (any? felt (erlePG(@ 
t 
“I8P—1) Fe 7) dr) dé. 
+f WPF E,7) ar) ag 


From the above calculation for the homogeneous problem and calculation as 
in the proof of Theorem 5.1, step (vi), we obtain the formula 


1 2 
it) = —ly-al?/(4t) g 
u(x, t) avai Inn oye y 
t 
+f er ely 2l?/(A(t-7)) dy dr. 


(2Va(t=7)) 


This function u(x,t) is a solution of the above inhomogeneous initial value 
problem provided 
@€C(R"), sup|9(x)| < co 
R” 


and if 


f € C(R” x [0,00)), M(r) := sup|f(y,7T)| < 0, O< T < oo. 
R. 


6.3. Maximum principle 
Let 2 c R” be a bounded domain. Set 


Dr = 2x(0,T), T>0, 
Sr {(x,t): (x,t) € Q x {0} or (a, t) € OQ x [0,T]}, 
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Figure 6.3: Notations to the maximum principle 


see Figure 6.3 


Theorem 6.2. Assume u € C(Dr), that Ut, Ux,2, exist and are continuous 
in Dr, and 
up—-Au<0 in Dr. 


Then 


max u(x,t) = maxu. 
Dr or 


Proof. Assume initially u— Au < 0 in Dr. Let ¢ > 0 be small and 
O0<e<T. Since u € C(Dr_z), there is an (x9, to) € Dr_- such that 
u(xo, to) = max u(z,t). 


Dr_< 


Case (i). Let (%o,to) € Dr—-. Hence, since Dr_- is open, uz(xo, to) = 0, 
Uz, (x9, to) =0,1=1,...,n and 


S- Uzn, (£0, to)GiGe <0 for all ¢ € R”. 
Lk=1 


The previous inequality implies that uz,2,(xo,to) <0 for each k. Thus we 
arrived at a contradiction to u,— Au <0 in Dr. 
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Case (ii). Assume (20, to) € 2x {T—e}. Then it follows as above Au < 0 in 
(xo, to), and from u(2, to) > u(xo,t), t < to, one concludes that u,(xo, to) > 
0. We arrived at a contradiction to u,— Au < 0 in Dp again. 


Summarizing, we have shown that 


max u(x,t) = max u(z,t). 
Dr_« T-e 


Thus there is an (x-,t-) € Sp_- such that 


u(@e,te) = max u(x,t). 
Dr_- 


Since u is continuous on Dy, we have 


lim max u(x,t) = max u(z,t). 
e—0 Dr-< Dr 


It follows that there is (%,t) € Sr such that 


u(#,t) = max u(z,t) 
Dr 


since Sp_- C Sy and Sr is compact. Thus, theorem is shown under the 
assumption u,— Au < 0 in Dr. Now assume uy — Au < 0 in Dr. Set 


v(a, t) := u(x,t) — kt, 
where k is a positive constant. Then 
vu, —- Av =u,— Au—k <0. 


From above we have 


maxu(x,t) = max(v(x,t) + kt) 
Dr Dr 
< maxu(2,t)+kT 
Dr 


= maxvo(z,t)+kT 
Sr 


IA 


max u(x,t) + kT, 
Sr 


Letting k — 0, we obtain 


max u(x,t) < max u(z,t). 
Dr Sr 
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Since Sy; C Dr, the theorem is shown. 


If we replace in the above theorem the bounded domain 2 by R”, then the 
result remains true provided we assume an additional growth assumption 
for u. More precisely, we have the following result which is a corollary of 
the previous theorem. Set for a fixed T,0<T < ow, 


Dr S160) ee RY Ot ST}, 
Proposition 6.2. Assume u € C(Dr), that uz, Us;2, exist and are contin- 


uous in Dr, 
uzp—-Au<0 in Dr, 


and additionally that u satisfies the growth condition 
u(x,t) < Mell, 
where M and a are positive constants. Then 
max u(x,t) = maxu. 
Dr Sr 
It follows immediately the 


Corollary. The initial value problem uw — Au = 0 in Dr, u(x,0) = f(x), 
x € R”, has a unique solution in the class defined by u € C(Dr), ut, Ux;x, 


exist and are continuous in Dp and |u(z,t)| < Mele, 

Proof of Proposition 6.2. See [10], pp. 217. We can assume that 4aT < 1, 
since the finite interval can be divided into finite many intervals of equal 
length 7 with 4ar < 1. Then we conclude successively for k& that 


u(x,t) < sup u(y, kr) < sup u(y, 0) 
yeER” yeER” 
for kr <t< (k+1)r, k=0,...,N—1, where N = T/r. 
There is an € > 0 such that 4a(T'+ €) < 1. Consider the comparison 
function 


vy(a,t): = u(x,t) — w(4n(T +6 —t)) 7/7 ele oP /AT +B) 
= u(a,t) —wK (ia, iy,T +e -t) 
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for fixed y € R” and for a constant js > 0. Since the heat kernel K (iz, iy, t) 
satisfies K; = AK,, we obtain 


O 
Apt — Ste = te — Aus 0. 


Set for a constant p > 0 
Dry = {(z,t): |jz—y| <p, 0<t < TH. 
Then we obtain from Theorem 6.2 that 
Uly, t) < max Uns 
ST,p 
where Sp, = Sp of Theorem 6.2 with Q = B,(y), see Figure 6.3. On the 
bottom of S',, we have, since wk > 0, 


Up(@,0) < u(x,0) < op f(z). 


On the cylinder part |x — y| = p,0 <t<T, of Sr, it is 


vu(z,t) < Me? — y(4n(T +6 —t))0? eh /AC te) 
< Merlllte!” — y(4a(T +6)? (A+) 
< sup f(z) 
zER” 


for all p > po(), po sufficiently large. We recall that 4a(T +) < 1. 
Summarizing, we have 


max v,(x,t) < sup f(z) 
ST,p zER” 


if p > po(). Thus 


vp(y,t) < max v,p(x,t) < sup f(z) 
ST, p zER” 


if p > po(s). Since 


ouly,t) = uly,t) —w(4e(T +e-t))"”? 
it follows 


u(y,t) — w(4n(T +e—t))~"/”? < sup f(z). 
zER” 
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Letting « — 0, we obtain the assertion of the proposition. 


The above maximum principle of Theorem 6.2 holds for a large class of 
parabolic differential operators, even for degenerate equations. Set 
n 
Lu = So ali(a,A) 
ij=l 
where a) € C(Dr) are real, a’? = a/*, and the matrix (a’’) is nonnegative, 
that is, 


> a (x, t)G¢; >0 for all ¢ € R", 
ij=l 
and (2,t) € Dr. 


Theorem 6.3. Assume u € C(Dr), that uz, Ur,2, etist and are continuous 
in Dr, and 

—Iu<0 in Dp. 
Then 


max u(x,t) = max wu. 
Dr oT 


Proof. (i) One proof is a consequence of the following lemma: Let A, B real, 
symmetric and nonnegative matrices. Nonnegative means that all eigenval- 
ues are nonnegative. Then trace (AB) = + jel qs bi; > 0, see an exercise. 


(ii) Another proof exploits transform to principle axis directly: Let U = 
(21,---;2n), where z is an orthonormal system of eigenvectors to the eigen- 
values \; of the matrix A = (a’/ (zg, to)). Set € = Un, « = UT (x — 29)y and 
v(y) = u(x + Uy, to), then 


n 


0O< So a’ T(a9, to)GiGj = Yi 


ij=l 
n 
= 2 
0 = y ies Vyrys hi 
ij=l i=l 


It follows A; > 0 and vy,y, <0 for all 7. Consequently 


n 


Soa’ I(x9, to)Uz, a; (xo, to) -y NVyy;, < 0. 


j=l 
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6.4 Initial-boundary value problem 


Consider the initial-boundary value problem for c = c(z,t) 


c¢ = DA&Ac inQ x (0,00) (6.5) 
c(z,0) = oo(x) rE 
a = 0 on 02 x (0,00). (6.7) 


Here is 2 C R”, n the exterior unit normal at the smooth parts of 00, Da 
positive constant and co(z) a given function. 


Remark. In application to diffusion problems, c(x,t) is the concentration 
of a substance in a solution, co(x) its initial concentration and D the coef- 
ficient of diffusion. 

First Fick’s rule says that w = DOc/0n, where w is the flow of the substance 
through the boundary 09. Thus according to the Neumann boundary con- 
dition (6.7), we assume that there is no flow through the boundary. 


6.4.1 Fourier’s method 


Separation of variables ansatz c(z,t) = v(a#)w(t) leads to the eigenvalue 
problem, see the arguments of Section 4.5, 


—Av = Av inQD (6.8) 
Ov 
a 0 on OQ, (6.9) 


and to the ordinary differential equation 
w'(t) + A\Dw(t) = 0. (6.10) 


Assume Q is bounded and OQ) sufficiently regular, then the eigenvalues 
of (6.8), (6.9) are countable and 


0= Ao < Ay S AQ <<... 9 Ow. 


Let v;(z) be a complete system of orthonormal (in L?(Q)) eigenfunctions. 
Solutions of (6.10) are 
w;(t) = Cre", 


where C; are arbitrary constants. 
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According to the superposition principle, 
N 
Cnn = LS Cye Pity; (a) 
j=0 
is a solution of the differential equation (6.8) and 
co 
c(z, t) = ye Cye Pi", (2), 
j=0 


with 


C; =| co(x)v;(x) da, 


is a formal solution of the initial-boundary value problem (6.5)-(6.7). 


Diffusion in a tube 


Consider a solution in a tube, see Figure 6.4. Assume the initial concentra- 


Figure 6.4: Diffusion in a tube 


tion co(21, 22,23) of the substrate in a solution is constant if x3 = const. 
It follows from a uniqueness result below that the solution of the initial- 
boundary value problem c(21, 72,73, t) is independent of x; and 22. 
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Set z = x3, then the above initial-boundary value problem reduces to 


Ch = Des 


c(z,0) = coz) 
cz, = 0, z=0, z=l. 


The (formal) solution is 


ioe) ahd 
c(z,t) = DF mc, * cos (Snz) ; 
n=0 
where 
1 l 
Co = al co(z) dz 
l 0 
l 
Cr, = al Co(z) cos (Fnz) dz, n>1 
l 0 


6.4.2 Uniqueness 


Sufficiently regular solutions of the initial-boundary value problem (6.5)- 
(6.7) are uniquely determined since from 


ce = DA&Ac inQ x (0,00) 
c(z,0) = 0 


Oc 
Aan 0 on 02 x (0,00). 


it follows that for each 7 > 0 


0 = i. ‘a (cee — D(Ac)c) dxdt 

= ail J 2 

= i | sail’ 2) dtd +D f i |Vac|” dadt 
>/ e(a,T) ax +d | [ |Vac|? dadt. 
2Jo a Jo 


6.5 Black-Scholes equation 


Solutions of the Black-Scholes equation define the value of a derivative, for 
example of a call or put option, which is based on an asset. An asset 
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can be a stock or a derivative again, for example. In principle, there are 
infinitely many such products, for example n-th derivatives. The Black- 
Scholes equation for the value V(.S,t) of a derivative is 


1 
Yt 57 5 Ves +rSVs—rV =0 ing, (6.11) 
where for a fixed T,0<T <«o, 
QS{GHeR? VS F< oo, 0StS Th, 


and o, r are positive constants. More precisely, 
a is the volatility of the underlying asset S, 
r is the guaranteed interest rate of a risk-free investment. 

If S(t) is the value of an asset at time t, then V(S(t),t) is the value 
of the derivative at time ¢, where V(S,t) is the solution of an appropriate 
initial-boundary value problem for the Black-Scholes equation, see below. 

The Black-Scholes equation follows from Ito’s Lemma under some as- 
sumptions on the random function associated to S(t), see [26], for example. 


Call option 


Here is V(S,t) := C(.S,t), where C(S,t) is the value of the (European) call 
option. In this case we have following side conditions to (6.11): 


C(S,T) = max{S — E,0} (6.12) 
C(0,t) = 0 (6.13) 
C(S,t) = S+o0(S) as S > o, uniformly in t, (6.14) 


where F and T are positive constants, E’ is the exercise price and T the 
expiry. 


Side condition (6.12) means that the value of the option has no value at 
time T if S(T) < E, 
condition (6.13) says that it makes no sense to buy assets if the value of 
the asset is zero, 
condition (6.14) means that we buy assets if its value becomes large, see 
Figure 6.5, where the side conditions are indicated. 


Theorem 6.4 (Black-Scholes formula for European call options). The so- 
lution C(S,t),0<S <o,0<t<T, of the initial-boundary value prob- 
lem (6.11)-(6.14) is explicitly known and is given by 


C(S,t) = SN(di) — EeT? 9 N(d), 
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C=max{S-E,0} 


S 
Figure 6.5: Side conditions for a call option 
where 
VES = / ane 
V2n —0oo 
4 = In($/E) + (r + o?/2)(T — t) 
, oVT—-t ; 
b = In(S/E) + (r — 07 /2)(T —t) 
: oVT—-t ; 
Proof. Substitutions 
x T 
S=Ee*, t=T- a2)’ C = Ev(a,7T) 
change equation (6.11) to 
Ur = Ure + (Kk — 1)vz — kv, (6.15) 
where a 
k= ; 
o?/2 
Initial condition (6.15) implies 
u(x, 0) = max{e* — 1,0}. (6.16) 


For a solution of (6.15) we make the ansatz 


y =e +Pru(a, 7), 
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where a and @ are constants which we will determine as follows. Inserting 
the ansatz into differential equation (6.15), we get 


But ur, = a7 ut 2aug + Ugg + (k 1)(au + uz) — ku. 


Set G = a? + (k —1)a—k and choose a such that 0 = 2a + (k — 1), then 
Ur = Ure. Thus 


v(a,7) =e FD) 4/4 (7), (6.17) 


where u(z,7) is a solution of the initial value problem 


Up = Ugr, —-EO<X<w, T>OD 
u(x,0) = uo(x), 


with 
uo (x) = max {elt+be/2 _ e(k-N2/2 0} 


A solution of this initial value problem is given by Poisson’s formula 
Dre ssf ug(s)e~@-9)"/(47) de, 
, 2/7T Joo 
Changing variable by g = (s — )/(V2T), we get 


1 Bice 
u(z,T) = =| uo(qV2r + x)e~7/? dq 


= ,-h, 
where 

1 ee 2 

Pt. ee i: elhtIV(ctaV27) 9/2 dg 
V2n J—a/(\/2r) 
1 ims 2 

fy nas ee / elk-W(etaV3) 9/2 gg, 
V2n J—a/(\/2r) 


An elementary calculation shows that 


= e(k+1)2/2+(k+1)?7/4 7 (d,) 
in elk 1a /2+(k-1)?7/4 (gy), 
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where 
d eg ied oe 
: V2r 2 
x 1 — 
dg = k—1)V2 
eh i 


N(dj) e*/2 ds, G=1, 2. 


Tale 


Combining the formula for u(x, 7), definition (6.17) of v(#,7) and the previ- 
ous settings 2 = In(S/E), rT = 0?(T — t)/2 and C = Ev(z,7), we get finally 
the formula of of Theorem 6.4. 


In general, the solution u of the initial value problem for the heat equa- 
tion is not uniquely defined, see for example [10], pp. 206. 


Uniqueness. The uniqueness follows from the growth assumption (6.14). 
Assume there are two solutions of (6.11), (6.12)-(6.14), then the difference 
W(S,t) satisfies the differential equation (6.11) and the side conditions 


W(S,T) =0, W(0,t) =0, W(S,t) = O(S) as S — co 
uniformly inO<t<T. 
From a maximum principle consideration, see an exercise, it follows that 


|W(S,t)| < cS on S>0,0<t<T. The constant c is independent on S 
and t. From the definition of u we see that 


1 
u(x,T) = =e" PWS, t), 
E 
where S$ = Ee”, t = T — 27/(07). Thus we have the growth property 
|u(a,7)| < Me", ceR, (6.18) 


with positive constants M and a. Then the solution of u; = ugy, in —co < 
z<oo,0<7 < 0°T/2, with the initial condition u(a,0) = 0 is uniquely 
defined in the class of functions satisfying the growth condition (6.18), see 
Proposition 6.2 of this chapter. That is, u(x, 7) = 0. 
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Put option 


Here is V(S,t) := P(S,t), where P(S,t) is the value of the (European) put 
option. In this case we have following side conditions to (6.11): 


P(S,T) = max{E— S,0} (6.19) 
P(0,t) = Ee(?-4 (6.20) 
P(S,t) = o(S) as Soo, uniformly in0<t < T. (6.21) 


Here E is the exercise price and T the expiry. 


Side condition (6.19) means that the value of the option has no value at 
time T if S(T) > E, 
condition (6.20) says that it makes no sense to sell assets if the value of the 
asset is zero, 
condition (6.21) means that it makes no sense to sell assets if its value 
becomes large. 


Theorem 6.5 (Black-Scholes formula for European put options). The solu- 
tion P(S,t),0<S<c,t< T of the initial-boundary value problem (6.11), 
(6.19)-(6.21) is explicitly known and is given by 


P(S,t) = Ee-™?—-9 N(—d2) — SN(—d1) 
where N(x), di, dg are the same as in Theorem 6.4. 


Proof. The formula for the put option follows by the same calculations as 
in the case of a call option or from the put-call parity 


C(S,t) — P(S,t) = S — Ee") 


and from 
N(x) + N(—2) = 1. 


Concerning the put-call parity see an exercise. See also [26], pp. 40, for a 
heuristic argument which leads to the formula for the put-call parity. 
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Exercises 
Show that the solution u(x,t) given by Poisson’s formula satisfies 
inf p(z) < u(a,t) < sup y(z) , 
zER" zeER” 


provided y(x) is continuous and bounded on R”. 


. Solve for given f(x) and yw € R the initial value problem 


UptUe + MUerer = O in RXR, 
u(z,0) = f(x). 


. Show by using Poisson’s formula: 


(i) Each function f € C([a,6]) can be approximated uniformly by a 
sequence f, € C™[a, }] . 

(ii) In (i) we can choose polynomials f,, (Weierstrass’s approximation 
theorem). 


ly-a/? 


Hint: Concerning (ii), replace the kernel K = exp(—“3j-—) by a se- 


_ ly-2/? 
qt 


quence of Taylor polynomials in the variable z = 


. Let u(x,t) be a positive solution of 


Ut = Urge, € > 0, 


where yw is a constant. Show that @ := —2yu,/u is a solution of 
Burger’s equation 
0: + O00, = HOr2, t > 0. 


. Assume uj(5,¢), ..., Un(s,t) are solutions of uz = uss. Show that []j;_, ux(re, ¢) 


is a solution of the heat equation u;— Au = 0 in R” x (0,00). 


. Let A, B are real, symmetric and nonnegative matrices. Nonnegative 


means that all eigenvalues are nonnegative. Prove that trace (AB) = 
aa at bj; 2 0. 

Hint: (i) Let U = (21,...,2n), where z is an orthonormal system of 
eigenvectors to the eigenvalues A; of the matrix B. Then 


Jr (0 ee 
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10. 


11. 


is a square root of B. We recall that 


dM 0 0 
epee ek : 
0 0 Ne 


(ii) trace (QR) =trace (RQ). 

(iii) Let jui(C), ... n(C) are the eigenvalues of a real symmetric n x n- 
matrix. Then trace C = )7/_, «(C), which follows from the funda- 
mental lemma of algebra: 


det (AI — C) ee (ohne eee oer Pome ee 
(A — pi) +... + (A= bn) 


DS iis ont Pee 


. Assume 2) is bounded, wu is a solution of the heat equation and wu sat- 


isfies the regularity assumptions of the maximum principle (Theorem 
6.2). Show that u achieves its maximum and its minimum on Sr. 


. Prove the following comparison principle: Assume 2 is bounded and 


u, v satisfy the regularity assumptions of the maximum principle. Then 


u—Au < uy—Av inDr 


u < vonsp 


imply that u < v in Dr. 


. Show that the comparison principle implies the maximum principle. 


Consider the boundary-initial value problem 


u—- Au = f(ax,t) in Dr 
u(a, t) (x,t) on Sp, 


where f, ¢ are given. 
Prove uniqueness in the class u, w, Uz,2, € C(Dr). 


Assume u, v1, v2 € C?(Dr) A C(Dr), and u is a solution of the 
previous boundary-initial value problem and v1, v2 satisfy 
(vi),- Av, < f(ax,t) < (v2), -— Avg in Dr 
VU < @< v2 on Sp. 
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13. 


14. 


15. 
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Show that (inclusion theorem) 


v1 (x,t) < u(a,t) < vo(z,t) on Dr. 


Show by using the comparison principle: let u be a sufficiently regular 
solution of 


u—-Au = 1 in Dr 
0 on Sr, 


U 


then 0<u(z,t)<t in Dp. 


Discuss the result of Theorem 6.3 for the case 
n n 
Lu= ye Gig (2, t)Ueje;, + Se b;(z,t)uz, + c(x, t)u(z, t). 
ij=l i 
Show that - 
ot = Se cnet sin(nx), 
n=1 
where 


2 TT 
C= - | f(x) sin(nz) da, 
T JO 
is a solution of the initial-boundary value problem 


Ut = Ugr, © € (0,7), t>0, 


f (2), 


S 
— 
=) 

eH et © 
So Ne 

I 

oO 


if f € C*(R) is odd with respect to 0 and 27-periodic. 


(i) Find the solution of the diffusion problem c; = Dc,, in0< z <I, 
0<t< oo, D=const. > 0, under the boundary conditions c,(z,t) = 
0 if z =0 and z =/ and with the given initial concentration 


co= const. if O<z<h 
c(2,0) = eo(2) =f : 0 if h<z<l. 


(ii) Calculate limy... c¢(z, t). 
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16. Prove the Black-Scholes Formel for an European put option. 


Hint: Put-call parity. 
17. Prove the put-call parity for European options 
C(S,t) — P(S,t) = S — Ee"T-9) 


by using the following uniqueness result: Assume W is a solution 
of (6.11) under the side conditions W(S,T) = 0, W(0,t) = 0 and 
W(S,t) = O(S) as S > o, uniformly on 0 <t¢< T. Then W(S,t) = 
0. 


18. Prove that a solution V (5, t) of the initial-boundary value problem (6.11) 
in Q under the side conditions (i) V(S,T) = 0, S > 0, (ii) V(0,¢) = 0, 
0<t<T, (iii) img. V(S,t) = 0 uniformly in 0 < t < T, is uniquely 
determined in the class C?(Q) 9 C(Q). 


19. Prove that a solution V (5, t) of the initial-boundary value problem (6.11) 
in 9, under the side conditions (i) V(S,T) = 0, S > 0, (ii) V(0,¢) = 0, 
0<t<T, (iii) V(S,t) = S+o0(S) as S > ov, uniformly on0<t<T, 
satisfies |V(S,t)| < cS for all S>Oand0<t<T. 


Chapter 7 


Elliptic equations of second 
order 


Here we consider linear elliptic equations of second order, mainly the Laplace 
equation 
Au = 0. 


Solutions of the Laplace equation are called potential functions or harmonic 
functions. The Laplace equation is called also potential equation. 
The general elliptic equation for a scalar function u(x), 2 € 2 C R", is 


n 


Lu:= S- a" (x) Ueray + — b) (x) ue, + c(x)u = f(x), 
j=l 


t,j=1 


where the matrix A = (a’’) is real, symmetric and positive definite. If A is 
a constant matrix, then a transform to principal axis and stretching of axis 


leads to 
n 


a = 
y a Ugix; = Av, 


t,j=1 


where u(y) := u(Ty), T stands for the above composition of mappings. 


7.1 Fundamental solution 


Here we consider particular solutions of the Laplace equation in R” of the 
type 
u(x) = f(jz— yl), 
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where y € R” is fixed and f is a function which we will determine such that 
u defines a solution if the Laplace equation. 
Set r= |x — y|, then 


1, \ bi Yi 
fry 


tee, = fn) = vi | ra (2 ee 


r 


S 
8 
I| 


Thus a solution of Au = 0 is given by 


_ qinr+c : n=2 
ora ee ee : n>3 


with constants cy, C2. 

Definition. Set r = |x — y|. The function 

Ge n>3 

is called singularity function associated to the Laplace equation. Here 


is Wy the area of the n-dimensional unit sphere which is given by w, = 
2n”/? /D(n/2), where 


lee) 
I(t) =| e ’p'' dp, t>0, 
0 
is the Gamma function. 


Definition. A function 


y(x,y) = s(r) + O(a, y) 


is called fundamental solution associated to the Laplace equation if ¢@ € 
C?(Q) and A,¢ = 0 for each fixed y € Q. 


Remark. The fundamental solution y satisfies for each fixed y € 2 the 
relation 


ap (x, y)A, P(x) dx = ®(y) for all BE C2(), 


7.2. REPRESENTATION FORMULA 177 


see an exercise. This formula follows from considerations similar to the next 
section. 
In the language of distribution, this relation can be written by definition 
as 
—Avy(,y) = O(a — y), 


where 6 is the Dirac distribution, which is called 6-function. 


7.2 Representation formula 


In the following we assume that 0, the function ¢ which appears in the 
definition of the fundamental solution and the potential function u consid- 
ered are sufficiently regular such that the following calculations make sense, 
see [6] for generalizations. This is the case if 2 is bounded, OQ is in C’, 
@ € C?(Q) for each fixed y € Q and u € C?(Q). 


ny 


Figure 7.1: Notations to Green’s identity 


Theorem 7.1. Let wu be a potential function and y a fundamental solution, 
then for each fired y € Q 


w= ff (nen G2 - ue BE?) as, 


Proof. Let B,(y) C Q be a ball. Set O,(y) = Q\ Bp(y). See Figure 7.2 for 
notations. From Green’s formula, for u, v € C?(Q), 


| (vAu — uAv) dx = | (os = uw) dS 
p(y) an,(y) \ On On 
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Figure 7.2: Notations to Theorem 7.1 


we obtain, if v is a fundamental solution and u a potential function, 


| (5 ux) dS = 0. 
IQp(y) On On 
Thus we have to consider 
| soe ds: = | nee ds +f poe dS 
IQp(y) on an On OB,(y) On 


if ee dS = | get ds +f get dS. 
OQp(y) dn 0Q an OBp(y) an 


We estimate the integrals over 0B,(y): 


(i) 


IA 


M |v| dS 
OBp(y) 


M | s(p) dS + Cunp" J, 
OBo(y) 


IA 


where 


Ne 2 WAG Sane aula pS 
Boo (y) 


C = Cly)= sup |¢(z,y)|. 
LEBpo (y) 
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From the definition of s(p) we get the estimate as p — 0 
Ou { O(p|Inp|) : n=2 
v— dS = 7.1 
oe On O(p) : n23. co 


(ii) Consider the case n > 3, then 


1 1 
i Rae ees uy as+ f 2? as 
ABy(y) On Wn JaByy) P ABy(y) On 


= af vas+0u") 


wy p” 1 


for an xo € OB,(y). 


Combining this estimate and (7.1), we obtain the representation formula of 
the theorem. 


Corollary. Set ¢ = 0 and r = |x — y| in the representation formula of 
Theorem 7.1, then 


—U 


nr— 
Ong Ony 


e 
ite 
= 
= 

I 


= (1 ee Tet) de, 0 =2, (7.2) 
Qn fale) 
1 1 Ou O(r?-”) 


uy) = ree (a5 a ) dSz, n> 3. (7.3) 


7.2.1 Conclusions from the representation formula 


Similar to the theory of functions of one complex variable, we obtain here 
results for harmonic functions from the representation formula, in particular 
from (7.2), (7.3). We recall that a function u is called harmonic if u € C?(Q) 
and Au=0 in 2. 


Proposition 7.1. Assume u is harmonic inQ. Then u€ C™(Q). 


Proof. Let Qo CC Q be a domain such that y € Qo. It follows from 
representation formulas (7.2), (7.3), where Q := Qo, that D!u(y) exist and 
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are continuous for all / since one can change differentiation with integration 
in right hand sides of the representation formulae. 


Remark. In fact, a function which is harmonic in 2 is even real analytic 
in Q, see an exercise. 


Proposition 7.2 (Mean value formula for harmonic functions). Assume u 
is harmonic in Q. Then for each By(x) CCQ 


reo) 
u(x) = ——— U dS,. 
©) = Tt fo oy CY 


Proof. Consider the case n > 3. The assertion follows from (7.3) where 
Q := B,(«) since r = p and 


1 1 
| = Ou dSy = —> | Os 
OBy(x) 1” ONy pr Jap, (a) My 


1 
= =a | Au dy 
p Bp(x) 
0. 


We recall that a domain 2 € R” is called connected if Q is not the union of 
two nonempty open subsets 21, Q2 such that Q,NQ2 = 9. A domain in R” 
is connected if and only if its path connected. 


Proposition 7.3 (Maximum principle). Assume u is harmonic in a con- 
nected domain and achieves its supremum or infimum inQ. Then u = const. 
inQ. 


Proof. Consider the case of the supremum. Let xo € 2 such that 


u(xo) = sup u(x) =: M. 
Q 


Set OQ, := {2 € OD: ule) = M} and Oo := {a € 1: ule) < M}. The 
set 1 is not empty since xp € 2). The set Q2 is open since wu € C20): 
Consequently, Q2 is empty if we can show that 2; is open. Let % € 01, then 
there is a py > 0 such that B,,(Z) C Q and u(x) = M for all x € B,, (2). 
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If not, then there exists p > 0 and & such that |x — Z| = p, 0 < p < po and 
u(%) < M. From the mean value formula, see Proposition 7.2, it follows 


1 M 
Wn p dBp (2) Wn Pp OBp(Z) 


which is a contradiction. Thus, the set Q2 is empty since Q, is open. 


Corollary. Assume 2 is connected and bounded, and u € C?(Q) NM C(Q) 
is harmonic in 2. Then wu achieves its minimum and its maximum on the 
boundary 02. 


Remark. The previous corollary fails if Q is not bounded as simple coun- 
terexamples show. 


7.3. Boundary value problems 
Assume 2 C R” is a connected domain. 


7.3.1 Dirichlet problem 


The Dirichlet problem (first boundary value problem) is to find a solution 
u € C2(Q) N C(Q) of 


Au = 0 inQ (7.4) 
u = ® on OQ, 


where ® is given and continuous on 02. 


Proposition 7.4. Assume Q is bounded, then a solution to the Dirichlet 
problem is uniquely determined. 


Proof. Maximum principle. 
Remark. The previous result fails if we take away in the boundary condi- 
tion (7.5) one point from the the boundary as the following example shows. 


Let 2 C R? be the domain 


Q = {x € By (0): x2 > 0}, 
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% 


Figure 7.3: Counterexample 


Assume u € C?(Q)9 C(Q \ {0}) is a solution of 


Au = 0 inQd 
0 on O02 \ {0}. 


U 


This problem has solutions u = 0 and u = Im(z +271), where z = 21 + ize. 
Another example see an exercise. 

In contrast to this behaviour of the Laplace equation, one has uniqueness 
if Au = 0 is replaced by the minimal surface equation 


0 Un, to) Une 
=0 
Oxy (atts | Ox2 (San) 
7.3.2 Neumann problem 


The Neumann problem (second boundary value problem) is to find a solution 
u € C2(Q) NCQ) of 


Pa = Dio (7.6) 
Ou 
as ® on OQ, (7.7) 


where ® is given and continuous on OQ. 


Proposition 7.5. Assume Q is bounded, then a solution to the Dirichlet 
problem is in the class u € C?(Q) uniquely determined up to a constant. 


Proof. Exercise. Hint: Multiply the differential equation Aw = 0 by w and 
integrate the result over 2. 
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Another proof under the weaker assumption u € C1(Q) 9 C?(Q) follows 
from the Hopf boundary point lemma, see Lecture Notes: Linear Elliptic 
Equations of Second Order, for example. 


7.3.3 Mixed boundary value problem 


The Mixed boundary value problem (third boundary value problem) is to 
find a solution u € C?(Q) N Cl(Q) of 


Map = On (7.8) 

a +hu = © on AQ, (7.9) 

where ® and h are given and continuous on 0.e ® and fA are given and 
continuous on OQ. 


Proposition 7.6. Assume Q is bounded and sufficiently regular, then a 
solution to the mized problem is uniquely determined in the class u € C?(Q) 
provided h(a) > 0 on OQ and h(x) > 0 for at least one point x € OQ. 
Proof. Exercise. Hint: Multiply the differential equation Aw = 0 by w and 
integrate the result over 2. 


7.4  Green’s function for A 


Theorem 7.1 says that each harmonic function satisfies 


ua) = (s.2) uly) u(y) he) dS, (7.10) 


Ony 


where 7(y, x) is a fundamental solution. In general, u does not satisfies the 
boundary condition in the above boundary value problems. Since 7 = s+ 4, 
see Section 7.2, where ¢ is an arbitrary harmonic function for each fixed 2, 
we try to find a @ such that u satisfies also the boundary condition. 
Consider the Dirichlet problem, then we look for a ¢ such that 


Vy,t) =0, ye oQ, rEQ. (7.11) 


Then 


u(x) = -{ Only, 2) u(y) dSy, reQ. 
an Ony 
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Suppose that u achieves its boundary values ® of the Dirichlet problem, 
then 


ney | CUD eG) ast (7.12) 


an Ony 


We claim that this function solves the Dirichlet problem (7.4), (7.5). 


A function y(y,x) which satisfies (7.11), and some additional assump- 
tions, is called Green’s function. More precisely, we define a Green function 
as follows. 


Definition. A function G(y,x), y, 7 € 0, x # y, is called Green function 
associated to 2 and to the Dirichlet problem (7.4), (7.5) if for fixed x € Q, 
that is we consider G(y,x) as a function of y, the following properties hold: 


(i) Giy,z) € C?(Q\ {z})NCO\ {z}) A,Gly,z) =0, r#y. 

(ii) Gly, 2) — (|e — yl) € C?(Q)/N CM). 

(iii) G(y,r) =O0if ye OQ, «Fy. 

Remark. We will see in the next section that a Green function exists at 


least for some domains of simple geometry. Concerning the existence of a 
Green function for more general domains see [13]. 


It is an interesting fact that we get from (i)-(iii) of the above definition two 
further important properties. We assume that 2 is bounded, sufficiently 
regular and connected. 


Proposition 7.7. A Green function has the following properties. In the 
case n = 2 we assume diam Q < 1. 


(A) Glw,y) =G(y,2) (symmetry). 
(B) 0<G(a,y)<s(le—yl), @, yEeQ, Fy. 
Proof. (A) Let a, 2@ €Q. Set By = B(x), i = 1, 2. We assume 


B; C Q and By NM By = O. Since G(y,a) and G(y, x) are harmonic in 
Q\ (By U B3) we obtain from Green’s identity, see Figure 7.4 for notations, 
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Figure 7.4: Proof of Proposition 7.7 


0 
Torco (At 9 esa) 
8(2\(BiUB2)) Ny 


The integral over OQ is zero because of property (iii) of a Green function, 
and 


) ) 

Gy, c«?)—Gy,2) -— G(y,¢?)—Gy, 2% as 

[.,, (eee grote) — Geral) z-Glu.al) Jas, 
> Gc, 2), 

[ (Gena G2) = G(y.2)32cty,0%)) as, 
OBz Ony Ony 
2h —G(2), a)) 

as p — 0. This follows by considerations as in the proof of Theorem 7.1. 


(B) Since 
G(y, x) a s(x => yl) + oly, x) 
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and G(y,x) = 0 if y € OQ and « € 2 we have for y € 0D 


o(y, x) = —s(|x — y}). 


From the definition of s(|z — y|) it follows that d(y,z) < 0 if y € ON. 
Thus, since A,@ = 0 in 2, the maximum-minimum principle implies that 
o(y, x) <0 for all y, x € Q. Consequently 


G(y, x) < s(|x =4 yl), Hn Q, x # y- 
It remains to show that 
Gly,z)>0, 2, yeQ, «Fy. 


Fix « € 2 and let B,(a) be a ball such that B,(x) C Q for all 0 < p < po. 
There is a sufficiently small pg > 0 such that for each p, 0 < p< po, 


Giy,z) >0 for ally € B(x), «4 y, 


see property (iii) of a Green function. Since 


A,Giy,z) = 0 inQ\ B,(z) 
0 if y € OB,(x) 
0 if ye 0a 


pee 

Ss 

aa 
Iv 


it follows from the maximum-minimum principle that 


Giy,x) >0 on 2 \ B,(2). 


7.4.1 Green’s function for a ball 


If 2 = Br(0) is a ball, then Green’s function is explicitly known. 

Let 2 = Bpr(0) be a ball in R” with radius R and the center at the 
origin. Let x, y € Br(O0) and let y’ the reflected point of y on the sphere 
OBR(0), that is, in particular |y||y’| = R?, see Figure 7.5 for notations. Set 


G(e,y) = s(r)— 8 (Fn), 


where s is the singularity function of Section 7.1, r = |x — y| and 


n n R2 2 
f= nab (m- Fu) - 
i=1 


i=1 
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Figure 7.5: Reflection on 0Br(0) 


This function G(x, y) satisfies (i)-(iii) of the definition of a Green function. 
We claim that 3 
u(x) = -{ ~—— G(x, y)® dS, 
ABR(0) ONy ‘ 


is a solution of the Dirichlet problem (7.4), (7.5). This formula is also true 
for a large class of domains 2 C R”, see [13]. 


Lemma. . al? 
0 1 R°-|az 

=e He) as) ly—2|" 
lyl=R ae A 


Ony 


Proof. Exercise. 


Set 


(7.13) 
which is called Poisson’s kernel. 


Theorem 7.2. Assume ® € C(O0Q). Then 
ule)=f  H(a.9)®lu) a5, 
OBR(0) 


is the solution of the first boundary value problem (7.4), (7.5) in the class 
C2(Q)NC(Q). 
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Proof. The proof follows from following properties of H(z, y): 
(i) H(x,y)eC™, lyl=R, lal <R, Fy, 

(ii) AcH(x,y)=0, lal <R, ly =R, 

(iii) JoBr@) Hie gy wdSyats ale Rk, 

(iv) H(az,y)>0, lyl=R, |a| < R, 


(v) Fix ¢ € OBr(0) and 6 > 0, then lim,_.¢ jz;<r A(x, y) = 0 
uniformly in y € OBr(0), |y—¢| > 6. 


(i), (iv) and (v) follow from the definition (7.13) of H and (ii) from (7.13) 
or from 


y — 0G») 


’ 
Ony 


ycOBrR(0) 


G harmonic and G(z, y) = G(y, 2). 
Property (iii) is a consequence of formula 


u(x) = i A(x, y)u(y) dSy, 


for each harmonic function u, see calculations to the representation formula 
above. We obtain (ii) if we set u = 1. 


It remains to show that u, given by Poisson’s formula, is in C(Br(0)) and 
that u achieves the prescribed boundary values. Fix ¢ € O0Br(0) and let 
x € Br(0). Then 


u(x) -8(Q) = [ acy 8) BW) ~ BO) 45 
= a 


where 


ee i H(«,y) (®(y) — ®(0) aS, 
OBR(0), |y—¢|<d 


Il = (x,y) (®(y) — ®(¢)) dSy. 
OBR(0), |y—<¢|26 
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For given (small) e > 0 there is a 6 = d(€) > 0 such that 
|O(y) — O(¢)| <e€ 


for all y € OBR(0) with |y— ¢| < 6. It follows || < € because of (iii) and 
(iv). 
Set M = maxgz,o) |¢|. From (v) we conclude that there is a 6’ > 0 such 


that 
€ 


2Mw,Rr-1 
if x and y satisfy |x —¢| < 6’, |y—¢| > 6, see Figure 7.6 for notations. Thus 


H (x,y) < 


Figure 7.6: Proof of Theorem 7.2 


|I2| < € and the inequality 


|u(a) — ®(¢)| < 2€ 


for « € Br(0) such that |x — ¢| < 6’ is shown. 


Remark. Define 6 € [0,7] through cosé = x- y/(|z||y|), then we write 
Poisson’s formula of Theorem 7.2 as 


2 2 
= 1 
u(x) = a | O(y 5 dS,,. 
wrR  Japa(o) ~ (|x|? +R? — 2|x|Rcosd)"/ 
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In the case n = 2 we can expand this integral in a power series with respect 
to p:= |x|/R if |z| < R, since 


R- |x|? _ 1p" 
|z| + R? — 2|2|Rcos6 p? — 2pcosd +1 
co 
= 142 s p” cos(nd), 
n=1 


see [16], pp. 18 for an easy proof of this formula, or [4], Vol. II, p. 246. 


7.4.2  Green’s function and conformal mapping 


For two-dimensional domains there is a beautiful connection between confor- 
mal mapping and Green’s function. Let w = f(z) be a conformal mapping 
from a sufficiently regular connected domain in R? onto the interior of the 
unit circle, see Figure 7.7. Then the Green function of Q is, see for exam- 


@) 
w=f(z) (w) 


—_—_—_—_=> 


Figure 7.7: Conformal mapping 


ple [16] or other text books about the theory of functions of one complex 
variable, 


1 


pe f(z) f Zo) 
21 


f(z) — f(0) 


b] 


G(z, 2) = 


where z = 41 +122, 20 = y1 + tYyo- 


7.5 Inhomogeneous equation 


Here we consider solutions u € C?(Q) 9 C(Q) of 
—Au = f(z) nQ? (7.14) 
0 on dQ, (7.15) 


U 
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where f is given. 

We need the following lemma concerning volume potentials. We assume 
that Q is bounded and sufficiently regular such that all the following integrals 
exist. See [6] for generalizations concerning these assumptions. 

Let for x € R”, n> 83, 


1 
V(z) =| LO) 2 dy 


and set in the two-dimensional case 


va) = f fin (—) dy. 


We recall that w,, = |OB,(0)]. 


Lemma. 


(i) Assume f € C(Q). Then V € C1(R") and 


) 1 
Ve(t) = [ tox (=) dy, ifn 2 3, 


vale) =f sodge,(m(eaG)) a Yam? 


(ii) If f ¢ C1(Q), then V € C?(Q) and 


> 
< 
| 


—(n—2)un f(x), cE, n> 3 
—20f (2), CEO, w=2: 


> 
~ 
| 


Proof. To simplify the presentation, we consider the case n = 3. 

(i) The first assertion follows since we can change differentiation with inte- 
gration since the differentiate integrand is weakly singular, see an exercise. 
(ii) We will differentiate at x € Q. Let B, be a fixed ball such that x € By, 
p sufficiently small such that B, C Q. Then, according to (i) and since we 


have the identity 
St ee 
Ox; \|x — y| dys \|a— y| 
which implies that 


10m (aa) * an PO Raa) tO 
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we obtain 


vate) = ff val i 

Tue? ae (=a) +f, an (Ea) 
= foo, 1 = (gaa) 

+f ca (fwjgaq) thelemy) 


p 


- Pe i — dy 


1 
ee Sy:y eerie vf, A mers dSy, 


where n is the exterior unit normal at OB,. It follows that the first and 
second integral is in C!(Q). The second integral is also in C!(Q) according 
to (i) and since f € C1(Q) by assumption. 

Because of A,.(|z — y|~') =0, « 4 y, it follows 


=f Sang a (Ea) @ 


Bp j= 


a) 1 
_- ae 3 Oni (i = i) ie 


Now we choose for B, a ball with the center at x, then 


AV=H]+h, 


where 


— esd 
4 i, p(x) 5 Ss wl = 


— -f Ty) dy. 
aB,(a) ns ’ 


We recall that n-(y— x) = p if y € OB, (a). It is I) = O(p) as p > 0 and 
for Iz we obtain from the mean value theorem of the integral calculus that 
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for a 7 € OB, (x) 


1 
Ip = ——xf y if dS. 
p? ( ) OB,(a) 
= —wnf (¥), 
which implies that lim,—9 2 = —wn f(a). 


In the following we assume that Green’s function exists for the domain Q, 
which is the case if 2 is a ball. 


Theorem 7.3. Assume f € C1(Q)M C(Q). Then 


wa) =f Glens ay 
is the solution of the inhomogeneous problem (7.14), (7.15). 


Proof. For simplicity of the presentation let n = 3. We will show that 
ula) = f Gle.nflw) dy 


is a solution of (7.4), (7.5). Since 


af 


G(a,y) = tale + $(z,y), 


where ¢ is a potential function with respect to x or y, we obtain from the 
above lemma that 


1 1 
Ak siz: wef foe tut fA oole atten a 


Ar |x — 


= =F), 


where x € 2. It remains to show that u achieves its boundary values. That 
is, for fixed xp € OQ we will prove that 


lim u(x) =0. 
r—x0, cEQ 


Set 
u(x) =h+lo, 
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where 


h(x) = | Cayrayay 
Q\Bp(xo) 

Ley) = : CGaiG) di. 
QNBp (xo) 


Let M = maxg|f(a)|. Since 


1 
G = ——_ 
(x,y) ed + O(2,y), 
we obtain, if x € B,(xo) NQ, 
M dy 
Jbl < mo tO) 

Ar QNBp(xo) |x ~~ y| 
M d 

ae + O(¢?) 
An Bop(x) |x _ yl 


= O(p") 

as p > 0. Consequently for given € there is a po = po(e) > 0 such that 
|I2| < : for all 0< p< po. 

For each fixed p, 0 < p < po, we have 


lm (x) =0 


r—x9, LEQ 


since G(x, y) = 0 if ye 2 \ B,(xo) and G(a, y) is uniformly continuous in 
z € Byjo(xo) NQ and y € O \ B,(xo), see Figure 7.8. 


Remark. For the proof of (ii) in the above lemma it is sufficient to assume 
that f is Hélder continuous. More precisely, let f € C*(Q), 0 < A < 1, then 
V € C?A(Q), see for example [9]. 
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7.6 


Figure 7.8: Proof of Theorem 7.3 


Exercises 


. Let y(x,y) be a fundamental solution to A, y € Q. Show that 


- f ve) A®(x) dx = ®(y) for all ®E C2(Q) . 


Hint: See the proof of the representation formula. 


. Show that |z|~! sin(k|z]) is a solution of the Helmholtz equation 


Au+t k?u = 0 in R” \ {0}. 


. Assume u € C?(Q), Q bounded and sufficiently regular, is a solution 


of 


w in Q 
0 on ON. 


Au 


U 


Show that u=0 in 2. 


. Let Ng = {z € R?: x1 >0, 0 < 22 < x tanas, O0<a<7. Show 


that 
u(x) = re" sin 2) 


is a harmonic function in Qg satisfying u = 0 on OQq,, provided k is 
an integer. Here (r,@) are polar coordinates with the center at (0,0). 
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5. Let u € C?(Q) be a solution of Au = 0 on the quadrangle Q = 


10. 


(0,1) x (0,1) satisfying the boundary conditions u(0, y) = u(1,y) =0 
for all y € [0,1] and uy(xz,0) = uy(x,1) = 0 for all x € [0,1]. Prove 


that u=0 in 2. 


. Let u € C?(R”) be a solution of Au = 0 in R” satisfying u € L?(R"), 


1 ine. te) de <0. 
Show that u = 0 in R”. 


Hint: Prove 


const. 
i \Vul? dx < Re i |u|? dx, 
Br(0) Bor(0) 


where c is a constant independent of R. 

To show this inequality, multiply the differential equation by ¢ := 7?u, 
where 7 € C! is a cut-off function with properties: 7 = 1 in Br(0), 
n = 0 in the exterior of Bor(0), 0 < 7 < 1, |Vn| < C/R. Integrate 
the product, apply integration by parts and use the formula 2ab < 
ea? + 132, e>0. 


. Show that a bounded harmonic function defined on R” must be a 


constant (a theorem of Liouville). 


. Assume u € C?(Bi(0)) A C(Bi(0) \ {(1,0)}) is a solution of 


Au = 0 in B,(0) 
wu = O on OB,(0) \ {(1,0)}. 


Show that there are at least two solutions. 
Hint: Consider 
1 — (x? + y?) 


MODS aaa 


. Assume 2 C R” is bounded and u, v € C?(Q)NC(Q) satisfy Au = Av 


and maxago |u — v| < € for given « > 0. Show that maxg|u — v| < e. 


Set Q = R” \ B,(0) and let wu € C?(Q) be a harmonic function in Q 
satisfying lim),),.. u(x) = 0. Prove that 


max |u| = max |u| . 
2 0a 


Hint: Apply the maximum principle to NM Br(0), R large. 
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11. 


12. 


13. 


14. 


Let Ng = {2 € R?: 2, > 0, 0 < ag < a tanas, O< a <7, 
Qe. = QM Br(0), and assume f is given and bounded on Qa, p. 
Show that for each solution u € C1(Qa,R) A C?(Qa,r) of Au = f in 
Qo R satisfying u= 0 on OQa,.R2M Br(0), holds: 


For given € > 0 there is a constant C(e) such that 
ju(x)| <C(6 |ale7 in Qa,p. 


Hint: (a) Comparison principle (a consequence from the maximum 
principle): Assume Q is bounded, u,v € C?(Q) A C(Q) satisfying 
—Au < —Av inQ and u< von OQ. Then u < vin 2. 


(b) An appropriate comparison function is 
v = Ara‘ sin(B(O+7)) , 
A, B, 7 appropriate constants, B, 7 positive. 
Let Q be the quadrangle (—1,1) x (—1,1) and u € C?(Q) N C(Q) a 


solution of the boundary value problem —Au = 1 in Q, u = 0 on O02. 
Find a lower and an upper bound for u(0, 0). 


Hint: Consider the comparison function v = A(a? + y”), A = const. 
Let u € C?(B,(0)) M C(Ba(0)) satisfying u > 0, Au = 0 in B,(0). 
Prove (Harnack’s inequality): 
a’ *(a— Cl) a’*(a+ ICI) 
(a+ |¢])"? o=|C)0= 
Hint: Use the formula (see Theorem 7.2) 


Wy) a2 lef 8 as 
(uv) i} ? 


QW, x|=a |x = y|” 


u(0) <ul) < u(0) . 


for y= ¢ and y = 0. 


Let ¢(0) be a 27-periodic C*-function with the Fourier series 


(0) = S— (an cos(nO) + by sin(n8)) . 
n=0 
Show that 
i= S- (an cos(n@) + bp sin(né)) r” 
n=0 


solves the Dirichlet problem in B,(0). 
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15. 


16. 


17. 


18. 


19. 
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Assume u € C?(Q) satisfies Au = 0 in Q. Let By(¢) be a ball such 
that its closure is in Q. Show that 


[Deu(¢)| <M (<2) i 


a 


where M = supzep,(c) |U(x)| and Yn = 2nWp—1/((n — 1)wn). 


Hint: Use the formula of Theorem 7.2, successively to the k th deriva- 
tives in balls with radius a(|a| — k)/m, k =0,1,...,m—1. 


Use the result of the previous exercise to show that u € C?(Q) satis- 
fying Au = 0 in 2 is real analytic in 2. 


Hint: Use Stirling’s formula 


nl =n"e” (van +O (=) 


as n — oo, to show that wu is in the class Cx,,(¢), where K = cM and 
r =a/(eyn). The constant c is the constant in the estimate n” < ce”n! 
which follows from Stirling’s formula. See Section 3.5 for the definition 
of a real analytic function. 


Assume Q is connected and u € C?(Q) is a solution of Au = 0 in Q. 
Prove that u = 0 in 2 if D°u(¢) = 0 for all a, for a point ¢ € 2. In 
particular, uw = 0 in 2 if u=0 in an open subset of 2. 


Let © = {(21,22,23) € R?: a3 > O}, which is a half-space of R°. 


Show that 
1 1 


dnje— yl Snje— Bl 


G(x, y) = 
where Y = (Y1, y2, —y3), is the Green function to Q. 


Let 2 = {(x1, 22,73) € R®: x?+23+23 < R?, x3 > 0}, which is half 
of a ball in R®. Show that 


1 R 
4njx—y|  4rly||a — y*| 
1 | R 


4n|c—y| | 4rlylla— 9’ 


G(x, y) 


where 9 = (41,92, —ys), y* = R’y/(lyl?) and y* = R?y/(ly|*), is the 
Green function to 2. 
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20. 


21. 


22. 


23. 


24. 


25. 


Let Q = {(21, 22,273) € R°?: x2 >0, x3 > 0}, which is a wedge in R?. 
Show that 
1 1 
ie—y ae 
ie a 
iae—yl * Sea 


G(x, y) 


where ¥ = (y1,y2,—ys), ¥’ = (y1,—y2,y3) and 7 = (y1,—ye, —ys), is 
the Green function to 2. 


Find Green’s function for the exterior of a disk, i. e., of the domain 
Q={reER?: |z| > R}. 


Find Green’s function for the angle domain Q = {z EC: 0< argz< 
at},0<a<rT. 


Find Green’s function for the slit domain Q = {z € C: 0<argz< 
27}. 


Let for a sufficiently regular domain 2 € R”, a ball or a quadrangle 
for example, 


F(a) = | Kau) dy, 


where K(x, y) is continuous in Q x Q where x 4 y, and which satisfies 
c 


IK <. 
Pus je yl? 
with a constants c and a, a <n. 


Show that F(x) is continuous on 2. 


Prove (i) of the lemma of Section 7.5. 


Hint: Consider the case n > 3. Fix a function 7 € C'(R) satisfying 
O0<n<10<7' < 2, n(t) = 0 fort < 1, n(t) = 1 fort > 2 and 
consider for € > 0 the regularized integral 


ve) = f sno 


=o 


where 7. = (|x — y|/e). Show that V. converges uniformly to V 
on compact subsets of R” as « — 0, and that OV.(x)/0x; converges 
uniformly on compact subsets of R” to 
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ase— 0. 


26. Consider the inhomogeneous Dirichlet problem —Au = f in Q, u = 
¢ on OQ. Transform this problem into a Dirichlet problem for the 
Laplace equation. 


Hint: Set u = w+ v, where w(x) = fo s(|x — yl) f(y) dy. 
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